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INTRODUCTION 

During  the  last  decade  considerable  progress  has  been  made  in 
understanding  the  self-diffusion  of  metal  atoms  on  well  defined  crystal 
planes  [1],  This  advance  has  occurred  through  the  use  of  the  field  ion 
microscope,  which  makes  it  possible  to  visualize  individual  atoms,  and  to 
prepare  highly  perfect  surfaces.  With  this  microscope  it  is  possible  to 
determine  by  direct  observation  the  mean  square  displacement  of  atoms  ad- 


sorbed at  a surface,  and  thus  to  derive  their  diffusion  coefficient  D 
from  Einstein's  relation,  which  for  one-dimensional  motion  assumes  the 
form  [2] 

<(AR)2)  = 2Dt  . (1) 

A considerable  fund  of  information  about  the  directionality  of  surface 
diffusion,  about  the  role  of  surface  structure,  and  about  the  energetics 
of  atomic  jumps  has  been  built  up  this  way.  Despite  that,  very  little  is 
known  about  the  actual  atomic  motions  involved  in  the  diffusion  process. 

It  is  by  now  customary  to  assume  that  in  moving  over  a surface 
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where  I^Cx)  is  the  modified  Bessel  function  of  order  n.  The  mean  square 
displacement  \ (Ar)“>  is  then 

<(AR)2>  * 2at22  . (3) 

The  validity  of  this  model  cannot  be  established  from  an  examination  of  the 
mean  square  displacement  alone.  This  requires  a detailed  comparison  of  the 
probability  density,  determined  from  experiments,  with  the  values  predicted 
by  Eq.  (2).  One  such  comparison  has  been  made  [4],  for  the  behavior  of 
Re  atoms  on  the  (211)  plane  of  tungsten.  Within  the  sizable  statistical 
error,  arising  from  the  limited  data  set,  it  appears  that  Eq.  (2)  is 
adequate  to  describe  the  actual  diffusion  process. 

It  still  remains  to  be  established,  however,  that  surface  diffusion 
involves  only  jumps  between  nearest  neighbors.  With  the  data  available  right 
now,  it  is  not  possible  to  exclude  the  possibility  that  occasionally,  an 
atom  might  carry  out  a jump  spanning  two  nearest  neighbor  spacings  instead 
of  just  one.  The  aim  of  this  thesis  is  to  examine  the  feasibility  of 
detecting  such  double  jumps  in  diffusion  measurements.  To  this  end  we 
idealize  the  one-dimensional  diffusion  of  an  adatom  as  involving  single  as 
well  as  double  jumps  at  random  moments,  the  former  occurring  at  a rate  la, 
the  latter  at  a rate  29.  The  statistical  consequences  of  this  model  are 
developed  in  Chapter  I;  Chapter  II, is  devoted  to  a detailed  discussion  of 
the  conditions  that  must  be  satisfied  in  an  experiment  designed  to  measure 
the  contributions  of  multiple  jumps.  Especially  important  in  these  con- 
siderations are  estimates  of  the  statistical  errors  involved  in  such 
determinations.  These  estimates  rely  on  a fairly  complicated  mathematical 
apparatus;  this  is  developed,  at  some  length,  in  separate  appendices. 
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I.  RELATION  BETWEEN  PROBABILITY  DENSITIES,  JUMP  RATES,  AND  MOMENTS 


Oddly  enough,  the  behavior  of  random  walks  involving  jumps  of 


different  lengths  has  been  neglected  in  the  literature.  Lakatos -Linde nberg 


and  Schuler  [5]  have  examined  the  statistics  of  a particle  making  transitions 


to  both  nearest  and  next -nearest  neighbor  sites;  however,  they  considered 


only  jumps  at  fixed  time  intervals,  a model  which  is  not  strictly  appropriate 


to  actual  diffusion  processes.  We  will  therefore  develop  here  the  conse 


quences  of  extending  the  usual  model  for  one-dimensional  diffusion,  in  which 


atoms  are  presumed  to  execute  jumps  after  random  time  intervals,  by  allowing 


displacements  both  between  nearest  and  next-nearest  neighbors.  Transitions 


in  a given  direction  are  presumed  to  occur  at  a rate  a to  the  nearest 


Motion  to  the 


right  and  to  the  left  is  equally  probable 


A.  Probability  Densities  and  Jump  Rates 


The  first  step  in  developing  our  model  is  to  find  the  probability 


designate  this  probability  as  p (t)  and  find  it  via  the  Kolmogorov 


equation,  which  for  this  system  assumes  the  form  [6] 


here  a is  the  rate  at  which  adatoms  execute  unit  displacements  and  3 is 


the  rate  at  which  adatoms  execute  double  displacements.  The  distance  l 


between  adjacent  sites  is  taken  to  be  unity 
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In  order  to  solve  the  Kolmogorov  equation  and  thereby  find  the 
probability  density  p^(t),  we  need  to  uncouple  the  time  variable  from  the 
position  variables.  This  is  accomplished  by  the  use  of  the  well  known 
generating  function  for  the  probability  density.  The  generating  function 
G(t,z)  is  defined  as 

G(t,z)  = E znp™(t)  • (5) 

n 

The  symbol  £ is  throughout  used  to  designate  summation  over  the  limits 
n 

n = -co  to  «,  that  is 

CO 

E = E . 

n n=-ao 


We  use  Eqs.  (4)  and  (5)  lo  derive  a simple  differential  equation  involving 
G(t,z).  Once  the  generating  function  has  been  obtained  from  this,  p™(t) 
can  be  found  using  the  definition  in  Eq.  (5). 

To  derive  the  differential  equation  for  G(t,z),  we  note  that  z 
is  time  independent;  it  therefore  follows  from  Eq.  (5)  that 


ac  , „dp>> 

a?  (t'z)  ‘ E 5 — 

n 


(6) 


Multiplying  the  Kolmogorov  equation  by  z and  summing  over  n,  we  find  that 


dG(t,z)  „ nr„  m m „ , a,  m m _ m i 

~H  ' Z Z [Sp„-2+“I’n-l-2<,1'fB)f,n+a,Pn+l+Sp„+2) 
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From  these  five  equations  above,  we  conclude  that 


n z 


(9) 


This  can  be  written  more  simply  as 


dOfti.2).  = Z znpin{g(z2  + z'2)+a(2  + z'I)-2(&  + S)]  . 

ot  n 

n 


(10) 


Recalling  the  definition  of  G(t,z)  in  Eq.  (5  ),  Eq. (10)  can  also  be  written 


as 


= G(t,z)^(z2  + z'2)  +a(z  + z‘1)-2(a  + 9)} 


(ID 


It  is  apparent  that  the  differential  equation  for  the  generating  function 
is  much  simpler  than  the  original  Kolmogorov  relation.  Equation  (11)  has 
the  simple  solution, 

G(t,z)  = C exp[  [3  (z2  + z 2)+a(z  + z i)-2  (ct  +g  ) ] t}  . 

The  constant  C can  be  found  from  the  boundary  conditions.  Given 
that  the  adatom  begins  its  random  journey  at  m,  we  know  that 


m , 


p (t  =0)  = 5 
rn  nm 


where  6 

rim 

is  Kronecker's  delta.  Therefore 

\ 

& 

j 

G(t=0,z)  = Z 5^  = zm 

n 

j 

( . | I 

and 

1 

j 

! i 

G(t,z) 

= zmexp{  [3  (z2  + z 2)  + or(z  -t-z  ^)-2  (ot  +3)]t]  . 
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In  order  to  determine  the  probability  density  p‘  (t),  we  must 

n 

expand  G(t,z)  in  a Laurent  series.  This  can  be  done  using  the  Schldmilch 


series  [7] 


exp  [ 2 (y  + y = £ ynin(x)  • 

n 


This  allows  us  to  express  G(t,z)  as  a series  involving  modified  Bessel 
functions  I of  order  n.  Equation  (12)  can  be  written  in  the  equivalent 


G(t,z)  = zmexp{£t  [z2  + z 2]}  exp  {art  [z  + z l]}  exp{-2  (a  + 3 )t} 


According  to  Eq.  (13);  however, 


exp{0t[z  + z'2]}  = 2 1 (23t)z  1, 

nl  1 

f -1  , n2 
explat[z+z  ]s  =2  1 (2at)z 

n2  n2 

From  the  definition  of  the  generating  function,  in  Eq.  (5),  it 


follows  that 


_ . _ m n 

G(t,z)  = 2 pn  z 

n 


(2n  +n  +m) 

= exp[-2(a+P)t]  2 I (2pt)I  (2at)z  z 


n ' n„ 
n^^  1 2 


Identifying  n as 


we  have 


n = 2n^  + n^  + m , 


n^  = (n-m)  - 2n^. 


From  Eq.  (15)  the  probability  density  p^(t)  can  now  be  written  as 

£(t>  - -'2<“  + S)tpk(2St)I(n.„).2t(20,t)  , 


n I 

Li  1 


.*f  / * * « t '■‘■W**  V-V-  * r 
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where  k has  been  substituted  for  n^,  the  dummy  variable.  When  only  unit 
displacements  occur,  3=0  and 


p“(t)  « e'2at  E I (0) 
n k k 


I(n-m)-2k('2at) 


m,..  a~20tt-  T 

Pn(t)  = 6 1 (n-m) 


Note  that  1^(0)  ^ 0 if  and  only  if  k = 0;  for  that  case  IQ(0)  = 1,  so  that 

(2«t)  (17) 

for  3=0.  This  agrees  with  the  standard  expression  for  an  adatom  dif- 
fusing via  unit  displacements.  When  only  double  displacements  are  allowed, 

P®(t)  - e‘20t  E L (23 1)  I, 


(n-m) -2k 


(0) 


= e”2^CI  (2^t) 

i(n-m)/2U  C> 


(18) 


for  a = 0;  here  n-m  must  be  an  even  integer. 


B.  Jump  Rates  and  Moments 

In  order  to  determine  the  rates  3 and  or  at  which  a diffusing 
adatom  makes  double  and  single  displacements,  we  must  reduce  the  above 
equations  into  explicit  expressions  for  observable  quantities.  Such 
quantities  are  the  moments  u/  of  the  displacements  executed  by  the  adatom 
in  the  one  dimensional  walk. 

Suppose  that  an  adatom  is  initially  at  position  m^.  We  allow  t 
seconds  for  the  adatom  to  perform  a random  walk  and  call  its  position  at 
the  end  of  the  walk  n^.  The  i-th  initial  and  final  positions  that  we 
observe  will  be  designated  by  nu  and  n^.  On  M occations  m.  and  n^  are 
recorded;  we  are  thus  able  to  define  M displacements  x^,  where 

(n^  - nt)  . 


X.  = 
i 


(19) 
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If  we  were  to  make  an  arbitrarily  large  number  of  observations  we  would  be 
able  to  find  p/,  the  r-th  moment  about  zero,  which  is  defined  as 


M. 


I 

r 


lira 

M-*° 


(20a) 


As  shown  in  Appendix  B,  the  r-th  moment  about  the  mean 


M- 
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lim  _1 
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i=l 


(*• 


“i> 


r 


(20b) 


is  often  a more  useful  quantity  than  p> ' . The  one  dimensional  random  walk 
is  synsnetric  about  the  adatom's  initial  position,  so  the  odd  moments  about 
zero  vanish  in  the  limit  of  an  infinite  number  of  observations.  Because 
the  mean  is  an  odd  moment  and  therefore  zero,  we  know  from  Eqs.  (20a) 
and  (20b)  that  for  our  example 


r 


(21) 


In  order  to  find  these  moments  in  terms  of  tne  jump  rates  a and 
3,  we  must  seek  a way  of  relating  the  moments  to  the  generating  functions 
given  by  Eqs.  (5)  and  (12),  that  is  to 

G(t , z ) = Z P™(t)zn 
n 

= zra  expt[3(z2+z‘2)+a(z  + z'1)-2(a+3)]t}  . 


From  Eq.  (20a)  we  know  that  the  mean  displacement  for  the  adatom  per- 
forming the  random  walk  is 


lim 

M-ta> 


A S 

M i-i 


(22) 
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We  also  know  that  the  probability  for  the  diffusing  adatom  to  occupy  site  n, 
having  started  at  site  m t seconds  earlier,  is  p^(t).  An  equivalent 
expression  for  the  mean  position  of  the  adatom  is  then 

- E n p“(t)  . (23) 

n 

The  similarity  between  the  relation  for  the  generating  function  and  for  the 
mean  causes  us  to  realize  that 


G(t,z)  = £ p®(t)nzn_1 


/ 3 . „ . _ mn 

(z  ^7)  G(t,z)  =*  Z n pn  z 

n 


(z  ^)  G(t,z)  I = Z n p”; 

| z=l  n 


Without  loss  of  generality,  we  define  the  initial  position  m of 
the  adatom  as  zero,  so  that  the  generating  function  in  Eq.  (12)  assumes  a 
simpler  form: 


G(t,z)  = exp{[3(z2+z  2)+a(z  + z 1)-2  (a +3  ) ]tl 


By  substituting  Eq.  (25)  into  Eq(24)  we  can  write  the  mean  as 
M>  ^ = z(^)  exp([0(z2  + z 2)+a(z  + z 1)-2  (a  + 3)  ]t) 

z*l 

= z[[3(2z-2z  2)+a(l-z  2)]t}  exp{[3(z2  + z 2)+a(z  + z ^)-2(a  + 3)]t} 


l([3(2-2)  + a(l-l)]t}  «xp[[p(l  + l)H»(l  + l)  *2  (a  + 0 ) ] t}  = 0 . 


Just  as  expected  from  symmetry,  ^ vanishes. 


Any  r-th  moment  can  be  found  by  the  same  route  that  lead  to 


Eq.  (24).  Notice  that 


. 3 .r  _ n m 

<z  1 z p„ 
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Z 
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r n m 
n z p 
*n 


(27) 


Evaluating  Eq.  (27)  at  z = 1 leads  to  the  well  known  relation 


(z  ^)r  G(t,z) 


_ r m 

, " s n pn 

z=l  n 


or 


M. ' = (z  Tl)r  G(t , z ) 


(28) 


(29) 


z=l 


We  now  have  accessible  any  r-th  moment  of  the  adatom  displacements  occurring 
during  a random  walk  lasting  t seconds.  The  derivation  of  moments  via 
Eq.  (29)  is  conceptually  simple;  the  actual  calculations  of  higher  moments 
are  tedious  and  are  therefore  reserved  for  Appendix  A. 

Knowing  the  moments  of  the  adatom  displacements  as  a function  of 
the  jump  rates  a and  P,  we  are  ready  for  the  problem  of  finding  a and  $ in 
terms  of  the  moments.  In  Appendix  A,  Eqs.  (25)  and  (29)  are  used  to  derive 
the  first  two  non-zero  moments 


u.'2  « 2(43  + a)t  (30a) 

^ = 7^2  + 30i£)2  - 12(23  + a)  t.  (30b) 

, . i 

For  the  case  where  no  double  jumps  occur  and  so  3=0,  the  second  moment  is 
M>2  - 2a t , (31) 


just  as  we  would  suspect  from  Eq.  (3)  for  1=1.  Equations  (30a)  and  (30b) 
are  linearly  independent;  hence  a and  3 can  be  expressed,  in  terms  of  p.' 


so  that  at  and  3t  can  be  found  in 


terras  of  moments  about  the  mean  as  well  as  about  zero 


C.  Moments  and  Probability  Densities 


From  Eq.  (16)  it  is  obvious  that  the  probability  density  p (t) 


for  a random  walk  involving  double  and  single  jumps  is  uniquely  characterized 


by  the  two  jump  rates  3 and  a,  and  by  the  diffusion  time  t.  We  can  reduce 


the  number  of  parameters  necessary  to  specify  the  distribution  by  repre 


seating  p (t)  as  a function  of  the  second  moment  '40  and  the  ratio  P/a 


Equation  (30a)  gives  the  second  moment  in  terms  of  a,  P,  and  t,  as 


(1  + 4P/a) 


Also 


Substituting  Eqs.  (34a),  (34b),  and  (34c)  into  Eq.  (16)  for  the  probability 


density  p (t)  yields  this  quantity  as  a function  of  u and  P/a 
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In  Appendix  D is  listed  the  FORTRAN  program  PPLOT  [8]  which 
calculates  and  plots  the  probability  p°  as  a function  of  lattice  position 
n for  various  values  of  and  3 /a.  In  Fig.  1 we  show  the  behavior  of  an 
adatom  executing  a one-dimensional  random  walk.  The  distribution  of 
distances  for  an  adatom  only  allowed  to  make  single  jumps  is  compared  with 
the  distribution  in  which  the  adatom  can  make  both  double  and  single  jumps 
as  well  as  that  in  which  only  double  jumps  are  allowed.  As  we  vary  the 
second  moment  from  10.  to  1.,  the  distribution  where  only  double  jumps  are 
allowed  begins  to  resemble  that  for  single  jumps  more  closely;  the  mixed 
distribution,  for  which  j3/a  = l.,  seems  to  become  less  like  the  distribution 
for  single  jumps.  This  is  the  sort  of  trend  that  we  are  looking  for.  In 
order  to  do  an  experiment  aimed  at  detecting  double  jumps  there  must  be 
differences  between  the  distributions  for  different  jump  rates  3/a.  The 
trends  in  Fig.  1 suggest  that  the  diffusion  experiment  should  be  done 
under  conditions  such  that  the  second  moment  of  the  displacement  distri- 
bution is  much  less  than  10.0.  Figure  1 also  indicates  that  the  second 
moment  should  be  larger  than  .1;  at  this  value  all  three  distributions 
look  very  much  the  same,  as  little  diffusion  has  occurred  for  any  of  the 
three  examples. 

In  Fig.  2 we  see  how  the  probability  distribution  p°  calculated 
by  PPLOT  compares  with  available  data.  The  experimental  distance  distri- 
bution is  taken  from  Stolt's  observations  on  one -dimensional  diffusion  of 
Re  atoms  over  W(211)  [4],  These  observations  were  not  made  with  the 
intention  of  looking  for  double  jumps,  so  it  is  not  surprising  that  the  data 
are  inadequate  to  show  whether  or  not  adatoms  make  occasional  double  jumps. 
The  experiment  that  can  answer  this  question  will  have  to  be  planned 
specifically  for  the  purpose  of  detecting  such  double  jumps. 


II.  PROBLEMS  IN  THE  DETERMINATION  OF  JUMP  RATES 


I I 


A.  Estimating  Jump  Rates  from  Diffusion  Data 

The  problem  now  is  to  determine,  by  experiment,  the  values  of 
the  jump  rates  a and  j3  for  single  and  double  jumps.  From  Eqs.  (32a), 
(32b),  and  (21),  we  know  that 

at  - £ t^2  + 3(n2)2  - n4]  I 

= i4  " 3^2)2  ' ^2]  * ' 

The  moments  entering  here  are  defined  by  Eq.  (20b)  as 


(36a) 


(36b) 


lim  1 


M = <Xi  • “i>r  * 
1=1 


where  x^  is  the  i-th  observed  adatom  displacement  and  p.|  is  the  mean  of  the 
displacements. 

The  moments  p.^  and  p^  in  Eqs.  (36)  are  quantities  established 
from  an  infinite  number  of  observations.  Such  an  infinitely  large  set  of 
observations  is  called  a population  or  parent  population;  quantities  such 
as  p.^  and  jj.  that  are  derived  from  a population  are  referred  to  as  popula- 
tion values.  Specifically,  p.^  and  p^  belong  to  a class  of  population 
values  called  parent  moments.  In  order  for  p.^  and  p.  to  be  of  any  use  to 
us,  we  must  learn  how  to  estimate  them  from  a finite  number  of  observations. 
Given  a finite  sample  of  observations  taken  from  the  parent  population,  we 
can  generate  statistics.  A statistic  is  defined  as  any  parameter  generated 
by  a finite  number  of  observations;  for  each  population  parameter  there  is 
usually  a corresponding  statistic.  The  statistic  that  corresponds  to  a 


parent  moment  is  called  a sample  moment.  The  r-th  sample  moment  about 
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zero,  m' , is  defined  as 
r 

1 M 

mr  ' M Z (Xi)C  * (38) 

Ui=l 

The  r-th  sample  moment  about  the  mean  is  defined  as 

*r  * 5 <V»i>C  • <39> 

1=1 

The  statistic  that  corresponds  to  the  second  parent  moment  about 
the  mean,  is  the  sample  moment  about  the  mean,  m_.  From  the  definitions 

of  ^ and  m^,  we  see  that  m^  is  a crude  estimate  of  p . If  m^  were  an 
unbiased  estimate  of  > the  average  or  expectation  value  of  would  be 
equal  to  p.^.  This,  however,  is  not  the  case.  In  Appendix  B we  find  that 
the  expectation  value  of  m^,  <m2>  *-s  actually 

<m2>  " ^ ^2  ‘ (40) 

A better  estimate  of  is  then 

' & "2  • <41> 

By  definition,  p^  is  the  unbiased  estimator  of  p^*  In  general,  the  statistic 

A 

6 is  the  unbiased  estimator  of  the  population  value  6 if  and  only  if  [9] 

<5  > = 6 . (42) 


Not  all  expectation  values  are  quite  as  simple  as  For 

example,  the  expectation  value  of  the  fourth  parent  moment  about  the  mean 
is 

- n ^ 4 L 12  18  1 

^4  ^L1  (M-l)  + [21  + [31  J 


. 4 

12  18 

(M-l) 

+ [2]  + 

(M-l)1  J (M-l) 

2 M 

i 3M  1 

(M-l) 

(M-l)[3]  ^ 

■ * m t -y*  %•*%*  v «rl 
- - ' .f  * . 1^.  ..  . ' 


where  M is  the  number  of  observations  and 


The  formal  procedure  for  finding  unbiased  estimators  and  expectation  values 


is  left  to  Appendix  3.  Here,  we  are  only  interested  in  giving  the  motivation 


for  such  formal  techniques.  In  Appendix  C we  employ  these  techniques  to 


find  the  expectation  values  of  at  and  Pt  as  given  by  Eqs.  (C8a)  and  (C8b) 


The  coefficients  on  the  right  hand  side,  according  to  Eqs.  (C9)  and  (CIO) 


We  can  substitute  Eqs.  (46)  and  (47)  into  Eqs.  (45a)  and  (45b)  to  obtain 


expressions  for  at  and  j3t.  That  is 
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A 2 
at  = — 


M 


ra„ 


3 (M-l)  2 

+ 


M/~  1 

2 V(M- 

K 


i) 


(M-l) 


[2] 


(M-l) 


my  m2 


1 + Trr-rr  + — ~ 


36 


\ 


(M_1)  (M-l)f2] 


(M 


.i)[3j;  m4 


(48a) 


A 1 
= ' 24  m2 


_ M A 1 


8 'V(M-D  - (M.1)  [2]  - (M.L)[3]y  m2 


ifi  + -Z_ 

24  V (M- 


24 


L)  (M-l)[2]  (M-l) 


36  N 
[3] ; m4  • 


(48b) 


The  diffusion  time  t can  generally  be  determined  with  arbitrarily  good 

A A A A 

accuracy.  Then  knowing  t,  at,  and  0t  gives  us  a and  0 directly,  as 

« A 
a = at/t 


5 A 

P = 0t/t 


(49a) 

(49b) 


We  now  know  how  to  estimate  the  jump  rates  a and  8 from  an 
hypothetical  data  set.  This  knowledge  is  of  little  value,  however,  unless 
we  can  also  estimate  the  random  error  incurred  in  making  the  measurements. 
The  next  problem  we  therefore  consider  is  that  of  finding  a measure  of  the 
random  errors  inherent  in  an  experiment. 


B.  Estimating  the  Random  Errors 

The  magnitude  of  the  random  error  in  the  random  variable  at  is 

indicated  by  the  variance,  defined  as 

A A A . 2v 

var  at  = ( (at  - (at/ ) ) . 


(50) 
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According  to  the  definition  of  an  unbiased  estimator  as  expressed  in 


Eq.  (42), 


(at)  = at  . 


equation  (50)  can  be  reduced  to  the  equivalent  form 


A A ? . 

var  at  = ( (at  - at)-) 


A 


A 


= ((at)  - 2at  at  + (at)  ) . 


The  expectation  value  of  a sum  is  equal  to  the  sum  of  the  expectation  values 

for  each  term  of  the  sum  so  that 

A A 2 A 2 

var  at  = ((at)  ) - (2at  at)  + ((at)  ) . 

Because  at  is  a population  value, 


(2at  at)  = 2a t (at)  = 2 (at)" 


and 


((at)2)  = (at)-  . 


I 


The  variance  of  at  is  thus  given  by  the  well  know  relation 

A , A 2\  2 

var  at  = ((at)  ) - (at) 

A 

Similarly,  the  variance  of  lit  is 

A A 2 2 

var  St  = ( (3t)  > - Qt)  . 


(51a) 


(51b) 


However,  the  variances  according  to  Eqs.  (51)  are  population  values,  and 

therefore  cannot  be  determined  from  a finite  experiment.  The  unbiased 
A A 

estimators  of  var  at  and  var  ,3t  though  are  random  variables  that  are 

accessible  through  a finite  number  of  observations.  Using  the  operator  U 

to  indicate  that  an  unbiased  estimator  is  being  taken,  we  write  the  unbiased 
A 

estimator  of  var  at  as 


A / A A 2 

U(var  at)  = U(((at)  ) - (at)  ) . 


(52) 


■ 'tr  *■ 


. 
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The  operator  U acts  linearly,  so  that 

U (var  at)  = U[((at)2>]  - u[(at)2]  . 

The  untiased  estimator  of  the  expectation  value  of  a random  variable  is 
equal  to  the  random  variable  itself  and  therefore 

U[<(at)2>]  = (at)2  . 


An  equivalent  expression  for  U(var  at)  is  then 

A A 2 2 

U (var  at)  ^ (at)  - U[(at)  ] 


(53a) 


Similarly,  the  unbiased  estimator  of  the  variance  of  gt  is 

A A 2 A 2 

U (var  gt)  = (0 1 ) - U[(gt)  ] . 


(53b) 


The  error  in  measuring  t is  negligible  compared  with  the  error 

A A 

in  a and  3 . That  is , 


U fvar  a]  _ Ufvar  at] 


(a)' 


(at  y 


- A 

Ufvar  3 ] ~ U f var  g 1 1 

* 2 a 2 

(g)  (pt)z 


(54a) 


(54b) 


Substituting  Eqs.  (45a)  and  (45b)  into  Eqs.  (53a)  and  (53b)  respectively,  we 
arrive  at 

Ufvar  at]  = (Rm2  + Sam2  + T^)2  - U[(at)2] 

= Ram2  + 2RaSam2  + 2RaTam4m2  + Sam2 


+ 2S  T m.m2  + T2m2  - U[(at)2] 
ot  a 4 2 a 4 


(55) 


where  the  coefficients  R^,  S^,  and  T have  already  been  defined  by  Eqs.  (46a)- 

2 

(4oc).  In  Appendix  C,  we  find  that  U[(at)  ] can  be  expressed  in  terms  of 


f '-XT  ' 

inrt^jiAn  iii  -|%r 


— - . . i iriiMiailill  • I . — -»■ — . - — 1-^.-  w . 
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i 
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f i 

I t 
1 


5 k 

Ih 


i 

i 

i 


parameters  <77,  ^42’  an<3  ^44’  known  as  '^-scaciscics , as 
U [ (ore ) 1 * 9 aZ2  ’ 9 k42  + 36  k44‘ 


As  shown  in  Appendix  C, 

2 2/  1, 

TT  + 

M 


k22  = a2*\  [2]  + J3]  + [4]  ) ( 


4 + 6 \ 


VMC2]  ^3]  m[4]  J 


16 

114 

432 

720  \ 

V mC2] 

- mC3]  * 

’ m[5] 

' x[6]J 

2^  1 

14 

Xi  * 4b 

37 

. 324 

. 540  ^ 

' W2] 

' M[3j 

m[4] 

V5] 

. l.p~f  4 , 38  . 144  240  "\ 

a^VM[3] + ^4i  - jsr^isTj 


&(- 


13 _54_  _ 90  N 


V M[5]  ' M[6]^ 


, f 1 

*44  = ®S‘MC 
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[2 

w i”  • 


M 


[31 


. 438  _ 4320  _ 27360  _ 103680  _ )AV±Q\. 


X 


£*] 


M 


[5] 


M 


[6] 


M 


[7] 


M 


[a]  ; 


m„a,M 
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256  2736  _ 13000  ^ 59120  ^ 120960  ^ 


W31  mW  ^m£5]  ^ v[6]  ^ J7]  ' j8]  J 


, JLf  8 . 208 

' 5 1f  ( r 3 7 ‘r  [4 ] 


2304 


+ m7M 


5 3 Vm[3] 

Zf.  _L 


M 

14 


M 


£5] 
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M^31  ^ 
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,£7J 
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mm 
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mm 
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34 
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m[73 
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[3]  ) 


(56) 


(57a) 


(57b) 


- - ' 


•T*  • 


By  substituting  Eq.  (56)  into  Eq.  (55),  we  can  estimate  the  variance  of 
A 

at  in  terms  of  sample  moments  about  the  mean  and  the  k-statistics . That 


is. 


U [var  at]  = ^m2  + 2 Vam2  + 2 Vara4ra2  + Sam2  + 2SaTam4m2 

^ J.  2 4 , 2 , 1 , 

+ Tara4  9 k22  9 k42  " 36  k44 


(58a) 


The  coefficients  R^,  S^,  and  TQ  are  given  by  Eqs.  (46a),  (46b),  and  (46c). 

A 

Similarly,  the  variance  of  j3t  can  be  estimated  by 

A 2 2 3 2 4 2 

U[var  Pt]  = Rgra2  + 2RgS^m2  + 2RpTgm^m2  + S^m2  + 2SgTgm^m2 

+ T8m4  " 576  k22  + 288  k42  " 576  k44  (58b) 


where  Rg , Sjj , and  Tg  are  given  by  Eqs.  (47a),  (47b),  and  (47c).  The 
details  are  available  in  Appendix  C. 

Equations  (58),  (57),  and  (54)  allow  us  to  estimate  the 

A 

variance  of  a and  ? from  a hypothetical  data  set,  consisting  of  M obser- 
vations of  adatom  positions.  We  do  not  yet  know,  however,  what  experi- 
mental conditions  yield  the  smallest  variances  for  a and  §;  nor  do  we 
know  how  many  observations  are  needed  to  insure  statistical  errors  small 
enough  to  make  or  and  p meaningful.  These  problems  are  addressed  in  the 
next  section. 


C.  Experimental  Conditions  Minimizing  Random  Errors 

To  minimize  the  random  error  in  determining  unbiased  estimators 
for  the  rates  of  single  and  double  jumps,  we  consider  the  variance  predicted 
for  different  combinations  of  the  rates  a and  3«  We  can  calculate  these 
variances  from  Eqs.  (51a)  and  (51b): 


¥ 


A z A 2\  2 

var  at  = ((at)  ) - (at) 

A A 2.  2 

var  0t  = ( (0t)  ) - (Pt)  . 

A 2 > A 2 

The  evaluation  of  ( (at)  ) and  ( ($t)  ) requires  techniques  developed  in 

Appendix  B,  and  is  detailed  in  Appendix  C.  It  is  convenient  to  express 

A A 

the  variances  of  at  and  pt  in  terms  of  the  cumulants  k^,  v.^,  and  Hg 

which  are  in  turn  defined  in  terms  of  the  parent  moments  p. P-g>  and 

in  Eqs.  (60a)- (60b).  That  is, 
o 

16K6K2  + mt.33.)  ,2  + U « H *2 

36  1 M + (M-l)  (M-l)  4 + }[2]  4 2 

24M (M+l ) 4 _ _6  . 8 (M  + 7 ) 

(m-d[3]  2 M (“-1)  42 


16  ..  . 16  (M  + 1)  2 , , . „ 2 

+ ' (M-l)  *2  } ' (at) 


(59a) 


16h,h. 


H 3At  = -L-  [ -1  + 1:6  2 + ..(M+33)  2 

r '3C  576  1 m + /vi- 1 \ + /vi_  n **•/.  + 


(M-l) 


(M-l)  4 


(M-l) 


[2]  h4H2 


24M (M  + 1) 
(M-l) 


4 ^6 

*2  ‘ M 


2 (M+  7)  *4*2 

(M1!) 


, , (M+l)  2 3 2 

M (M-l)  2 * (Pt)  - 


(59b) 


The  properties  of  cumulants  required  for  this  work  are  discussed  in  Appendix 
B.  Cumulants  are  related  in  a simple  way  to  parent  moments.  Thus, 


r i 

i 

■ i 


*4  = ^4 


(60a) 


(60b) 
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K6  = "6  * 15^2  + 30i2 


(60c) 


*8  = ^8  ' 28V2  ‘ 354  4 + 42CV2  - 63aV 


(60d) 


A A 

Though  the  exact  expressions  for  var  at  and  var  pt  are  very  complex, 


the  corresponding  approximations,  to  first  order  in  1/M  are  fairly  simple. 
In  Section  4 of  Appendix  C,  Eqs.  (C20)  and  (C22),  we  show  that  these 
approximations  are  given  by 


var  = 3k  ^8  ' 12>W  " 1W,A  + 17  A + 136m'41j'; 

- 72P-3  - 8^6  + 16M-4  - 32p,2  ) 


(61a) 


var  Pt  576M  ^8 


12M.6^2-  2 n2  + 54^-  45^-  2u6 


3 2 

+ 16m*4P-2  - 18p-2  + m-4  - 2p.2  ) . 


Both  of  these  approximations  are  within  6%  of  the  corresponding  exact 


(61b) 


variances  for  any  value  of  3/a  ranging  from  0 to  1000,  as  long  as  M > 100 
and  n2  < 10. 

We  can  find  the  parent  moments  in  terms  of  a,  3,  and  t from  the 
generating  function,  using  Eqs.  (25)  and  (29).  These  tedious  calculations 


' f I 


are  preformed  in  Appendix  A and  yield  the  following: 


U2  - 2 (43 t + at) 


l*4  * ?y,2  + 3(p2)  " 12<2et  + at) 

p.  = 1283 1 + 2at  + 6 0 (43 1 + at)2  + 7203t(43t  + at) 
0 

+ 120 (43t  + at)3 


(62a) 


(62b) 


(62c) 


= 5123t  + 2at  +64,512  (3t)2  + 252(at)2  + 12,096  (at)  (3 1) 

O 

+ 20,160  (43t  + at)2+l680(43  +at)3  + 1680(4Pt+at)4. 


(62d) 


_ i* * j * - ^ a _ si  i 
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The  aim  of  all  this  is  to  find  the  random  error  inherent  in 
/s  * A A 

estimates  of  the  jump  rates  a and  3.  The  variances  of  at  and  0t  are  given 
by  Eqs.  (59),  (60),  and  (62).  Inasmuch  as  the  error  in  the  diffusion  time 

A A 

t is  insignificant,  the  error  in  measuring  a and  3 is  identical  to  that  of 

A 0 

at  and  Pt;  that  is. 


var  a 

(a)2 

A 

var  3 

~ 2 
<P) 


A 

var  at 
2 

(at)4 

A 

gar  gt 
A 2 

(Pt) 


(63a) 


(63b) 


We  recall  from  Eqs.  (34a)  and  (34b),  that  at  and  pt  can  be  expressed  in 
terms  of  p,  and  B/a  as 


at  = r 


2 (1+40 /a) 


(64a) 


1 ^2  3/0 
Pt  - T — 


2 (1+40/a)  * 


(64b) 


The  relative  error  j /a  and  ctq/3  in  a and  9 is  defined  as 

a p 


'“a  _ .var  a.  1/2 

a ~ -2  ' 

a 


(65a) 


I 

I 

1 

1 


a0  _ ,var  3,1/2 

s ( §2 


(65b) 


Equations  (51),  (59),  (60),  and  (62)-(65)  allow  us  to  calculate 
exactly  the  relative  error  in  a and  8 as  a function  of  for  given  values  of 

M and  3 /a.  These  calculations  were  performed  by  the  FORTRAN  program  KPLOT 
listed  in  Appendix  D.  Figure  3 shows  o^/a  and  J^/3  for  3/a  = .1  and  .2  for 
M = 100,  1000,  and  10,000.  From  Fig.  3 and  Tables  1D-3D  in  Appendix  D,  we 


M = 100 


! 


see  that  o^/a  and  O'jj/P  achieve  their  minimum  values  ir.  the  neighborhood  of 

p.?  = .237  for  the  values  of  M and  9/2  given  above.  It  is  for  this  value 

of  the  second  moment  that  measurements  of  the  jump  rates  a and  9 are  the  most 

advantageous.  Using  the  FORTRAN  program  MPLOT  to  calculate  a /a  and  o /3  as 

a 9 

functions  of  M for  = .237  and  3/a  = «1  and  .2,  we  find  that  only  after  a 
few  thousand  observations  have  been  taken  does  the  relative  error  drop  to 
407,  as  shown  in  Fig.  4.  For  M > 5000,  further  increases  in  the  number  of 
observations  have  little  effct.  Precise  determinations  of  the  rate  of  double 
jumps  therefore  appear  to  be  difficult;  however,  it  should  be  possible  in 
experiments  of  reasonable  length  to  ascertain  semiquantitatively  if  double 
jumps  participate  to  any  significant  extent  in  the  diffusion  of  single  atoms 
over  a surface. 


Number  of  Observations  M/IOOO 


Figure  Relative  error  in.  jump  rates  as  a function  of  M 
number  of  adatom  diffusion  intervals  observed. 


APPENDIX  A 


CALCULATION  OF  PARENT  MOMENTS 


I.  General  Approach 


The  aim  of  this  appendix  is  to  derive  relations  for  the  moments  of 


the  parent  population  of  adatom  displacements  in  terms  of  the  jump  rates  3 


finding  the  parent  moments  about  the  mean  p, 


Equations  (21) 


(25),  and  (29)  give  us  an  algorithm  for  deriving  the  k-th  parent  moment 


about  the  mean  u 


or  equivalently 


Though  the  calculations  for  p.  and  p.,  are  easy  enough  to  perform 


finding  p.  and  p.  requires  a bit  of  forethought  and  careful  bookkeeping.  In 


order  to  facilitate  the  bookkeeping,  we  introduce  the  definitions 


F=F(t,z)=  (0(z2-t-z2)+a(z  + z 1)-2(a  + 0)Jc 


F " 5J  F(t,z) 

F - (^)2  F(t,z) 

F(k)  ■ (^)k  F(t ,z)  . 


(A7a) 


(A7b  ) 


(A7c) 


The  procedures  for  utilizing  Eqs.  (A2)-(A7)  to  find  P-,  p.  , p.  , and  pQ  are 

Z 4 o o 

outlined  in  the  next  section. 


2.  Outline  for  Moment  Calculations 

a.  Find  T(t,z)  through  TgCt.z)  in  terms  of  G(t,z)  and  G^8)(t,z) 
using  the  recursion  relation  (t,z)  = (z  ■—)  F^(t ,z) . 

b.  Find  G(t,z)  through  Gv  (t,z)  in  terms  of  F(t,z)  through 

F^(t,z)  using  the  recursion  relation  G^k+L\t,z)  = (S/dz)G^  (t,z). 

• (8) 

c.  Find  F(t,z)  through  F (t,z)  in  terms  of  t,z,a,  and  0 using 
the  recursion  relation  F^k+1^(t,z)  = (d/oz)F^  (t,z) . 

d.  Evaluate  F(t,l)  through  F^8\t,l)  in  terms  of  a,  0,  and  t. 

e.  Evaluate  G(t,l)  through  G^(t,l)  in  terms  of  F(t,l) 

through  F^(t,l). 

f.  Evaluate  f^(t,l)  for  k= 2,4,6,  and  8 first  in  terms  of 

F (t , 1)  through  Fv  (t , 1)  and  then  in  terms  of  a,  P,  and  t. 

g.  Keeping  in  mind  that  p,  = T.  (t,l),  find  p , p , p and  p0  in 


2 9 r~  / 9 r J 

4 6 


terms  of  a,  3,  and  t. 
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3.  Determination  of  F ^ (t,z)  through  F ^ (t ,z) 

We  use  equation  (A4)  to  generate  the  following  recursion  relation 


rk+l(t,S)  - (s  JJ)  rk(t,«)  . 


(A8) 


This  recursion  formula  will  be  used  to  find  ~^(t,z)  through  Fg(t,z)  in 

(Q\ 

terms  of  G(t,z)  through  G*1  '( t,z ).  For  the  purposes  of  these  calculations 
we  introduce  the  symbols 


G00  * G(k>(t,z) 


rk  = Fk(t,z)  . 


Using  these,  we  now  write  the  F^'s  as 

ri  ■ It5  6 “ Z‘G 


r2  = (Z  l?)  ri  = (Z  ll)  (z®} 


(A9a) 


(A9b) 


(AlOa) 


= z (G  + zG  ) 


= zG  + z G 


(AlOb) 


; } 


‘3 


F = (z  f^)  (zG  + Z"G) 


= zG  + z2G  + 2z2G  + z3G(3) 


• 2"  3 (3) 

= zG  + 3z  G + z ’ 


(AlOc) 


I 


r,  = 


3 • 2“  3 (3) 

(z  ~)  (zG  + 3z  G + zV 


• 2"  2-  3 (SI 

= (zG  + z G)  + (6z  G + 3zJGv  ') 

+ (3z3G(3)  + zV4>) 

= zG  + 7z2G  + 6z3G^3'>  + z4G^4') 


(AlOd) 
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= (z  |j)[,6+7z2G  + 6«3GTn  ' + .V4)] 

= (zG  + z2G)  + 7(2z2G  + z3G(3)) 

+ 6(3zV3)  +zV4))  + (4zV4)  +z5G(5)) 

= zG  + 15z2G  + 25z3G(3)  + 10z4G(4)  + z3G(5) 


(AlOe) 


T6  = (z  |^)(zG  + L5z2G  + 25z3G(3)  + 10z4G(4)  + z5G(5)) 

= (zG  + z2G)  + 15  (2z  2(5  + z3G(3)) 

+ 25(3z3G(3)  + zV4>)  + 10(4z4G(4)  + z5G(5)) 


+ (5zV5)  + z6G(6)) 


= zG  + 31z^G  + 90z3g(3)  + 65z*G 
+ 15z5G(5)  + z6G(6) 


.4„(4) 


(AlOf) 


Ty  = (2  ^-)(zG  + 31z2G+90z3G(3) +65z4G(4) +15z5G(5) +z6G(6)) 


. 2--  2-  3 (3) 

(2G  + z G)  + 31(2z^G  + zJG(  J) 


+ 90(3z3G(3)  + z4G(4))  + 65(4z4G(4)  + z5G(5)) 
+ 15(5z5G(5)+z6G(6>)  + (6z6G<6)+27G(7)) 

= zG  + 63z2G  + 301z3G(3) 

+ 350z4G(4)  + 140z3G(5)  + 21z6G(6)  + z7G(7) 


(AlOg) 


« I 
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- 3*2"  3 

8 = (z  f^)(zG  + 63z  G + 301z  ; 


+ 350z4G(4)  + 140z5G(5)  + 21z6G(6)  + z?G(7)) 


• 2-  2-  3 C3-) 

= (zG  + z G)  + 63(2z  G + z G^  ') 


+ 301(3z3G(3)  f z4G(4))  + 350(4z4G(4)  + z^G^) 


+ 140 (5z  5G(5)  4-  z6G(6))  + 21  (6  z V6>  + z7G(7)) 


+ (7z7G(?)  H-  z8G(8>) 


zG  + 127z2G  + 966z3G(3) 


+ 1701z4G(4)  + 1050z5G(5)  + 266z6G(6) 


+ 28z7G(7)  + z8G(8) 


This  completes  the  first  stage  of  the  derivation. 


(8) 


4.  Relations  for  G(t,z)  through  Gv  ;(t,z) 


We  use  Eq.  (A5c)  to  generate  the  recursion  relation 


.(8) 


(AlOh) 


(All) 


This  recursion  formula  will  be  used  to  find  G(t,z)  through  G^  ' in  terms  of 
(8) 


F(t,z)  through  F (t,z).  For  the  purpose  of  these  calculations  we  simplify 
the  symbols  to  read 


F = F(t  ,z) 
,(k)  = -(k), 


(Al2a) 

F'"'  = F (t  ,z ) . (Al2b) 

Using  Eqs.  (A2),  (A5c),  and  (A6),  we  can  now  write  the  functions  as 

(A13) 


„(k)  .3  k 

G = ^ exP(F) 


The  recursion  relation  given  by  Eq.  (All)  implies  that  the  first  eight 
derivatives  of  G(t,z)  can  be  written  in  terms  of  the  first  eight  derivatives 
of  F(t,z)  as 


34 

• 3 

G = (exp  F)  = F exp  F (Al4a) 

G = |^-  (F  exp  F)  = F exp  F + F(F  exp  F) 

• • 9 

G = F exp  F + (F)  exp  F (Al4b) 

G(3)  = |^-  (F  exp  F + (F)2exp  F) 

(3}  ••  • 

= [F v 'exp  F + F(F  exp  F) ] 

+ [ (2FF)exp  F + (F)2(F  exp  F)] 

= Fv  'exp  F + FF  exp  F + 2FF  exp  F + (F)  exp  F 
C3  \ •• . .3 

= F ' 'exp  F + 3FF  exp  F + (F)  exp  F (Al4c) 

Gk  j ^ (F'  'exp  F + 3FF  exp  F + (F)  exp  F) 

= [F(4)exp  F + F(3;)(F  exp  F)] 

+ 3[F^  'F  exp  F + FF  exp  F + FF(F  exp  F)] 

+ [3(F)  F exp  F + (F)  (F  exp  F)] 

= F(4:>exp  F + F(3)F  exp  F + 3F(3;>F  exp  F 

+ 3(F)2exp  F + 3F(F)2exp  F + 3(F)2F  exp  F + (F)4exp  F 

= (FV ’ + 4F^  ;F  + 3(F)“+  6F(F)  + (F)  ) exp  F (Al4d) 

G(5)  = [F(5)  + 4(F(4)F  + F(3)F)  + 3(2FF(3)) 

(3  V . 2 • ••  2 *3" 

+ 6 ' (F)  +2F(F)  )+  4(F)  F]  exp  F 

(TlV  .2  ”2*  • 3 • S 

+ [FV  ;F  + 4FV  ;(F)  + 3(F)  F + 6F(F)  + (F)  ] exp  F 

= tF(5)  +4F(4)F  + 4F(3)‘f  + 6FF(3)  + 6F(3)  (F)2 

. ••  2 . 3-  (4).  (3)  . 2 

+ 12F (F)  + 4(F)  F + FWF  + 4FV  ' (F) 

+ 3(F)2F  + 6F (F)3  + (F)5}  exp  F 


(A14e 


= IF ' '+  5FV  'F  + 10F ' 'f  + 10F k ' (F) 
+ 15(F)2F  + 10F(F)3  + (F)5}  exp  F 


= {f(6^  + 5(F(J)F  + F(4)  F) 

+ L0(F(4)F  + (F^3))2)  + 10(F(4)(F)2+  2F(3)FF) 

(3)...  ..  3 (3)  . 3 ...  2- 

+ 15(2FV  'FF  + (F)  ) + 10 (F ' '(F)  + F3(F)  F)  + 5(F) 

r (5)*  (41  • 2 (3)"»  (31  • 3 

+ [f'‘  'F  + 5Fk  ' (F)  + 10Fk  'FF  + 10Fk  '(F) 

+ 15  (FF)2  + 10F(F)4  + (F)6}  exp  F 

= iF^6)  + 5F  ^ ^F  + 5F(4)F  + 10F^4)F 

+ 10(F(3))2  + 10F^4)(F)2  + 20F(3)FF  + 30F(3)FF 

+ 15(F)  + lOF*-  '(F)  + 30 (FF)  + 5F(F) 

+ f(5:)f  + 5f^4)(f)2  + iof^3)ff  + iof(3;)(f)3 


+ 15 (FF)2  + 10F(F)4  + (F)6}  exp  F 

= {f(G)+  6F(5)F  + 15F^4)F  + 10(F(3))2+  15F(4)(F)2 

(3  v ••  3 (31  *3  •••  2 

+ 60F 'J'FF  + 15(F)  + 20FV  '(F)  + 45(FF) 

+ 15F(F)4  + (F)6 } exp  F 


= [F(7>  + 6(F(6:>F  + F(5)F) 

+ 15(F^5)F  + f^F^)  +10(2F^3)F(4))  + 15(F(5^  (F)2+ 

+ 60(f(4)ff  + f(3)f(3)f  + F^3)FF) 

+ 15(3(F)2F(3))  + 20(F(4)(F)3  + F(3)3(F)2F) 

+ 45  (2)  (FF)  (F  ''J'F  + (F)  ) + 15  (F  ^ '(F)  + F4(F)  F) 

+ 6(F)5F}  exp  F 

+ 1f(6)F  + 6F(5)(F)^+  15F(4)FF  + 10(F<‘3))2F 
+ 15F(4)(F)3  + 60F(3)F(F)2  + 15(F)3F 
+ 20F(3)(F)4+  45(F)2(F)3+  15F(f)5  + (F)?}  exp  F 


4"  ■> 

FJ  exp  F 


(Ai4f 


2FFF(4)) 


■*(7) 


of 
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(6)4  * a-(5)' 

f r On  i 


,,.(5). 

Lo  r r 


L3F(4)F(3)  * 20F(4V3) 


+•  13F(;5)(F)2  + 3 OF  ^FF  +■  60F'4,FF 


6C(F(3))2F 


(3)  ••  7 
60F^  ; (F) 


+ 45(F)2F(3)  +■  20F(4)(F)3  + 60Fv3)F(F)2 
+ 15F(3)(F)4 

(4)..  • 

+ L5F^  ;FF  + 10 ( 


,>»^2  ,i,3  . ^ „ (6 ) i 


i3)„  .2  ••  3* 

90FV  ;F(F)  +■  90(F)  F 

6F(J)(F)2 


60(F)  (?)  + 6(F)  F + Fu  'F 

(3)  (4)  ■ 3 (3)..  • 7 ••  3- 

Fk  ')  F +•  15F ' ;(F)  + oOF'1  ;F(F)  + 13(F)  F 

(3)  • 4 ..7.3  ••  • 3 • 7 , 

20Fw;(F)  + 43(F)"(F)  + L5F(F)  + (F)  } ex?  F 

[F<‘/')  + 7F^F  + 21F<‘5')r  + 35F^F^  -t-  21F<'5'(F)2 

I05F('4)FF  t 70(F(3))2F  + 105F(3)(F)‘ 


+ 210F(3)F(F)2^  105(F)3F  + 35F(3)(F)4 
+•  21F(F)3  + (F)2  } exp  F. 


(41)  • 3 
r 35FV  ; (F) 

..?.3 

105 (F)  (F) 


(Al4g) 


It  is  now  convenient  to  define  the  function  H-,,  so  that 


• (ft) 

5.  Relations  for  F(t,z)  through  F y(t,z) 

• (S') 

In  order  to  find  F(t,z)  through  F (t,z)  in  terras  of  t , z,  ar, 
and  3,  we  need  only  take  the  first  eight  derivatives  of  F(t,z).  The 
function  F(t,z)  is  defined  by  Eq.  (A6)  as 


F(t,z)  = [S(z2  + z"2)  +ot  U + Z_1)  - 2(a  + g)]t  . 

As  a way  of  simplifying  Eq.  (A6)  we  define 

A 5 art 

(Al5a) 

B * St  . 

(Al5b) 

These  definitions  result  in  a new  expression  for  F(t,z): 

F (t , z)  = 

(A16) 

B[z2+z  “]  + A[z  + z *■]  -2(A+B)  . 

The  first  eight  derivatives  of  F(t,z)  with  respect  to  z 

now  are  given  by 

F = 3 [2z  - 2z"3]  + A[l-z”2] 

(A17a) 

F - B(2  +6z"4)  + 2Az'3 

(A17b) 

F(3)  = -24Bz‘5  - 6A z"4 

* -6  [4Bz  3 + Az  4] 

(Al7c) 

F(4)  = 24[5Bz'6  + Az"5] 

(Al7d) 

F(5>  = -120[6Bz'7  + Az'6] 

(Al7e) 

Fr6)  = 720 [7Bz'8  + Az'7] 

(A17£) 

F(7)  = -5040 [8Bz~9  + Az‘8] 

(Al7g) 

F(8)  = 40,320[9Bz*10  + Az'9] 

(A17h) 

* 
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• (8) 

6.  Evaluating  F(t,z)  through  F (t,z)  at  z = 1 

Sections  3 through  6 give  us  all  of  the  information  necessary  to 

find  T (t,z)  through  i Q(t,z)  in  terms  of  t,  z,  a,  and  3.  In  order  to  find 
1 o 

p..  , u , and  y.Q  in  terms  of  t,  a,  and  8;  we  only  need  to  find  T (t,z), 
r4(t,z),  r6  (t,z),  and  Tg(t,z)  at  z = l.  Because  all  the  r^'s  are  in  terms 
of  the  G^'s,  and  an  the  G^'s  are  in  terms  of  the  F^'s,  the  F^'s 
must  first  be  evaluated  at  z = 1. 

Substituting  z = 1 in  Eqs.  (Al7a)- (A17h)  and  (A16),  we  find  that 


I 


* 

! 

I 

1 

i 


F (t , 1)  = b*2  + A* 2 - 2(A  + B)  = 0 

(A18) 

F(t , 1)  = B(2-2)  + A(l-l)  = 0 

(A19a) 

F(t , 1)  = B[2  +6]  + 2A  = 8B  + 2A 

(Al9b) 

F(3)(t,l)  = -6  [4B  + A] 

(Al9c) 

F(4)(t,l)  = 24  [5B  + A] 

(Al9d) 

F(5)(t,l)  = -120  [6B  + A] 

(Al9e) 

F(6)(t,l)  = 720[7B  + A] 

(Al9f ) 

F(7)(t,l)  = -5040 [8B  + A] 

(Al9g) 

F(8)  (t,l)  = 40,320[9B  +A]  . 

(Al9h) 

• (8) 

7.  Evaluating  G(t,z)  through  G (t,z)  at  z = 1 

* (8) 

Rather  than  immediately  expressing  G(t,l)  through  G (t,l)  in 

terms  of  a,  3,  and  t,  we  choose  to  first  express  these  functions  in  terms 
* (8) 

of  F (t , 1)  through  Fv  y(t,l).  We  find  that  Eqs.  (A14a)  through  (A14h)  are 
considerably  simplified  by  making  the  substitutions  F(t,l)  = 0,  F(t,l)  = 0, 
and  z = 1.  That  is , 

G(t , 1)  = exp  F (t , 1)  = 1 (A20) 

G(t , 1)  = F (t , 1)  * 0 (A21a) 


HI  I 
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G(t,l)  = F (t , 1) 

(A21b) 

g(3) (t , 1)  = f(3)  (t,l) 

(A21c) 

G(4)  =F(4>+3(F)2 

(A2 Id) 

g(5)  = f(5>  + iof<3)f 

(A21e) 

c(6)  = F(6)+  15F(4)‘f  + 10(F(3))2  + 15(F)3 

(A21f ) 

G(?)  = F(?)+  21F(5)F  + 35F(4)F(3)+  105F(3)(F)2 

(A21g) 

F^  + 7F^6)  F + 21F(6)  F + 21F^5)F(3)+  35F(5)F(3)  + 

35  (F (4))2 

+ 105F(4)  (F)2  + 70(F(3))2F  + 105F(4)  (F)2  +210(F(3))2F  + 105  (F)4 
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G + 

31G  + 

90G(3)+65G(4) 

+ 

15G<5)  + 

g(6) 

0 + 

3 IF  + 

90F(3)  + 65 

(F 

(4) 

+ 3(F)2) 

15  (F 

,(5)  + 

lOF^F)  + 

(F 

(6). 

b 15F^F 

+ 10(F(3))2 

31F 

+ 90F 

(3)  + 65F(4) 

+ 

195(F)2  + 15F(5)  + 150F 

F(6) 

+ 15F^F  + 10  (F 

(3)}2 

+ 15(F)3 

(3); 


= 31  (8B  + 2 A)  + 90[-6(4B+A)]  + 65[24(5B+A)] 

+ 195  (8B  + 2A)2  + 15  [-120(6B  + A)  ] + 150  [-6  (4B  +A)  (8B  + 2A)] 

+ 720 (7B  + A)  + 15[24(5B+A)(8B  + 2A)]  + 10[-6(4B+A)  ]2+ 15(8B  + 2A) 
= 62  (4B  + A)  - 540  (4b  + A)  + 156  0(5B+A)  + 780(4B+A)2 

- 1800(6B  + A)  - 1800(4B+A)2  + 720(7B  + A)+  720(5B+A)(4B+A) 

+ 360  (4B  + A)2  + 120(4B  + A)3 

= 62(4B+A)  - 540(4B,+  A)  + 1560 (5B  + A)  - 1800 (6B  + A)  + 720 (7B  + A) 

+ 780 (4B  + A)2  - 1800(4B  + A)2  +360(4B  + A)2  + 720(5B  + A)  (4B  + A) 

+ 120(4B+A)3 

= -478 (4B  + A)  + (7800B+  1560A  - 10800B  - 1800A  + 5040B  + 720A) 

- 660(4B  + A)2  + 720(5B +A)  (4B +A)  + 120(4b+A)3 

* -1912B  - 478A  + 2040B  + 480A  - 660(4B  + A)2  + 720(4B  + A)2 


+ 720B(4B  +A)  + 120  (4B  + A) 

- 128B +2A  + 60(4B  + A)2  + 720B(4B+A)  + 120(4B+A)3 


(A22c) 


r (t,l)  - G+ 127G-I-966G(3)  + 1701G(4)  + 1050G(5) +266G(6) +28G(?) +G(8) 
- L27F  + 966F(3)  + 1701  (F(4)  + 3(F)3)  + 1050 (F(5)  + 10F(3)F) 

+ 266  (F(6)  + 15F(4)F  + 10(F(3))2  + 15(F)3) 

+ 28(F(7)  +21F(5)F  + 35F^4)F(3^  + 105F(3)(F)3) 

+ F(8^ +28F('6)  F +56F(5)F(3) +35(F^4))  2 
+ 210F(4)  (F)2  +280(F(3))2F  + 105(F)4 


. . J _ 


B 


- 
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= 127F  + 966F('3)  + 1701F('4)  + 5103(F)2 
+ 105CF^  + lO.SOOF^F 

+ 266F(6:>  + 3990F(4)F  + 2660(F(3))2  + 3990(F)3 
+ 28F(7)  +588F(5)F  + 980F(4)F(3)  +2940f(3)  (kT 
+ F(8)+28F<6)F  + 56F(5)F(3)+35(F(4))2 

+ 210F(4)  (F)2  +280(F(3))2F  + 105(F)4 
= 127  [2  (4B  + A)]  + 966  [-6  (4B  + A)] 

+ 1701  [24  (5B  + A)  ] + 5103[2(4B+A)]2 
+ 1050  [-120  (6  B + A)  ] + 10,500  [-6(43  + A)]  [2  (4B  + A)  ] 

+ 266  [720 (7B  + A) ] + 3990 [24 (5B  + A) ] [2 (4B  + A) ] 

+ 2660 [ -6 (4b  + A) ] 2 + 3990 [2 (4B  + A) ]3 
+ 28  [-5040(8B +A)]  + 588[-120(6B+A)]  [2(4B  + A)] 

+ 980  [24  (5B  + A)  ] [-6(4B  + A)J  + 2940  [-6  (4B  + A)  ] [2(4B  + A)]2 
+ [40,320  (9B  + A) ] + 28 [720 (7B  + A) ] [2 (4B  + A) ] 

+ 56 [-120(6B  + A)] [-6 (4b  + A)]  + 35 [24 (5B  + A) ]2 
+ 210[24(5B+A)]  [2(4B  + A)]2  + 280  [-6  (4B  + A)  ] 2 [2  (4B  + A)  ] 

+ 105 [2 (4B  + A) ]4 

= B [127 (2)4  - 966  (6)4  + 1701(24)5  - 1050(120)6  + 266  (720)7 

- 28(5040)8  + 40,320(9)  ] +A[127(2)  - 966  (6)  + 1701(24) 

- 1050(120)  + 266(720)  - 28(5040)  +40,320] 

+ B2 [5103(4)  (16)  - 10,500(12)  (16)  + 3990(43)  (20)  + 2660(36)  (16) 

- 588(240)  (24)  - 980(144)  (20)  + 28(1440)  (28)  + 56 (720)  (24) 

+ 35  (24)2  (25) ] 

+ A2 [5103 (4)  - 10,500(12)  + 3990(48) 

+ 2660(36)  - 588(240)  - 980  (144) 

+ 28(1440)  + 56(720)  + 35  (24)2] 
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; 


I 

* ■ 


t L 


+ AB [5 103 (4) (8)  - 10,500(12) (8)  + 3990(48) (9) 

+ 2660(36) (8)  - 588(240) (10)  - 980(144)  (9) 

+ 28(1440)  (11)  + 56  (720)  (10)  + 35(24)2(10)  ] 

+ (4b  + A) 3 [3990 (8)  - 2940(24)  + 280(72)] 

+ 210(96)  (5B  + A)  (4B  + A)2 
+ 105(16)  (4B  + A)4 
= 512B  + 2A 
+ 64,512  B2  + 252A2 
+ 12.096AB  - 18,480(4B +A)3 
+ 20 , 160 (5B  + A)  (4B  + A)2  + 1680  (4B  + A)4  . 

Note  that 

-18,480(4B  + A)3  + 20,160[(5B+A)  (4B +A)2] 

= -18,480(4B  + A)3  + 20,160[B(4B+A)2  + (4B  + A)3] 

= 1680(4B+A)3  + 20,160B(4B+A)2. 

This  makes  it  possible  to  write  r (t,l)  as 

O 

rg(t,l)  = 512B+2A  + 64.512B2  + 252A2  + 12.096AB 

+ 20, 16 OB (4B  + A)2  + 1680(4b  + A)3  + 1680 (4B  + A)4.  (A22d) 

This  completes  the  derivation  of  r2(t,l),  r^(t,l),  ^(t.l),  and  rg,(t,l). 

9.  The  Parent  Moments  in  Terms  of  a,  8,  and  t 

We  know  from  Eq.  (A3)  that  = Tk(t,l);  Eqs.  (A3),  (A15),  and 

(A22)  therefore  allow  us  to  find  the  even  moments  of  the  parent  population 

in  terms  of  the  jump  rates  or  and  |3  and  the  time  interval  t,  as 


1 


■> 


I 


'I  I 


M>2  = 2(43  + at) 

^4  = 7^2  + 3 (^2)2  - 12  (20 1 + at) 


(A23a) 

il 

(A23b)  j 

il 


I 


— 


s.  = 1283t + 2at +60(43t +at)2  + 7203t(43t +at) 
o 

+ 120  (48 1 +ort)3  (A23c) 


= 5123 1 + 2at  +64,512  (8 1)2  + 252  (at )2  + 12 , 096  (at)  (3 1) 

O 

+ 20,160(43t +at)2  + 1680(43t +at)3  + 1680(4pt +at)4  . (A23d) 

I u 


This  completes  the  derivation  of  the  moments  expressed  in  terms  of  the 
jump  rates  for  single  and  double  jumps. 
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APPENDIX  B 

UNBIASED  ESTIMATORS  OF  MOMENTS  ABOUT  THE  MEAN 


This  appendix  is  meant  to  serve  as  a tutorial  for  the  statistical 
techniques  employed  in  this  thesis.  The  major  source  of  this  material  is 
Kendall  and  Stuart's  book  on  advanced  statistics  [9],  The  reader  is 
assumed  to  have  a knowledge  of  statistics  at  the  level  of  Bevington's 
book  on  data  analysis  [10].  Topics  covered  here  include  power  sums,  augmented 
symmetric  functions,  cumulants,  and  k-statistics. 


The  Second  Moment 


As  an  introduction  to  the  statistical  techniques  used  in  finding 
the  variances  of  the  jump  rates  or  and  g,  we  will  consider  the  problem  of 
estimating  the  second  parent  moment  about  the  mean  given  a sample  con- 


sisting of  M observations.  It  is  well  known  that  the  unbiased  estimator 


of  u is 


M 


^2  = M-l  m2  * 


(Bl) 


Here  m2  is  the  second  moment  for  a sample  consisting  of  M observations. 


If  Eq.  (Bl)  correctly  estimates  then  Eq.  (42)  implies  that  the  expectation 


value  of  p,  is  ^2*  That  is, 


<|i2>  = n2  • 


(B2) 


Because  the  number  of  observations  M is  not  a random  variable,  we  know  that 


(B3) 


- - r *•  * « 1 %*■ 

* - «■  - * „ j *.  4 ‘ \ 

■ r»  1 V 'at r . 


J 
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Equation  (31)  can  be  verified  by  finding  the  expectation  value  of  m^. 

From  Eqs.  (38)  and  (39)  we  find  that  can  be  expressed  in  terms  of  the 
M observations  x^  comprising  the  sample  as 

1 M 1 M 2 

ra2  = q 1 (xi  - S z xi>  • (B4> 

4 11  i-i  1 " j-i  J 


Here  and  throughout  this  Appendix,  we  will  omit  the  bounds  on  the  summation 
symbol  and  assume  that  every  index  is  summed  from  1 to  M.  Equation  (34) 
can  then  be  written  equivalently  as 


m2  = M Z^xi  ' E Xj//M)2  • 

We  can  reduce  to  a simpler  form  by  expanding  Eq.  (B5) : 

m2  = M Z^Xi  * + (ZXj/M)2] 

= m Zxi  ’ “T  dxi)(2x.)  + (Ex  /M)2 
M~  J J 

= — Ex2  - — (Ex  )2 

M i ^ (Z  i' 


(B5) 


(B6 ) 


Because  the  expectation  value  of  a sum  is  the  sum  of  the  expectation  values, 


1 

I 

I 

I 

I 


OO  = < £ Ex2)  - < 4r  (Ex . )2> 


M i 


M 


= m 1 (xi}  ‘ *7  <ffxi>2>  • 


(B7) 


We  already  know  that 


(x2> 


M 


2 ’ 


the  first  term  of  Eq.  (B7)  is  therefore  obvious.  In  order  to  write  the 
second  expression  on  the  right  in  terms  of  parent  moments,  we  must  first 


a .1 


*■  A:  m 4 a.-*  - * f s »• 

M VBt  ktiRM a 


observe  that 


2 ? 

(Ex .)  = Ex.  + E x.x.  . 

1 1 L J 


(B8) 


Because  i^j  for  each  uncorrelated  product  x^ , assigning  M different  values 
to  i would  leave  only  (M-l)  values  to  be  assigned  to  j.  The  second  term  of 
Eq.  (B8)  therefore  consists  of  a sum  over  M(H-l)  uncorrelated  products 
x^.  By  substituting  Eq.  (B8)  into  Eq.  (B7)  we  find  that 


<m0> 


M Z ' 4 + 2 x x ) 

" L M"  1 i^j  L J 


^2  " \ (M<xJ>  + M(M-l)<x  x >] 
M J 


In- 

M ^2 


mi 

M 


/ 

\ 


x.x . > . 
i J 


(B9) 


The  expectation  value  of  the  product  of  statistically  uncorrelated  factors 
is  the  product  of  the  expectation  values  of  those  factors.  Therefore, 


and 


= ( x L>  ( x ^ ) 

(m  ) = m-.U.  u,  . mi  w_  \ 

{ 2?  M ^2  M UD'Xj' 


= mi  u.  . 

M ^2  M *1*1 


mi 

M 


o*;  - ) . 


(BIO) 


(Bll) 


From  Eq.  (20b)  we  know  that 

M>2  * < (x£  - p.p2> 

= <x2  - 2x^1  + u.^2) 

= (x2)  - (2x  n')  + <M.'2)  . 


I 

I 


: 


T 

i. 


I 

I 
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The  parent  mean  is  not  a statistic,  so  that 


/ \ 2 n I / \ ,2 

42  = (xt  - 2y.|<xi)  + 

= - 2p.'p,^  + 


■ *i  ■ “l2 


(B12) 


We  substitute  Eq.  (B12)  into  Eq.  (Bll)  and  find  that 

(m  > = SSkJL 1 p. 
v r m 2 ’ 

which  agrees  with  Eq.  (B3).  We  have  proven  Eq.  (BL),  though  not  in  the 
most  elegant  manner.  In  the  next  section,  we  introduce  techniques  that 
simplify  such  calculations. 

2.  Augmented  Symmetric  Functions  and  Power  Sums 


We  noted  that  by  writing  m„  in  terms  of  Ex.  and  E x.x., 

2 1 .,.11 


we 


could  easily  find  (nu)  in  terms  of  parent  moments.  From  Eqs.  (B6)  and 
(B8),  we  know  that 


Zx2  - £x2  + Z x.x.) 

M ’ ' i . , . l i 

i^j 


"2  " M -i  m2 


= -^4^-  Zx2  - ~~  Z x.x.  . 

M2  1 M2  i*j  1 J 


(B13) 


The  sums  Ex.  and  E x.x.  belong  to  a class  of  statistics  known  as  augmented 
X i^j  1 J 

symmetric  functions.  In  the  notation  associated  with  this  class  of 


statistics,  these  two  sums  are  represented  by 


[2]  = Z x‘ 


(B14a) 


[1  ] = Z XX  . 


h ■ i 


. -1  1 . Mi.  f ),  % | . - j 


(Bl4b) 
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In  general,  augmented  symmetric  functions  are  defined  as 

P1  p2  Ps 

fplp2  '*•  Ps]  S E Xi  xi  •**  Xv 


(B15) 


The  indices  i through  v are  summed  from  1 to  M,  though  in  no  term  of  the  sum 
are  any  two  or  more  of  these  indices  allowed  to  have  the  same  value.  As 
an  example  of  this  notation,  we  will  compute  the  augmented  symmetric  function 
[221]  from  the  four  element  data  set:  x^  = -1,  x ^ = 3,  x^  = -1,  and  x^  = 0. 
We  can  express  [221]  as 

[221]  = I x^xk 

where  i^j , j^k,  and  i^k  for  every  term  of  the  sum.  The  number  of  terms 
belonging  to  [221]  is  equal  to  the  number  of  permutations  of  4 objects  taken 
3 at  a time.  This  number  is 


= 4(4-1)...  (4-3  + 1)  = 4.3.2  = 24. 

These  24  permutations  and  the  corresponding  term  of  each  is  given  in  Table 
Bl  below.  There  we  find  that  [221]  has  a value  of  -30  for  this  data  set. 

There  is  a more  elegant  notation  for  expressing  augmented  symnetric 
functions , than  that  shown  in  Eq.  (B15).  For  example,  we  could  write  [221] 
as  [2?'l]  or  [333  2222  11]  as  [3^2^1^].  In  general,  augmented  symmetric 
functions  can  be  defined  by  this  notation  as 
X1  X2  Xs  _ 

[PL  P2  •••  Ps  I = 

P1  P1  p2  P2  Ps  Ps 

Z(x  •••  x,  ) (x  ...  x,  )...(x  ...  x s)  . (B16) 

a u e n m r 

where  indicates  the  number  of  powers  p^  in  each  term  of  the  sum.  No 
two  or  more  of  the  indices  a through  r,  as  they  are  summed  from  1 to  M, 


• 4 ”*  **  V-  v V- 

1 b « 


r 
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TABLE  Bl  - The  augmented  symmetric  function  [211]  computed 
from  the  data  set  [-1,  3,  -1,  0 ]. 


i 

j 

k 

x . 

l 

x . 

J 

*k 

2 2 
x.x . 
i J 

1 

2 

3 

-1 

3 

-1 

-9 

1 

2 

4 

-1 

3 

0 

0 

1 

3 

2 

-1 

-1 

3 

3 

1 

3 

4 

-1 

-1 

0 

0 

1 

4 

0 

-1 

0 

3 

0 

1 

4 

3 

-1 

0 

-1 

0 

2 

1 

3 

3 

-1 

-1 

-9 

2 

1 

4 

3 

-1 

0 

0 

2 

3 

1 

3 

-1 

-1 

-9 

2 

3 

4 

3 

-1 

0 

0 

2 

4 

1 

3 

0 

-1 

0 

2 

4 

3 

3 

0 

-1 

0 

3 

1 

2 

-1 

-1 

3 

3 

3 

1 

4 

-1 

-1 

0 

0 

3 

2 

1 

-1 

3 

-1 

-9 

3 

2 

4 

-1 

3 

0 

0 

3 

4 

1 

-1 

0 

-1 

0 

3 

4 

2 

-1 

0 

3 

0 

4 

1 

2 

0 

-1 

3 

0 

4 

1 

3 

0 

-1 

-1 

0 

4 

2 

1 

0 

3 

-1 

0 

4 

2 

3 

0 

3 

-1 

0 

4 

3 

1 

0 

-1 

-1 

0 

4 

3 

2 

0 

-1 

3 

0 

-30 


1 

I 


i ~ ■ i ' i i ’ i mli 


v X' 


j&jlj 
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may  be  equal  for  any  one  term  of  the  sum.  Specific  examples  of  this  notation 


[2 12  ] = I x2x,x 

a b c 


where 


r,,3. . p 6 2 2 2 
[62  1]  = Z xaxbxcxdxe 

.,2,3.  ^ 77222 

I7  2 ] - S xaxbxcxdxe 

^2 

Each  term  of  [p^  p2 


pE  Vv  •••  xs 


p ] is  a product  of  p factors 
s 


(B17) 


X1  X2  Xs 

Each  term  of  [p,  p,  ...  p ] then  corresponds  to  an  element  of  the  set 

of  permutations  generated  by  M objects  taken  p at  a time.  The  number  of  such 

^1  ^2  X s 

permutations  and  therefore  the  number  of  terms  in  [p  p ...  p ] is 

i 2.  S 


M^P  s M(M-l)  ...  (M-p  + 1)  . 


It  follows  that  any  augmented  syimietric  function  characterized  by  M 


(B13) 


observations  and  a given  value  of  p can  be  expressed  as  a sum  of  MLHJ  terms. 

The  expectation  value  of  such  an  augmented  symmetric  function  is 

X.  X,  X 

/ t 1 2 slN 

([Pi  p2  ...  ps  ]) 

P1  P1  p2  p2  Ps  pc  , 

= (.E(x^  • • »3Cj  ) (Xq  • • ) . • • (x^  • • .x^  )/ 

r„i  Pi  Pi  p?  p?  P„  p 

= MlPJ<(x1  ...Xj  )(x  ...xr  )...(xu  ...xv  s)>  . 

Each  of  the  p subfactors  is  statistically  independent  so  that 
X,  X X 

<[pl  P2  ...Ps  S]>  - 

r„i  pi  pi  p?  P?  Pc  P 

M[0j«x.  1>...(xj  1>)«xq  2>...<xr  2>)...«xu  s>...(xv  s>) 


= ^ (P- ' • • • M» T ) (M* 1 ...u*  ) • . . (jx 1 ...p.  ) 

P1  P1  P2  P2  Ps  Ps 


f * v-  » • t—  * « 

. * *.  * . . 


It  is  Theorem  (BI9)  for  the  expectation  value  that  makes  augmented  symmetric 


functions  so  useful.  Returning  to  the  example  discussed  at  the  beginning  of 


the  section,  we  can  write  m in  terms  of  augmented  symmetric  functions  by 


substituting  Eqs.  (Bl4a)  and  (Bl4b)  in  Eq.  (B13),  ieading  to 


There  are  easier  ways  of  finding  m„  in  terms  of  augmented  symmetric  functions 


than  that  employed  above.  Moments  can  be  expressed  in  terms  of  power  sum 


statistics  which  are  in  turn  tabulated  in  terms  of 


functions.  Examples  of  power  sums  are 


(B21b) 


The  r-th  power  sum  can  be  defined  as 


(B22b) 


as  can  be  seen  from  Eq.  (38).  From  Eq.  (B6),  we  know  that 
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We  can  evaluate  (1)“  in  terms  of  augmented  symmetric  functions  by  sub- 
stituting Eqs.  (Bl4a),  (B14b),  and  (B21a)  in  Eq.  (B8) : 

(l)2  = [2]  + [l2]  . 


(B24) 


By  comparing  the  definitions  for  augmented  symmetric  functions  and  power 
sums,  Eqs.  (Bl5)  and  (B22) , we  can  see  that 

[r]  = (r)  ; (B25) 

[2]  = (2)  . (B26) 

We  substitute  Eqs.  (B24)  and  (B26)  into  Eq.  (B23)  and  find  that 

m _ 1 [2]  _1_  ([2]  + [l2]) 

2 "M  'm2 


= (M-l)  [2]  _ _1_  [1“] 


M2 


(B27) 


M 


Equation  (B19)  shows  us  how  to  find  the  expectation  values  of  [2]  and 

[l2], 

<[2])  = M ^ 

([12]>  = M(M-l)  M.[2  . 


The  expectation  value  of  m^  is  then 


(B28a) 

(B28b) 


<», > - t21  . i it). 


M 


M2 


(M-l)  < f 2 1 > _L  < [1  ]> 


M 


M2 


sttm  u.  . mzii  u.z 

2 M2  1 


M2 


= &-1)  , _ 2 

M ^2  ) * 
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We  arrive  at  the  same  result  for  <m > as  given  by  Eq.  (Bll). 

Equations  (B24)  and  (B26)  were  derived  only  for  heuristic  pur- 
poses. Ordinarily,  if  we  wanted  to  find  a power  sum  or  product  of  power 
sums  in  terms  of  augmented  symmetric  functions,  or  vice  versa,  we  would 
consult  a table  such  as  Table  B2  below. 


TABLE  B2  - Augmented  symmetric  functions  and 
power  sums  of  weight  2. 


[2] 

9 

[11 

(2) 

I 

-1 

a)2 

1 

I 

This  table  implies  that 

(2)  = [2] 

(1) 2  = [2]  + [l2] 

[2]  = (2) 

[l2]  = -(2)  + (l)2. 

To  express  a given  augmented  symmetric  function  in  terms  of 
power  sums,  read  the  corresponding  column  from  top  to  bottom,  down  to 
and  including  the  diagonal  element.  To  express  a given  product  of  power 
sums  in  terms  of  augmented  symmetric  functions,  read  the  corresponding  row 
from  left  to  right,  up  to  and  including  the  diagonal  element.  Note  that 
in  Table  B2  and  elsewhere,  each  diagonal  element  "1"  is  written  as  I: 
"capital"  one. 


(1 

u 
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Table  B2  is  said  to  be  a table  of  weight  2.  The  weight  of  a table 


is  the  same  as  the  weight  of  the  augmented  symmetric  functions  in  that  table. 

X,  X„  Xc 

12  s 

For  the  augmented  symmetric  function  [p^  . . . pg  J,  the  weight  p is 


defined  as 


P = + P2^2  + *•*  * 


(B29) 


The  weight  of  both  [2]  and  [1“]  is  2;  clearly  then  the  weight  of  Table  B2 
must  be  2. 


3.  Direct  Calculation  of 

So  far,  we  have  only  confirmed  a stated  relation,  i.e.  that 
^2  = Mm^/CM-l).  We  still  require  a routine  procedure  for  arriving  at  the 


unbiased  estimator  of  an  arbitrary  moment.  For  this,  must  be  reduced 


to  a sum  of  population  values  that  have  easily  derived  unbiased  estimators. 
The  most  obvious  class  of  population  values  possessing  this  property  are 
the  parent  moments  about  zero.  From  Eq.  (B19)  we  can  deduce  that  the 
unbiased  estimator  of  a product  of  parent  moments  about  zero  is 


X X 

U(p.'  v •••  M-  * ) = 


X X X r , 

[pl  P2  •••  Ps  S ] /M  P 


(B30) 


where 


p = + x2  + ...xs- 


From  Eq.  (B12),  we  know  that  can  be  expressed  in  terms  of  parent  moments 
about  zero  as 


m<2=  ^ - <^[)2 


Because  the  unbiased  estimator  of  a sum  is  equal  to  the  sum  of  the  unbiased 
estimators  for  each  term,  we  know  that 


— 


• * * fj  ry**  *-■  * v 


Equation  (B30)  allows  us  to  conclude  that 


-in  . nil 

M113  M[2] 

111  . IL2!.. 

M M(M-l)  ' 


Table  B2  implies  that 


[2]  = (2) 


so  that 


[l2]  = -(2)  + (l)2, 


111  . r -(2)  + (D2  ] 

M 1 M(M-l)  J 


= m . . nu2 

M M(M-l)  M(M-l) 


We  recall  from  Eq.  (B22)  that 

(r)  = Mm’  . 
r 

Equation  (B32)  can  then  be  written  in  the  equivalent  form 
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2 

Equation  (323)  implies  that  m^  = m^  - (mp  ; therefore 


^2  = M-L  m2 


There  are  other  ways  of  achieving  this  result.  First,  however,  we  must 
introduce  additional  functions  - curaulants  and  k-statistics  - which  reduce 
the  derivation  of  unbiased  estimators  to  a routine. 


4.  Cumulants 

We  have  seen  that  the  parent  moments  about  zero  are  both  simply 
defined  and  easy  to  use  for  finding  expectation  values  and  unbiased 
estimators.  There  is  a major  drawback  however,  in  their  application;  the 
value  of  any  moment  about  zero  is  dependent  on  the  choice  of  the  origin. 
Because  many  calculations  involve  quantities  that  are  independent  of  the 
choice  of  origin,  there  is  a need  for  an  origin- independent  population  value. 

We  can  gain  some  insight  into  this  problem  by  examining  the 
r-th  parent  moment  about  zero  p/  for  a continuous  probability  density 
function  f (x) , 


|j/  = j xr  f (x)dx 


(B33) 


Also  corresponding  to  f(x)  is  the  characteristic  function  0(t)  which  is 
defined  as 


0(t)  s j eltX  f(x)dx  . 
.00 


(B34) 


The  argument  of  the  exponential  is  complex,  so  that  the  integral  can  be 
defined  for  a wide  variety  of  functions.  We  can  expand  the  exponential 
in  Eq.  (B34)  and  obtain  a different  though  equivalent  form  for  0(t), 


f 


0(t)  - J*  (1+  + f<x)dx> 


CO  00 

= (J  f(x)dx)  + (J  x f(x)dx) 

-CO  -CO 


+ (J*  X2f (x)dx)  + ... 


(B35) 


by  substituting  Eq.  (333)  into  Eq.  (B35),  we  find  that 


>(t)  = + M-'  -^y1-  + ...  • 


(B36) 


I.  I 


We  can  use  the  characteristic  function  to  generate  u/  by  employing  the 
simple  algorithm 


Ur  = (_i  dt^  0(t) 


(B37) 


We  have  previously  remarked  that  one  of  the  disadvantages  of 
working  with  moments  about  zero  is  that  the  value  of  such  moments  depends 
on  the  choice  of  origin.  If,  for  instance,  we  were  to  augment  the  value 
of  the  origin  by  c,  the  corresponding  moment  p/ (c)  would  be 


^r(c)  = (_i  dF^  J eit(x‘c)f(x)dx| 


. . d . r -itc  r>  itx  . . . , 
(-i  J?)  e ^ e f (x)dx 


- (-i  ^)r  e'itC  0(t) 


(B38) 


In  general,  p.^(c)  promises  to  be  very  different  from  pi/.  Suppose  however, 
that  we  define  a new  population  value  by  the  cumulant  generating 
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function  log  0(t)  so  that 


Kr  = (-i  ^)r  log  0(t) 


(B39) 


t=0 


If  we  were  to  change  the  origin  by  c again,  but  this  time  generate  H^(c), 
we  would  find  that 

CO 

u . . . d xr  , it  (x-c) , . 

h^(c)  = (-1  — ) log(j  e dx) 


. d ,r,  . -itc  p itx 


= (-1  ■££)  log(e 


e dx) 


t=0 


t=0 


= (-i  ^-)r(*itc  + log  0 (t ) ) 


= [(-i^  )r-itc  +V\tmQ 


1 t=0 


or 


= -c  6.  + H 

lr  r 

= k for  r > 1 , 
r 


(B40) 


where  is  Kronecker's  delta.  These  population  values,  called  cumulants, 
are  invariant  with  respect  to  the  choice  of  origin  except  for  the  first, 
which  is  decreased  by  c. 

Cumulants  can  also  be  related  to  the  cumulant  generating  function 
by 

2 . r 

log  0(t)  = -^T1  *2  + ...  -^r1-  Hr  +...  , (B41) 


this  is  simply  a solution  to  Eq.  (B39).  We  can  substitute  Eq.  (B36)  into 
Eq.  (B41)  and  see  that  cumulants  are  related  to  moments  about  zero  through 


or  equivalently 


An  explicit  expression  for  the  r-th  cumulant  * in  terms  of  moments  about 


zero  can  be  found  by  expanding  the  logarithm  in  Eq.  (B43)  and  picking  out 


where  the  sum  is  over  the  partition  of  the  integer  r such  that 


sion  for  the  r-th  parent  moment  about  zero  in  terms  of  cumulants  can  be 


found  by  expanding  Eq.  (B44) 


t 

1 /“'  *1 

/m-'  \ 

, ,V(-l)p"l(P-l>!  1 pl 

\ I 

^2 

f 

i 

w 

r i-  V \Vj 

i 1 

F 

wi 
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and  then  picking  out  the  terms  in  the  exponential  expansion  which,  when 
multiplied  together,  give  a power  tr.  We  therefore  have 


/ H ^ i X X.  - 

• ^ r!  ( P1  1 _?2  2 

k2l  •••  l Pi!  P7'  ' 


(B47) 


/ \ / 


where  the  sum  is  over  the  partition  of  the  integer  r such  that  p ^ + 


* * * ~ t • 


By  choosing  the  origin  of  the  distribution  so  that  M-j  = 0,  the 
cumulants  can  be  written  in  terms  of  moments  about  the  mean.  For  example 

*4  = ^4  * 4^4'  - 3^2  + 12P^'2  - 6M-J4 
2 

= Up  - 3U 

4 "*2 

Equation  (B47)  can  be  used  to  show  that  =M-|.  If  we  assume  that  nj  = 0, 
Che  expressions  for  parent  moments  about  the  mean,  in  terms  of  cumulants, 
are  often  considerably  simpler  than  the  corresponding  expressions  for  parent 
moments  about  zero.  For  example,  for  the  fourth  moment  about  zero,  we  have 

^4  = *4  +4k3“i  + 3K2+6K2Ki  + *14; 
for  the  moment  about  mean,  however, 


L.  l 


“4  ■ \ + • 


Therefore  from  Eq.  (B47),  we  can  write  the  second  through  eighth  parent 
moments  about  the  mean  in  terms  of  cumulants  as 


' 
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n2  = s 


^3  = *3 


H,  » *,  + 3*~ 
4 4 2 


M*.  = *e  + 10^.»h 
do  3 2 

2 3 

Pg  = *g  15^4^2  + 10x3  + 15*2 


ll  7 = *7  + 21k5h2  + 35*4*3  + 105*3*2 

M*0  = Kg  + 2®'"’5X?  + + 35*4  + ^ ^’4*2 

2 4 

+ 280.  + lOS*^  * 


Similarly,  Eq.  (B45)  can  be  used  to  show  that 


(B48a) 
(B48b) 
(B-  ic ) 
(B48d) 


/»>/  O - \ 


(B48f) 


(48g) 


k.  = M- 
2 2 


* = u 

3 ^3 


* = d,  - 3p.; 

4 4 ^ 


*3  - (J'j  “ lOpgi^T 


K6  = ^6  15W  1(^3  + 3(Xi2 


*7  = M-7  - 2hi5.u2  - 35d.4p-3  + 21ft43p.2 


*8  ^8 


28|J.  iju  - 56(J.  p.,  “ 35M-?  + 

6 2 5 3 *+  4 2 


+ 560u3^2  - 630M2 


(349a) 

(B49b) 

(B49c) 

(B49d) 

(B49e) 

(B49f) 

(B49g) 


5.  k-statistics 

So  far,  we  have  only  shown  how  moments  are  related  to  cumulants. 
We  have  not  yet  shown  how  cumulants  can  be  used  to  find  the  unbiased 
estimators  of  parent  moments.  Because  every  parent  moment  can  be  expressed 


1 


b 


%**•  t*~ 

^ * a ^ 
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as  a sum  of  cumulant  products,  a way  of  finding  the  unbiased  estimator  of 
a cumulant  product  could  facilitate  the  search  for  an  unbiased  estimator 
of  a parent  moment.  Equation  (B45)  gives  us  an  expression  for  a single 
cumulant.  The  equivalent  expression  for  a product  of  cumulants  x x x ... 


x x ...  =r;s:  ... 
r s 


r A i /**'  \ x,/ 1*  \ x, 
- :/  pi\  Lf  2 

— 1 V X2 ’ Pr  ' \ P2 ’ 


(B50) 


the  symbol  tt  indicates  that  we  are  taking  the  product  of  the  sums  corre- 
sponding to  the  cumulants  that  comprise  the  cumulant  product  x x ...  . 

In  Section  10,  we  find  that  this  product  can  also  be  expressed  as 


-D1  (P'^hKp-i):}  — / _ 
tti(x1:x2:  ...)j  1 ; \ P 


(B51) 


where  for  the  index  of  the  first  cumulant  product,  namely  r,  we  have  a 
partition  of  the  integer  r such  that  p^X^  + p?\2  + . . . = r and  X^  + X?  -I-  . . . = p , 
There  are  similar  partitions  of  the  indices  s,t  ...,  and  inside  the  large 
summation  the  symbols  tt  and  I are  used  to  denote  multiplications  and  sums 
to  cover  all  of  r,  s,  t ...,  while  the  large  summation  sign  denotes  that 


we  must  sum  over  all  partitions  [11], 


Now  that  we  know  how  to  represent  a cumulant  product  as  a sum 

of  moment  products,  we  can  use  Eq.  (B30)  to  find  the  unbiased  estimator 

of  x x . . . ; 
r s 


"Up,  I)  (p?:>  ...} 


(-1)^ (*° ' - ) rrf(o-l)13  t"(P^  * ?2 

ncxii  \2[  •••>}  ‘ ; ' , ^ pq  ; 


^1  /V*  2 

na>l:)  (p2:)  ...} 


We  call  the  unbiased  estimator  of  * * ...  the  k-statistic  k ....  ^ha 

r s rs 


krs  " U(*r*s  •••> 


rrf  (a-l)'.:- 

^ n\L  1 X2:  ...} 


ttCcp,  :)  '■(p2I) 


X , \ - 

t"(p,  i?,  2. 

„izh 


where  the  sums  and  products  are  taken  over  tha  same  partitions  and  indie 
as  in  Eq.  (351). 

To  illustrate , let  us  work  out  k^ ^ . The  partitions  of  3 are  3 
3 7 

and  L , and  of  2 they  are  2 and  l-.  We  shall  therefore  have  terms  in  [3 
[3l“],  [2“1],  [212],  and  [L^].  The  first  term  is 


C-n^oi.o:  3:.  2: 


3: .2:  r 3 . 2 1 r 32 1 

i:*t:  (3 : ) 1 (2 : ) 1 m[2^ 


This  is  followed  by 


and  by 


and 


(-n0+1o: 

■ 1 1 

3:2: 

[3.12] 

[312] 

i:  2i 

(31). (I!)2 

m[1+2] 

(-D1+0i: 

.01 

31.21 

[21.2]  3 

1111. 1! 

(2 : ) 1 (1: ) 1. 

(2 1 ) 1 M[2+11  1 

l term,  however,  consists 

of  two  parts,  name 

<-i>2+02: 

.01 

31.21 

[l3. 2]  _ 

2 r 2 13 1 

3!.  11 

(1 1 )3 .2 1 

m[3+1] 

c-n1+1u 

.1: 

31.2! 

[21. I2] 

3 r 2 13 ] 

1:11.2: 

2:11. (i : )2 

m[2+2^ 

" hM 

last  term 

is 

(- 1 ) 2+12 1 

.i: 

31.21 

[l3_.l2l 

2[15] 

3121 

(l!)3.  (1!) 

2 m[3+2] 

we  have 

that 

[11] 

mi  _ iiiiL 

. 3f22lJ 

5[2.1_3I  . 

2[15] 

M[2^  M 

[3] 

[3] 

Ml  J 

m[4] 

* 

M 


[3] 


_ 5f21J] 


Fortunately,  many  k-statistics  have  been  tabulated  in  terms  of 
augmented  symmetric  functions.  These  tables  are  arranged  according  to 
the  order  of  the  k-statistics  and  the  augmented  synsnetric  functions  in  the 
table.  The  order  of  a k-statistic  is  equal  to  the  sum  of  its  indices; 
the  order  of  an  augmented  symmetric  function  equals  its  weight  as  defined 
by  Eq.  (B29).  For  example  k^  is  of  order  3 and  k321  is  of  order  six. 


¥ 


The  first  order  table  is 


kj_  = [1]/M  . 


(B54) 


The  second  order  table  of  k-statistics  and  augmented  symmetric  functions 
are  shown  in  Table  B3. 

TABLE  B3  - k-statistics  and  augmented  symmetric 
functions  of  weight  2. 


kll 

k2 

ri2]/M™ 

I 

-1 

[2]/M 

1 

I 

The  conventions  for  reading  Table  B3  are  the  same  as  reading  Table  B2. 
Thus 

[12]/M[2]  = k„ 


[2]/M  = ku  + k2 

k11  = [12]/M[2] 
k2  = -[12]/M[2]  + [2]/M 


We  have  now  acquired  most  of  the  tools  needed  to  find  the 
unbiased  estimator  of  a moment.  Next,  we  will  state  an  algorithm  for 
using  these  tools  and  apply  the  algorithm  to  two  examples. 
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6.  Algorithm  Bl;  Finding  the  Unbiased  Estimator  of  a Population  Value 
At  this  point,  we  choose  to  outline  the  preferred  method  for 
calculating  the  unbiased  estimator  of  a population  value. 

a.  Find  the  population  value  in  terms  of  cumulant  products. 

If  the  population  value  is  in  terms  of  parent  moments,  this  can  be  done 
by  employing  Eqs.  (B48a) - (B48g) . 

b.  Convert  each  cumulant  product  to  its  corresponding  k- 
statistic,  in  accord  with  Eq.  (B53 a). 

c.  Using  the  k-statistic/augmented  symmetric  function  table  of 
the  same  order  as  the  k-statistic,  find  each  k-statistic  as  a sum  of 
augmented  symmetric  functions. 

d.  Using  the  augmented  syrranetric  function/power  sum  table  of 
order  r,  find  each  augmented  symmetric  function  as  a sum  of  power  sum 
products . 

e.  Replace  each  power  sum  (r)  by  a power  sum  about  the  mean 

s , defined  as 
r 

M 

sr  = £ (xi  - m[)r  = Mmr  . (B55) 

i=l 

This  substitution  is  justified  by  the  fact  that  the  value  of  any 
k-statistic  is  independent  of  the  choice  of  the  origin,  so  the  con- 
venient choice  is  one  such  that  m|  = 0 and  therefore  s^ = 0.  The  resulting 
unbiased  estimator  can  then  be  expressed  in  terms  of  sample  moments  about 


the  mean. 


yr 


Ik 


7.  Algorithm  B1  as  Applied  to  u? 

a.  Find  p,  in  terms  of  cumulants  by  Eq.  (B48a) , u2  = K2> 

b.  Find  the  k-statistic  corresponding  to  *2>  According  to 
Eq.  (B53a) , h?  = and  therefore  D*2  = k^. 

c.  Find  k?  in  terms  of  augmented  symmetric  functions  via  a 
2nd  order  table  such  as  Table  B3.  We  thus  find 

k2  = -[12]/M[2]  + [2]/M. 

Therefore 

4,  - k2  = [2]/M  - [12]/(M(M-D)  • 

d.  Find  the  augmented  symmetric  functions  in  terms  of  power 
sums.  From  Table  B2,  we  know  that 

[2]  = (2) 

[l2]  = -(2)  + (l)2  . 


The  unbiased  estimator  ^2  can  thus  be  expressed  in  terms  of  power  sums 


P2  = ^ ([2]  - [ l2 ] / (M- 1) ) 

= jj  t (2)  - ((l)2-  (2) ) / (M-l) ] 

= ^ C (2) [1  + 1/ (M-l) ] - (l)2/ (M-l)} 


e.  Replace  each  power  sum  by  a power  sum  about  the  mean  as 


defined  by  Eq.  (B55), 


However  because  s^  = 0,  this  simplifies  to 

If  M •) 

^2  = M S2  M-T  ^ * 
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f.  Express  the  unbiased  estimator  - in  terms  of  sample  moments 


about  the  mean.  That  is 


- _ _M 2 _ _M_ 

“ M-l  M ~ M-l  m2 


8.  Algorithm  Bl  as  Applied  to  M^. 


a.  Equation  (B48c)  enables  us  to  find  in  terms  of  cumulants  as 


^4  = + ^2  ‘ 


b.  Because  of  Eq.  (B53a),we  can  express  in  terms  of  k-statistics 


“4  " k4  + 3k22  ' 

c.  These  two  k-statistics  are  both  of  fourth  order;  we  can  find 
each  of  them  in  terms  of  augmented  symmetric  functions  via  a fourth  order 
table  such  as  Table  B4  below. 

TABLE  34  - k-statistics  and  augmented  symmetric  functions  of 
weight  4. 


r 


> . 


T 
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From  34  we  find  that 


k4  = -6[I4]/Mt4]  + 12 [212] /Mt3]  - 3[22]/M[2]  - 4[31]/Mt2] 


+ [4]/M 


k2?  - [14]/M[4]  - 2[212]/M[3J  + [22]/M[2]  . 


Therefore,  4^  can  be  expressed  in  terms  of  augmented  syranetric  functions 


as 


5 = 6_[  l4]  12  [212]  3 [2  2 ] 413 1J  J4.1 

4 m[4J  m[3]  m[2]  m[2]  m 


+ 3f  i£l  . H2±!i  + iri_  } 

1 M[4]  M[3]  J ' 


3[14].  | 6[212I  4_j3_l]  14J 

M ™ M[3]  Ml2]  M ' 


d.  Table  B5  gives  the  relation  between  augmented  symmetric 
functions  and  power  sums  of  weight  4. 


TABLE  B5  - Augmented  symnetric  functions  and  power  sums  of 
weight  4. 


In  Section  2 of  this  appendix,  we  learned  how  to  use  a table  of  this  kind 
for  finding  augmented  symmetric  functions  in  terms  of  power  sum  products. 
We  therefore  know  that  Table  B5  implies  that 


[l4]  = -6(4)  + 8(3)  (1)  + 3 (2) 2 - 6(2)  (l)2  + (l)4 

[212]  = 2(4)  - 2(3)  (1)  - (2) 2 + (2)(1)2 
[31]  = -(4)  + (3)  (1) 

[4]  = (4)  . 


We  can  now  find  in  terms  of  power  sum  products,  as 

= -3  [-6  (4)  + 8 (3)  (1)  + 3 (2)2  - 6 (2)  (l)2  + (1)4]/M[4] 
+ 6 [2  (4)  - 2 (3)  (1)  - (2) 2 + (2)(1)2]/M[3] 


4 [- (4)  + (3)(1)]/M[2]  + (4)/M 


= (4) 


+ (2) 


i_  + J^_  + 08_ 
I2]  mI3]  m[4]. 


(4) 


ri  + j 

Lm  m 

1 1 6 9 

L M[31  ' M^4^ 

ri  + j 
LM  M 


+ (3) (1) 


12 


24 


im(2]  h™  m[4] 


+ (2)  (1)" 


_6 18 

lm[3] 


.mi 


4_  + -JL_  + _iL- 

™ M<31  Ml41 


-4(3)  (1) 


-L-  + _L_+  _L_ 

lm[2]  m[3]  m[4] 


3 (2) 


2 3 

? 

1 3 1 

Lm[3]  +M[4^. 

+ 6(2)(1) 

1 

'“a 

+ 

r—% 

CO 

' ‘"x 

t 

3 (1) 


e.  The  expression  for  can  be  greatly  simplified  by  substituting 
power  sums  about  the  mean  for  the  power  sums;  without  loss  of  generality 
we  can  assume  that  s^  = 0 and  obtain  the  simplified  expression 


•r  f 4 


-ill!  ' 


y»«»  ■ 


— - ■ 


%•**  v 5^ 
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" 1 . 4 12  18 

o 2 

'2  3 

— 

iJj  — Q 

4 4 

1 m[3]  1 m[4]J 

' 3S2 

M^3^  M[4]J 

f.  In  terras  of  sample  moments  about  the  mean,  £,  now  appear 


^4  = ”4 


1 + 


12 


18 


(M"L)  (M-l)f2^  (M-l)[31 


- 3m„ 


2M 


3M 


L(M-I) 


[2] 


(M-l) 


[3] 


1 

I 

I 

I 
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9.  Algorithm  B2 : Another  Method  for  Evaluating  k-statistics  f 12 1 

We  have  already  shown  how  to  find  a k-statistic  of  order  m in 
terras  of  augmented  symnetric  functions  from  a table  of  order  m.  The 
purpose  of  this  section  is  to  show  how  to  evaluate  a k-statistic  of  order 
m,  given  a table  of  a larger  order  n.  This  method  is  outlined  below. 

a.  We  want  to  find  the  m-th  order  k-statistic  k ...  from  on 

rs 

n-th  order  table.  This  is  accomplished  by  first  finding  the  n-th  order 

k-statistic  k 1 (n-m)  in  terms  of  augmented  symmetric  functions. 

For  instance,  if  we  wanted  to  find  k?  from  a fourth  order  table,  we  would 

first  have  to  find  k or  equivalently  k 

Z 2lz 

b.  Notice  that  the  k-statistic  k . (n-m)  can  be  represented 

r s • • • i 

as  a sum  of  symnetric  functions  which  all  contain  unity  to  at  least  the 
(n-m)th  power.  For  example,  Table  34  implies  that  can  be  represented 

in  terms  of  augmented  symmetric  functions  as 

k2ii  = -[14]/Mt4]  + [212]/M[3J. 

In  this  case,  (n-m)  is  equal  to  2;  we  can  see  that  each  augmented  symmetric 
function  contains  unity  to  the  second  power.  It  is  also  true  that  the 
argument  u for  all  M^U|,s  in  the  expansion  for  k .(n-m)  is  as  least 


tfSMfl 
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as  large  as  (n-m).  For  the  case  of  k^j,  arguments  of  and  are 


both  obviously  greater  than  two. 

c.  In  order  to  find  k from  k 


(n-m),  we  simply  reduce 


rs . . . ' rs . . . 1 

all  powers  of  one  in  the  associated  augmented  symmetric  functions  and  all 
arguments  of  by  (n-m).  From  k^?J  we  could  then  find  that  k?  is 

kn(2-2)  = k2 


- -[L4-2]/M14-2)  + (21<2-2)]/m‘3-21 
.-[12)/M[2)  + (2]/M  , 


as  can  be  verified  from  Table  B3. 

We  can  prove  this  algorithm  by  considering  Eq.  (B53a),  the 

definition  of  the  k-statistic  k , 

rs. . . 


k = U (x.  k . . . ) . 

rs. . . r s 


(B56) 


It  follows  that 


krs...l(n-m)  = U(Vs  > 


(B57) 


Because  we  also  know  that 


krs...l(n"m)  = U(Vs-  '^(n'm)) 


(B58) 


From  Eq.  (B51),  however. 


Vs"^l 


, (n-m)  _ 


5 •(.  1 \ /M1  \ 

\£(P  1)  „!  /„_i  ill  tr  I Pi  ■ *-  ■ ' 2 


••s'  TTJ  2l\  III  i u.(n 

2_i  t C7'- 1*^2*  -••/a  _ ? I * / \ P2  • j 1 


-m) 


(B59) 
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We  can  use  Eq.  (B30)  to  find  k (n-m)  from  Eq.  (B53b).  That  is, 

C S • • • 1 


krs...l(n'ra) 


(-i)Z(p~L)  tuco-d:] 

. ttU1:x2:  ...} 


r__r  f ^1  *2  , (n-m) 

ITT(P1  P,  •••! 


n[(PlI)  1(p2:)  2...}  M[i:p+(n'm)] 


However,  we  know  from  Eq.  (B53b)  that  k can  be  represented  as 

rs. . . 


(B60) 


rs. 


E(P-D 


rU=-UL). 


r!s! 


N 

[TT  (Pl  P2  ...)}] 


n W ...}  f x2  , JZp] 

1 2 TT  {(Pll)  i(p2!)  -...}  M 


Notice  that  the  coefficient  corresponding  to  each  term 
[rr(p  “p  2...)]/M[Z0]  in  k is  identical  to  the  coefficient  corre- 

• — — r S . . . 

sponding  to  each  term  [tt  (p  ^p  ‘ 2 (n-ra)  ] /M^Zp  + ^n_ra^  in  k (n-m)  . 

I Z ) l rs. . . I 

For  example,  the  coefficient  corresponding  to  each  term  [2]/M^  in  k„,  is 

i. 

2 f3] 

identical  to  the  coefficient  corresponding  to  each  term  [21  ] /M  in 
k211‘  This,  of  course,  is  the  same  principle  used  in  the  algorithm  out- 
lined above. 


10.  Partitions  as  Employed  in  Defining  Cumulants  and  k-statistics 
We  have  already  noted  in  Eq.  (B45),  that 


K = r ; \ (~I)P"1(P-l)'[M,pl  ' ^ 
r _ xL:x2:...  \ pLi  , 


X.  /Uf'  \ X 


p"r2 


_2  j 

P2 


where  the  sum  is  over  the  partition  of  the  integer  r such  that 

Pl^l  + + • • • = r anc*  ^"l  + ^2  + * ’ * = p * ^ examP-'-e  clarify  the 


term  partition.  Consider  the  integer  4;  the  set  of  combinations  belonging 


to  the  partition  of  4 according  to  the  above  prescription  is 


{4 , 3 + 1 , 2+2,  2 + 1 + 1,  1+1  + 1+ l}  . 


This  partition  can  be  represented  more  elegantly  as 


(4,  3+1,  2x2,  2+1x2,  1x4}  , 


or  better  yet  as 


r / 1 -,1,1  ^2  0 1 ,2  1^1 

14,  3 1,  2,21,1], 


Every  positive  integer  has  a corresponding  partition.  One  possible  combi- 


nation belonging  to  the  partition  of  r is 


X1  X2  Xt 

P1  p2  * * * pt 


where,  of  course, 


?1XL  + P2X2  + •”  Pt^t  = r * 


We  could  be  more  specific  and  consider  the  j-th  combination  of  the  partition 


of  r: 


Xlj  X2j 
Plj  P2j 


'lPlj'1J  • 

1 J 

The  quantities  p_  and  \ respectively  will  be  called  the  i-th  base  and 
the  i-th  power  of  the  j-th  combination  of  the  partition  of  r.  To  each 
j-th  combination  there  corresponds  a sum  of  its  powers  p : 


Pj  = \XJ  + x2J  + ...  xtj 


* 
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or 


= I X.  . 

i LJ 


Using  our  new  definitions,  can  be  represented  as 


A = 

r 


p -1 

^rl(-l)  J (p --1) : 
, X . . 

3 "Up,,): 


n .> 
i ii 


or  equivalently  as 


X. 


h = C.  tt  (p.1  ) 

r — , j ..  p . . 

j 1 ^ 


1J 


where 


C.  = 
J 


Pj"1 

r: (-1)  J (p 


jf<V: 


As  an  example,  we  compute  ^ using  Eq.  (B61),  by  first  noticing  that  the 
partition  of  3 is 

{3 , 2 + 1,  1 + 1 + 1} 


or 


£3 1 , 21!1,  I3} 


The  first  term  of  is  then 

Pi'1 

X 3i  (-1)  1 (p  -1)1 
Ci  , 


TT  (u.< 


rj  t(pn)'  (Xtl)!} 


(HI  ) 

pil 


’il 


3u.-Dlla.-ui. 

(3:>  (1!) 


= U ' 
3 


(B61) 


(B62a) 


(B62b) 
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This  is  followed  by  the  expression  for  the  term  corresponding  to  the  com 
bination  2^1^  and  j = 2 : 


c?  n (U'  ) 

2 i pi2 


p2-^ 

\.2  3i  (-1)  (p2-i> : 


X, 


n i(Pi2)l  i2(X.2)l}  1 


" ^2> 


Xi2 


= 3J,.(i1)2_'112-i,i: 

[(2!)1(i:)K(i:)1(i:)}  2 1 

= -3U^[  • 


The  third  term,  which  corresponds  to  1 , is 


tt  (4  ) 

3 i pi3 


p3“i 

x..  3 : (-i)  (p,-i>:  x. 

TT  (ji  ' ) 


n {(pi3).'"i3(Xi3)i}  1 Pi3 


iU-l)(3~I}  (3-1)1  ^t)3 

£ (11  )3  C3! )}  1 


= 24 


,3 


We  add  these  3 expressions  together  and  find  that 


= 4, 


- 


+ 24 


,3 


In  Eq.  (B61)  we  represented  a single  cumulant  in  terms  of  parent  moments 
When  dealing  with  cumulant  products,  however,  more  than  one  partition  is 
is  required.  For  example,  in  order  to  express  x3*2  in  terms  of  parent 


(B63) 
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£2,  1 + 1} 


or 


{21,  l2} 


The  term  corresponding  to  the  combination  2 is 


d-1) 


(2 : ) (1: ) 


= n 


2 • 


2 . 


The  term  corresponding  to  1 is 


m.-uil'.'-Vi).;  2 

(i:y (2D 


= -v>. 


Therefore 


*2  '“2  -Ui2 


From  Eqs.  (B63),  and  (B67),  we  know  that 


*3*2  = - 2U  'f{  + 2^[3)(U'  - ^2) 


= a 'a  ' - ■ ■ ■ • 2 - 1.  ' ' J-  'it  '■*  _ T..  '5 


H^J-2  '^3^1  ‘ ^2  M'i  + 5M'2P'i  ‘ 2Ui‘ 


(B6  7) 


(B68) 


Now  that  we  know  a^a 2 in  terms  of  parent  moments,  we  can  use  Eqs.  (B30) 


and  (B53a)  to  find  the  k-statistic  k^2 . That  is, 


k32  = 


U<H^)  - U<H^[2)  - 3U(JJ!^2y.[)  + 5U(^^[3)  - 2U^[5) 

[32]  [3 12]  3[22L1  . 5 f 2 13 1 2fl5] 

m[2]  ' M f 3 ] ' m[3]  T m[4]  ' m[5]  * 


(B69) 


In  order  to  express  a given  k-statistic  in  closed  form,  we  must  first 
change  Eqs.  (B66a)  and  (B66b)  from  a product  of  sums  to  a sum  of  products. 
In  order  to  accomplish  this,  we  must  change  our  indexing  system  somewhat. 
For  example,  it  is  clear  from  Eqs.  (B64b)  and  (B66b)  that  we  could  repre- 


sent and  ^2  as 


X3  C 11^*  3 + C21ii2M'l  + C31^1 

*2  = C12*2  + C2^i2 


and  ^2^2  3S 


*3*2  = cuc 


12^3^2  + C11C22P'3^1 


+ C21C12M'2P‘l  *U2  + C21 C22^  2^*  1*^  1 
+ C31C12y*  1 ‘^2  + C31C22^  1 '^1 


(B70) 


In  Eq.  (B70),  the  first  index  of  C^m  labels  the  j-th  combination  of  the 
m-th  partition.  We  could,  however,  allow  j to  lable  the  j-th  combination 
pair  of  the  partitions  3 and  2.  The  partitions  of  3 and  2 are 


(3,  2 + 1,  l + l + l}  and  [2,  1 + 1} 
{3L,  21!1,  I3}  and  [2l,  l2}  . 


or 


The  corresponding  set  of  combination  pairs  is 

£3L*2L,  3^- l2 , 2111»21,  2^1^* l2 , 13.2L,  13-12}, 

If  j is  the  index  for  combination  pairs,  would  be  represented  as 


*3*2  C11C12M'3’M,2  + C21C22^3*M,1 


. , , 2 


+ C31C32M'^'i'|J’2  + C41C42y’2Uru,[2 


! I 1 ^ , I , . _ , , 1 3 , 1 2 


+ C51C52?1  -M,2  + C6 1C62^  1 ,U,i 


For  a product  of  three  cumulants,  the  equivalent  sum  could  be  over 
combination  triplets;  in  general,  a cumulant  product  could  be  expressed 
as  a sum  over  a set  of  combination  multiplets.  Using  this  new  notation, 
any  cumulant  product  can  be  expressed  as 


X X 

rl  r2. 


! 


X 


= > j C . TT  <Jj,  ’ ) LJm 

j m \ 1 ijm 


\ 


Because  multiplication  is  associative 


X X 

rl  r2... 


I 


j \ ra 


C . 
jm 


II  )Xi> 

m,i  rijm 


K K 

rl  r2 . . 


V/—  r -1)! 

Zi  li-2 — I i 


n 


x. 


17  Up,j:  ljma ):}/  m,i>  ‘ljm/ 


ijm 


ijm 


• M. ' 

1 P< 


.X.  . 

\ ijm 


(B71) 


(B72) 


(B73) 


or  employing  Wishart's  notation  as  in  Eq.  (B51),  Section  5 


i 

I 
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s2(P-Dt 


f'j'  I 1. 

hh  = r's’  ^ (-irVK'*'n  t(o-l):}  77  J P1  J2 

rs...  r-s-v  nT(W  ...>j  !•  p2: 


P,  ^ 


From  Eq.  (B73),  we  can  easily  find  that  the  k-statistic  corresponding  to 
is 


rl  r2  ’ * ’ 


r ir2 ' 


— r i(-l)Pjm  (p.  -1)1  \ [n.  p..  ^ijm] 

\ i ii  m r ,jm  I m,i  i]m 


ki L1E_ 


m X- • / [ E p.  ] 

\ “tt  t(p..  )•  1Jm(x.  . )•}  / M » Jm 
V -•  rnm  wnm'  J 


t ljm  ijm 

because  we  know  from  Eq.  (B30)  that 


/ - — \ . \ " . ^ijm 

u i n o*1  ) ljra  = ■ q*-1-  ■ 

5 m>i  Pijm  / M[i  p jm 


Using  Wishart's  notation,  the  k statistic  k can  be  expressed  as 

r S • • « 


rs. , 


yc-i)Z(p'1}  TTj(o-l):} 

i ™ • • • • ^ 


IT  * S * 

X • O • . » i 


X1  X2 
l^(Pl  P2  ■ 


)] 


r X1  X2 

" C(PX!)  (P20...}  M 


1 

T 

X 

I 


Let  us,  for  example  consider  computing  k^  > t*le  unbiased 
estimator  of  *^2 » terms  of  augmented  symmetric  functions.  We  apply 
Eq.  (B74)  to  this  problem  and  let  r^  = 3 and  r2=2.  Therefore 


32 


Pjl'l 

31  (-D  J (Pn-1)! 


P i2-l 
2.'(-l)  ^ 


n n C(P1J2):''1J2(XiJ2)l}  M1'11 ' "J2' 


^ X . .. 

[TT  p..,  ii\.  , ^2] 
i Hijl  Pij2  J 

[p.,  +p7j 


(B74) 


(B75) 


We  recall  that  the  partitions  of  3 and  2 are 


(31,  21!1,  I3}  and  U\  1*}  , 
and  the  set  of  combination  pairs  that  Eq.(B75)  is  summed  over  is  therefore 

U1^1,  31.!2,  zhl-2lt  21!1.!2,  l3-2l,  13-12}  . 


1 ,2- 


The  term  corresponding  to  3 *2  is 


Pi i "1  P 12~2 

3i(-D  11  (pu-i)i  2!(-i)  (p12-p: 


[ pill  Pil2  ] 
i 


1 t(PtU>^1U(Xl11):1  i t(pil2);Xil2("il2):l  M 


[P11  + P 12^ 


C-1)(1~1)C 
(3 1 ) 1 (11) 


: (-n  (I~:L)' 

(2 1 ) 1 (10 


■1)1  f3  » 2 ] 


[1  + 1] 


= [32  ] /M1 


..  „ ,1.2  . 

The  term  corresponding  to  J l is 


P^i*l  P22~^ 

3 1 (-1)  <P21-1)1  21  (-1)  CP22~L>1 

'V  ^ • o ? 

— : , \ i i2 1 ,,  -vii  rr  f /»  i i ■LZ,i  t\ 


Xi21  \22 

[ 1 Pi  9 1 ?i 9? 


. pi21  ‘ i22  ] 


Up,9i>:  i21(^i9i>:}  " M 


tp21  + p22^ 


-i^q 

(31)1(l!) 


(102  (21) 


,1  2. 
> ‘1  J 
,[1  + 2] 


3 


The  term  corresponding  to  is 


The  term  for  the  combination  pair  1 '2  is 


The  last  term,  which  corresponds  to  the  combination  pair  1 -1 
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I 


We  can  define  ?(t)  a log  0(t)  and  make  the  change  of  variables  x=it,  to 


find  that 


+ 'r  = (^)r  exp  [p  (x)  ] 


\ - <s>r  »<*>| 


(B76a) 


(B76b) 


We  would  expect  then  that  G (t,l)  can  be  expressed  in  terms  of  F (t,l) 
in  the  same  way  that  is  expressed  in  terms  of  For  instance,  from 
the  fact  that  can  be  expressed  in  terms  of  cumulants  as 


= H4+4K3H1  + 3h2+6^2^  + 


we  know  that 


G(4)(l,t)  = F(4)(l,t)  + 4F(3)(l,t)F(1)(l,t)  + 3[F(2)(l,t)]2 


+ 6F(2)(l,t)[F(1)(l,t)]2  + [F(1)(l,t)]4  , 


as  is  confirmed  by  Eq.  (Al4d).  Therefore,  the  method  used  in  expressing 

(r  ) 

moments  in  terms  of  cumulants  can  also  be  used  to  express  G (l,t)  in 
terms  of  F ^ (l,t) . 
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APPENDIX  C 

VARIANCES  AND  UNBIASED  ESTIMATORS  OF  at  AND  3t 

1.  Introduction 

The  purpose  of  this  section  is  to  give  detailed  derivations  of 

quantities  mentioned  in  the  main  body  of  the  paper.  Here  we  derive  the 

A A 

unbiased  estimators  of  at  and  pt,  the  variances  of  at  and  pt,  and  the 
unbiased  estimators  of  these  variances. 

2.  Unbiased  Estimators  of  at  and  3t 

Equation  (36a)  gives  at  in  terms  of  parent  moments 

at  = i[4^2  + 3(p2)2  - n4]  . 

Algorithm  Bl  given  in  Section  6 of  Appendix  B shows  how  to  find  the  unbiased 
estimator  of  a population  value  such  as  at.  First,  we  find  at  in  terms  of 
cumulant  products.  Equations  (B48a)  and  (348c)  give  the  parent  moments 
1^2  and  p>4  in  terms  of  cumulant  s: 

2 2 

P4  = h4  + 3*2  . 

Therefore,  at  can  be  expressed  as 

at  = i[4(*2)  + 3(*2)2  - (k4  + 3*2)] 

= \[^2  - k4]  . (Cl) 

The  unbiased  estimator  of  at  is  then,  by  Eq.  (B53a), 

A lr/  A A 
at  = -[4  h2  - h4] 

= |-[4k2  - k4]  . (C2 ) 
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Because  k2  and  are  of  orders  2 and  4 respectively,  tables  of  orders  2 
and  4 are  necessary  for  finding  the  k-statistics  in  terms  of  augmented 
symmetric  functions.  Using  Tables  B3  and  B4,  we  find  that 


- (- 


. 1I21L  . MlLI  + iil  ^ 
[2]  [2]  m 'j 

Ml  J Ml  J 


From  the  second  and  fourth  order  augmented  symmetric  function/power  sum 
Tables  B2  and  B5,  and  Algorithm  Bl,  we  see  that 


” ' 6 1 


<"S2)  + M (S2)  + ('6S4  + 3s2  > 


^3T  <2S4  ■ s2>  + ('s4  + s2>  + (-S4>  ‘ IT  5 


-If  i, 

' 6 m[2)  m " h'4J  m[4]  ‘ m[3!  m13!  ' Ml2)  ' M[2>'  M 1 
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Similarly,  j3t  can  be  found  in  terms  of  the  sample  moments.  From  Eq.  (36b), 
we  know  that 

Bt  ■ 2l  [“4  - ^ ' 


In  terms  of  cumulants,  3t  is 


^ Ih4  + 3^  - 3*2  . k2] 


Again  we  use  Eq.  (53a)  to  find  the  expectation  values  of  the  cumulants 


We  already  know  k,  and  k in  terms  of  augmented  symmetric  functions  from 


calculating  at,  so  that 


Tables  B2  and  B5,  and  algorithm  Bl  show  that  j3t  can  be  expressed  as 


1 


i 

T 

1 


i i 


i 

i 


39 


A 


3t  =24  »• 


6(-6s4  + 3s-)  _i2_  2 

+ r 1 1 (2  S , S ^ ) 


M 


[4] 


M 


[3] 


M 


3 , 2. 

T2l('S4  + S2> 


4(“s4)  (s4)  (-s2)  (s2) 

“t  w "t  rol  ” , . } 


M 


f2] 


M 


M 


[2] 


M 


.if  fix  _i_  + _J_  + 24 

” OA  >•  Vw  + r ~ i T-  r^i  T 


36  ) s. 


24  M m[2]  m12]  m[3]  m[4] 


+ ( 


12 


[2]  ' J3]  - J4] 


M 

1 r ,1 


M 


M 


18  . 2 . 1 . 1 , 

} S2  ‘ ( [2]  + M S2 


M 


24 


36 


"24  1 (M  + „[2]  + v,[3 ] + [4]  ' ‘““4 


) Mm, 


M 


M 


M1- 


3 . 12  . 18  , ^2  ,1  . 1 


m[2)  ' Mm  T M[41  > H‘,I2  - < m[2]  + M > *2  5 


_1_  f , . , 7 24  36 

= 24  + (M-l)  + -v  [ 2 ] + ,„_1X[3]  } m4 


(M-l)1 


(M-l) 


14  6 

3M(  -rr-^r  + + 


> 2 _M_  •, 


<■«-»  ' (M-l)  I21  (M-l)131  2 O'"1’  2 


A A 

We  can  also  express  at  and  3t,  given  in  Eqs.  (C4)  and  (C7) , as 


A 

2 

at 

= R 

+ 

S ITU 

+ T m 

a 2 

a 2 

a ‘ 

A 

2 

Pt 

= R3m2 

+ 

S3m2 

+ Tgin 

where 


2M 


a 3 (M-l) 


S = ? ( -^L“  + 4 


a 2 v (M-l)  (m-1)[2^  (M-l)[3^ 


T = .1  (1  + -I—  + —Ik.... 

Ta  6 (L  + (M-l)  /M_i  >.  [2] 


36 


(M-l)1 


(M-l) 


[3] 


■ 'tsanki 


(C7) 


(C8a) 

(C8b) 


(C9a) 


(C9b) 


(C9c) 
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A A 

3.  Variances  of  at  and  8 1 


Equation  (51a)  gives  the  variance 

A A 2> 

var  at  = {(at)  ) - (at) 


We  already  know  at  from  Eq.  (C2);  the  squar 


It  is  not  immediately  obvious  how  we  can  fi 
A 2 

(at)  from  Eq.  (Cll).  However,  the  k-statis 
can  be  reduced  to  a more  pliable  state  by  r 
a sum  of  individual  k-statistics.  Consulti 
Appendix  F,  we  find  that  [11] 


2 k8  . 16k62  . 48k53  . (M+33' 


R?  = ' 24 (M-l) 


£ (1  + 


. M t + * 

8 (M-l)[2] 


We  can  now  write 


Performing  the  proper  multiplications,  it  appears  that 


144  M 


By  Eqs.  (B53a)  and  (42),  however 


According  to  Eqs.  (B49a) - (B49g) , the  odd  cumulants  involve  moment  products 


that  must  include  at  least  one  odd  moment  as  a factor.  Because  all  odd 


parent  moments  as  defined  by  Eq.  (29)  are  zero,  odd  cumulants  are  also 


we  know  that 
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2 _ A 
7at  = var  at 


-i-  [ Hi  + — V2+  -01+33)  2 72M 

*»  m i \ ' /w  1 \ k . *r 


36  1 M (M-l)  (M-l)  4 


(M-l) 


[2]  4 2 


24M(M+1)  4 _ __6  8 (M  + 7 ) , K 
(M-l)[3]  2 M <M“D  *2 


. 16K4  . 16(M+1)  2 i . .2 

+ M + (M-l)  K2  } - ‘ 


(C15) 


Equation  (51b)  gives  the  variances  of  gt  as 

A A 9 9 

var  Pt  = < (Pt)  > - (pt)  • 

, A 

We  already  know  Pt  from  Eq.  (C6);  the  square  of  this  quantity  is 

«V  - 5^  [kj  - 2k4k2  + k*  ] . 

Substituting  Eqs.  (C12a) - (C12 c)  into  Eq.  (C16),  we  find  that 


(C16) 


(pt)2  = -1-  { Hi  + -16  k + 48k53  (M  + 33)  72M 

576  1 M (M-l)  k62  (M-l)  (M-l)  k44  + [2]  k422 

(M-l) 

+ -1-44M  ,,  , 24M(M+1)  2k6  2(M  + 7) 

m r i - TT”  - rr — k . _ 


(M-l) 


[2]  332  (M-l)[3]  2222  ’ ~ ' (M-l)  k42 


-12-  k + — , -QL+D  k , 

(M-l)  *33  + M (M-l)  k22  1 


By  Eqs  (B53a)  and  (51b),  the  variance  of  Pt  is  then 


(C17) 


mm 


f 

. I 
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2 A 

a 3 var  gt 
Pc 


±ri+  16H6H2  (M_t3_3)_  k2 

576  1 M (M-l)  (M-l)  4 


72M 


(M-l) 


[2] 


2 

v2 


. 24MCM+1)  4 ^6  2 (M+  7)  ^4 

(m-d[3]  2 M (M'1}  4 2 M 


(C18) 


A A 

4.  First  Order  Approximations  for  var  at  and  var  ?t 


A 


From  Eqs.  (Cl)  and  (C15)  we  know  that  the  variance  of  at  can 
be  written  in  terms  of  cumulants  as 

2 


A 

var  at 


1,  J3  , 16k6*2  (M+33)  2 
36  L M (M-l)  (M-l)  4 


72M 

(M-l) 


[2]  *4*2 


. 24M(M+1)  4 f^6  8 (M  + 7 ) ^4  16(M+1)  2 . 

+ (m.1}[3]  K2  " M ' (M-l)  V2+  M + (M-l)  *2  J 


- I i (4h2  - »4)]2 


1 Kg  16k6k2  (M+  33)  2 . 72M  2 

= 36  [ M + (M-l)  + (M-l)  *4  + (M_l)  [2]  *4*2 

* 24M  (M  +1)  4 f^6  8 (M  + 7 ) 1^4  16.ffl+  1)  2 

+ /m1,[3]  *2  " M " (M-l)  H4K2  + M + (M-l)  2 

(M-l) 


- 16k.2  + ®K2H4  “ k4  3 • 


u 

0 
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To  first  order  in  1/M,  var  at  is  then 


^ L 2 2 4 

var  at  =.  (Kg  + 16k^k2  + 33k4  + 72k^k7  + 24k2  - 8k ^ 


- 56k. K.  + 16 h + 16k^  ) 

4 2 4 2 


We  can  substitute  Eqs.  (60a)- (60d)  into  Eq.  (C19)  and  find 


(C19) 


A 


var  at  in  terms  of  parent  moments: 


Var^-3k  [^8 


2a^2  - 35M.J  + 420^4M-2  - 630^2  ) 


+ 16  (^  - 15h4u,2  + 3 Op.8 ) p2  + 33  (^  - 3P2)2 

+ ?2(^4  - 3p2)  ^22  + 24^2  - 8fcfi  - 15n^2  + 30H8  ) 

- 56  (M-4  - 3P2 ) p2  + 16(H4  - 3p2)  + 16P*] 

■ 355  l“s  - 28V2  ■ 3*l  + 42(V2  ' 63l>12 


+ 16HgM-2 


240^2  + 480^2 


+ 33^  2 - 198^2  + 297^2+  72^2  - 216^ 

+ 24^2  - 8pfi  + 120^  - 24(^2 

- 56m-4|J.2  + 168p8  + 16p-4  - 48M-2  + 16p2  J 

= 36M^  ^8  " 12m‘6M'2  " 2y,4  ' 18*VJ,2  + 171^2  + 136m,4^2 

- 72M-2  - 8h6  + 16p>4  - 32p*  ] . 


(C20) 


H 


96 


r 


+ 2^2  - 2(h6  - 15m-4^2  + 30m3) 


- 14(n4  - 3^2  + (H4  - 3^2)  + 4 2 J 


5^5  '“8  - 2V2  ' 35>\  + 42<Vz  ' 630i2 


+ 16^2  - 240m. 4M-2  + 480m 2 


+ 33M4  - 198m 4M 2 + 2 9 7m 2 + 72m4M*  - 216m2 


4 3 

+ 24m2  - 2m6  + 30m4M2  - 60m 2 


14Va2  + 42P<2  + ^4  ' 3u2  + ^2 ^ 


'sTS  ‘“a  - l2V2  - 24  4 54^2  - 4 34 2 


- 2Mg  + 16m4M>2  - 18m2  + M<4  ■ 2m2  ]• 


(C22) 


A A 

5.  Unbiased  Estimators  of  var  at  and  var  3t 


Equation  (53a)  gives  the  unbiased  estimator  of  the  variance  of 


A 

at  as 


U (var  at)  = (at)2  - U[(at)2] 


From  Eq.  (Cl),  we  know  that 


at  = - [4k2  - k4] 


(«)2  = i U6H2  - SVl  + *2J  . 


(C23) 


Using  Eq.  (B53a)  to  derive  the  unbiased  estimator  of  (at)  , we  find  that 


U[(ot)2]  = ± [16k22  - 8k42  + k44]  . 


(C24) 


.1 a. 


■ a m «--y*  '•  *r  *4*  *• 

V-  . m..  •%  ■ 


i 

1 

I 


""" 
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Now  we  need  only  square  the  right  hand  side  of  Eq.  (C8a)  to  obtain  (at)  : 
(<£t)2  = (Ram2)2  + 2(R^)(S^  + T„mA)  + (S^m2  + T^)2 


a 2 a 2 a 4' 


a 2 a 4 


2 2 3 2 4 

R m + 2R  S m + 2R  T m.  m + S m„ 
a2  aa2  aa42  a2 


2 2 2 

+ 2S  T m. m + T m.  . 
a a 4 2 a 4 


(C25) 


T* 


I 

I 


The  coefficients  R^,  S^,  and  are  given  by  Eqs.  (C9a)-(C9c).  Me 
combine  Eqs.  (53a),  (C24),  and  (C25)  and  obtain 


/Y 


5 U (var  at) 


«• 

2 2 3 

„„  4 

= R m_  + 2R  S m„ 

+ 2R  T n. ra_  + S m~  + 

__ 

a 2 a a 2 

a a 4 2 a 2 

2 2 4 , 

2 1 

• • 

+ V>4  * 9 k22  + 

9 *42  " 36  44  * 

5aTam4m2 


(C26 ) 


A 


Similarly,  the  unbiased  estimator  of  var  0t  can  be  found  from  Eq.  (53b): 

A A 2 2 

U (var  3t)  = (0t)  - U[(Pt)  ] 

It  follows  from  Eq.  (C5)  that 

3t  = A [H4  " H2] 


(Pt)  ”575  tK^  “ ^K4H2  + * 2 ^ * 


(C27) 

2 

Using  Eq.  (B53a)  to  derive  the  unbiased  estimator  of  (j5t)  , we  note  that 

(C23) 


U[(Pt)2]  = [k44  - 2k 


42  + k22]  * 


A 2 


Squaring  the  right  hand  side  of  Eq.  (C8b)  it  follows  that  (pt)  is 


A 


,2  2 


-.2  4 


(8t)  = Rgm2  + 2R?S3m2  + 

+ 2SgTpm4m2  + l|m2  . 


(C29) 
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The  coefficients  Rg , Sg,  and  Tg  are  given  by  Eqs . (ClOa)-ClOc) . We  combine 
Eqs.  (53b),  (C28) , and  (C29),  to  find 

/^\  A 

Cfgt  = U(var  8t) 


= ^2  + 2R8S8m2  + 2R8T8m4m2  + Slm2  + 2S8TPm4m2 

2 2 1 k42  k44 

T3m4  ' 576  22  288  " 576 


(C30) 


6.  Evaluation  of  the  k-statistics  ,_k^2 , and 

2 

Equations  (C26)  and  (C30)  give  the  unbiased  estimators  of  c 

2 

and  a in  terms  of  sample  moments  and  k-statistics.  Here  we  find  the 
Pc 

relevant  k-statistics  in  terms  of  sample  moments  also,  so  that  the 
unbiased  estimators  of  the  variances  can  be  easily  evaluated. 

First  we  calculate  the  k-statistic  k ^ • This  k-statistic  is  of 
order  4;  therefore  we  can  find  it  in  terms  of  augmented  symmetric 
functions  via  Table  B4,  which  is  also  of  fourth  order.  That  is, 


k _ nh. . 21212]  i£ i 

*****  **  r / 1 Ml  — r 1 • 


22  m[3] 


M 


[2] 


(C31) 


Applying  Algorithm  Bl  to  Eq.  (C31)  and  Table  B5,  we  arrive  at 
k,,  = (-6s,  + 3s2)/M[4]  - 2 (2s.  - s?)/M[3]  + (-s.  + s2)/m[2] 


- 2.. 2 , 1 _2 . 3 . M , 1 ^ 4 ^ 6 ..  (C32) 

- m2M  ( pj  + pj  + . . ) -ra^M  ( p,  + p,  + ) 

ML  J Ml  J Ml  J Ml  J Ml  J J 


We  can  express  the  sixth  order  k-statistic  k^2  in  terms  of 
augmented  symmetric  functions  by  consulting  the  sixth  order  table  in 
Appendix  E and  find  that 


100 


, 540  324  75  12  _4 9 _J _J 

2‘  Ml«l  m!51  m[M  + m[3]  m[3]  + J3]  „[2]  > 


M 


M J M 


M 


M 


M ML 


2 240  144  30  8 4 

3V61  Mi5]  k[41  m'4]  «!3)) 


3,  90 

+ s_(  - 


54 


15 


2 ' m[6]  ' m(5]  ' m[4]  m[3] 


m < - L-  16  114  432  720  . 

= 6 M[2]  ' M[3]  ' M[4]  * M[5]  ' m[6] 


14  . 87 


324  , 540 


m4m2M"(  ^,[2]  + J3]  + u[4]  + „[5]  + „[6 ] 


M 


M 


M 


♦ m2^  4 


, 38  144  240 

i r / l ~i~  r i-  i r s 1 ) 


3“  V3]  m[4]  \m[5]  m[6] 


^ 3 3 _3 15  54  90 

' ^ ( ' -J31  ‘ m[4]  ’ m[5]  • m[6] 


M 


M 


(C33) 


The  k-statistic  is  of  eighth  order.  In  Appendix  E we  have 
a twelth  order  table;  an  eighth  order  table  is  not  available.  We  there- 
fore use  Algorithm  32  from  Section  9 of  Appendix  B to  find  k, , in  terms 
of  augmented  symmetric  functions.  This  procedure  gives 


k44  = 36[18]/M[81  - 144[216]/M[71  + 180 [22 l4] /M[6 ] - 72 [2312] /M[5] 
+ 9[24]/M[4]  + 48(315]/M[6]  - 96[3213]/M[5^  + 24 [3221] /Mt4] 

+ 16[3212]/M[41  - 12[414]/'rt[5lf24[4212]/M[4]  - 6[422]/M[3^ 

- 8[431]/M[3]  + [42]/M[2]  . 


Applying  Algorithm  B1  to  the  augmented  symmetric  function/power  sum  tables 
of  weight  8 given  in  Appendix  E,  yields 


* f -■r  +.  <*  * 

r ifi  \M  i jflftMH 


36  (-5040s Q + 3360s.  s + ?688scsQ  + 1260s?  - 1260s. at  - 1120s:s_ 
8 62  53  4 42  32 


+ 105s2)/M[8] 


- 144  (720s8  -480s6s2  - 384s5s3  - 180s2  + 180s4s2  + 160s2s2  - 15s2)/M[7] 


+ 180(-120sg  + 84s6s2+64s5s3  + 30s^  - 33s4s2  - 28s2s2  +3s2)/M[6] 


2.-  -2  • '-2«  ..A/m [5] 


72  (24sg  - 20s6s2  - 12s5s3  - 6s4  + 9s4s2  +6s3s2  - s2)/M 


2 24  f41 

+ 9 ("6sg  + 8s6s2  + 3s4  - 6s4s2  + s2)/Ml  J 


+ 48(-120sg  + 60s6s2+64s5s3  + 30s4-  15s4s2  - 20s2s2)/M[6] 


96(24sg  - 12s6s2  - 14s5s3  - 6s2  + 3s4s2  + 5s2s2)/M[5] 


+ 24(-6s6  + 4SgS2  + 4S5S3  + s2  - s4s2  - 2s2s2)/M^ 


+ 16  ("6Sg  + SgS2  + ^s5s3  + 2s2  - s2s  )/M^ 


- 12 (24sg  - 12s6s2  - 8s5s3  - 6s2  + 3s4s2  )/M[5] 


+ 24(-6sg  + 3s6s2  + 2s3s3  + 2s2  - s4s2)/M^ 


‘ 6<2s8  ‘ 2s6S2  ' 3 4 + S4S2)/m[3] 


- 8(2Sg  - s5s3  - s4)/M[3] 


+ ('s8  + s4>/mI2] 


=*  (-181440s.  + 120960s, s + 96 768s. s.  + 45360s2  - 45360s.  s2 
8 62  53  4 42 

- 40320s2s2  + 3780s2)/M[8] 

+ (-103680s.  + 69120s, s_  + 55296s, s.  + 25920s2  - 25920s. s2 
8 62  53  4 42 

-23040s2s  + 2160s4)/Mt7] 


102 


+ (-2 1600s  + 15120s  s + 11520scs_  + 5400s.2  - 5940s. s2  - 5040sJso 
8 62  53  4 42  32 

+ 540s4  - 5760s,  + 2880s£s,  + 3072s, s,  + 1440s2  - 720s. s 2 
2 o 62  5 3 4 42 

- 960S2s2)/m[6J 

+ (- 1728sQ  + 1440s, s,  + 864s, s,  + 432s2  - 648s, si  - 432s2s„ 

o 62  53  4 42  32 

+ 72s4  - 2304s  + 1152s, s,  + 1344s, s,  + 576s2  - 288s. s2 

2 8 62  53  4 42 

- 480s2s  - 288s,  + 144s,  s,  + 96s,s,  + 72s.2  - 36s. s2)/M[5J 

32  o 62  53  4 4 2 

+ (-54s,  + 72s, s + 27s2  - 54s. s2  + 9s4  - 144s,  + 96s, s,  + 96s,s, 

8 62  4 42  2 8 62  53 

+ 24s2  - 24s, s2  - 48s2s  - 96sQ  + 16s, s,  + 64s, s + 32sf 

4 42  32  8 62  53  4 

2 22  F41 

- ^s^s^  - 144sg  + ^^S6S2  + ^S5S3  + ” 24s^s2)/ML  1 

2 2 2 f 3 1 

+ (-12s£  + ^^S6S2  + ^S4  ” ^S4S2  " 1^S8  + ®S5S3  + J 

+ (-Sg  + s2)/m^2^ 

= (-181440s,  -h  120960s,  s.  + 96768scs,  + 45360s.2  - 45360s.  s2 
o 62  53  4 42 

- 40320s2s2  + 3780s4)/Mf83 

+ (-103680s,  + 69120s, s,  + 55296s, s,  + 25920s.2  - 25920s, s2 
8 62  53  4 42 

- 23040s2s2  + 2160s2)/M[7] 


+ (-27360s,  + 18000s, s,  + 14592s, s,  + 6840s2  - 6660s. s2 
o 0 2 j 3 4 42 

- 6000s2s2  -b  540s4)/M[6] 

+ (-4320s  + 2736s, s,  + 2304s, s.  + 1080s2  - 972s. s2 
8 62  53  4 42 

-912s2s2  + 72s2)/M[5] 


+ (-438s,  + 256s, s + 208s, s,  + 131s. 

o 6 2 5 3 4 

+ 9s2)/Mf4J 


102s4s2  - 64s2s2 


9 r 


I 


L03 


+ (-28s 8 + 12s6s2  + 8s5s3  + 14s2  - 6s4s2)/lJ31 
+ (-s8  + s2)/M[2] 


1 2 8 4 3 8 4 320  2 7 360 

mVL{'~  mI2]  " m[3]  ' m[4]  ‘ M [ 5 ] ' m[6] 


103680  181440 


M 


[7] 


M 


[8] 


JL.  12  . 256  . 2736  . 18000  . 69120  . 120960  . 

C r 1 ^ r a i re  1 ^ r * i **  mi  ^ r / 


+ m6m2iC(  JM  - J4]  - J5]  T m[6] 


M 


M 


M 


[7] 


2.  8 . 208 
+ mcm,M  ( r,i  + r , 1 +■ 


2304  . 14592  . 55296 


96768 


5 3 


M 


131  m[4] 


M 
131 


[5] 


M 


[6] 


M 


[7] 


M 


[8] 


22  1 14  131  , 1080  . 6840  . 25920  . 45360  . 

“4M  ( m[2]  + 13]  + „[4]  + [5]  + [6]  + [7]  + [8]  > 

n ci  ci  n n n n 


2„3  . 6 102  972  6660 

®4m2?r(  ' .,[3]  ' „f4]  _ m[5]  ' m[6]  " 


25920  45360 


Ml 


M 


M 


,[7] 


+ m3m2 


2 JS,  64  912  6000 

( “ fAl  ■ 


Ml 

23040  40320 


M 


[8] 


M 


[4]  M[5] 


M 


[6] 


M 


[7] 


M 


[8] 


4 4 . 9 . 72  ^ 540  . 2160  ^ 3780  , 

“2  ( Mf4]  M[5]  M[6]  Mf7]  M[8]  * 

pi  n n ci  w 


(C34) 


This  completes  the  present  task.  Expressions  for  the  terms  k^t  ^42’ 
and  k/t/|  appearing  in  the  estimates  of  the  variance  for  the  jump  rate  are 
now  available  in  terms  of  sample  moments,  allowing  an  immediate  estimate 
of  the  errors  from  the  experimental  results  themselves. 


I 
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APPENDIX  D 

COMPUTER  PROGRAMS  AND  TABLES 


' 


PPLOT 

THE  PURPOSE  OP  THE  P POGRAM  PPLOT  13  TO  CALCULATE  ANO  PLOT 
Th€  PROBABILITY  THAT  AH  AO A TOM  STARTING  AT  TN£  ORIGIN  HOPS 
TO  A POSITION  »I»  LATTICE  3PACINGS  TO  T*E  RIGHT,  THIS 
PROBABILITY  IS  ASSESSED  *0*  VARIOUS  VALUES  OP  THE  JUMP 
PaTe  HaTIO  8/A  ANO  THE  SECONO  MOMENT  *2, 

IMPLICIT  PEAL(M) 

DIMENSION  80 A (5) , PuNC (25,5) 

DIMENSION  PI  A (501 ,BI*(50) ,SE3IA(30),S£3IP(50) 

COMMON  P UnC 
TYPE  <?5 
ACCEPT  100, N 
IP  (N'EQ'01  STOP 

SET  THE.  NUM0EP  OP  TE»MS  «N«  EQUAL  TO  ZERO  IN  OPOER  TO 

terminate  execution  op  pplot. 

ACCEPT  ll0,0OA(n 
accept  ita,anA(2) 

ACCEPT  ita,8DA(3T 
TYPE  140 
ACCEPT  113. M? 

TYPE  US 
OU  S3  J ■ I . 3 
DO  53  NSTEP«0,a« 

A«M2/(2.*a.«80A(jn 

BaBO  A (J} • A 

CALL  IC(8.*AtRICA} 

CALL  IC (2 , «8 , RICB) 

CALL  INuE(2.*A,N,RICA,PlAi 
CALL  INU€C2.*0.N.RIC«.»I83 
NP4aN*4 
NT2a2*N 
DO  10  tai,N 
8ESlA(NP4-n  aPTAfU 
8E3I A (I *Np4 ) ap I A f 1 1 
RESIB(NP4-n  aPIB(I) 

86SI8 ( I+NP4] apIS(I) 

8E3IACNP4) «PICA 
•ESIB (NPa) aPICB 
NOOOa2*(N/2l«i 
N£ VN«2» (M/3) 

TSTPTga (?«NEVN) /2 
IST0P2a(2*NCVN1 /J 
ISTRTl a ( j vNODOI /? 

I3T0Pla(l*N0P01 /2 
08LP2a«. 

00  20  PaISTPT2,IST0P2 
08LP8a08LP2aEXP(-2'.a»l  aBESlB  (K*NP4) 
ia€XP(«2,aAl a8ESIA(NSTEP-2aP*Np«) 

CONTINUE 

PUNC (N3TEP. J1 aOBLP? 
continue 

00  60  Ia3.6 

type  130, I.punc (I, 1 ), punc (I, 23 .punc (I, 11 

POPMAT(»  GIVE  the  * OP  TERMS,  ANO  THE  3 VALUES  P OR  3/A*J 

PURMATCI23 

POPMAT (Gi2*6) 

PU»MAT(2H1  ) 


T * • v ■ ' 

•—  — ■ - 


orioonooo  oooor» 


P-i 
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FROM  COPY  FURN1£>H£D  10  DOC  

130  fqrmat(;h  ,I2,3h  ,3CEi2.3,2H  1J 

1«0  ^USMiTC'  GIVE  THE  VALUE  OF  THE  SECOND  MOMENT') 

lia  FORM  a T ( 1 m , G 1 2 . a ) 

caul  °uur 

SO  TO  5 
ENO 

THE  3u8H0uTIN£  "TC"  CALCULATES  A ZEROTH  ORDER  HQOIEISO 
9E33EL  ffUNCTIQN  FOP  a GIVEN  ARGUMENT  *X"  ano  ASSIGNS  THE 
VAuuC  QE  THE  FUNCTION  to  » R I c ». 

subroutine  iccx.ric) 

9ICbA<J3(X) 

IF  (aIC-3, 73} 1,1,2 

1 Z«*.X*7,iilltlE-2 

BICaf(CCl4.5913£-3*Z»3.a078aE«2)  • Z *2  . b59T  32E- 1 ) a Z*  l . 2067  A9E0 ) «Z 
l*3.0«R«»42Ea)  .Z*3.5l5fc23Eaj 
RETURN 

2 Z.3,7^/Ric 

9ic»  E*p (sin /1U»T  fRin  *(( ( Cf  C C C3.’2377E-3*Z-1.6a7fc33E-2) *Z 
U2,b35537c-a;  .Z-2.35T706E-2)  »Z*5'.  18281E-3)  «Z-1 .57E65Eb3)  *Z 
2a2.25319f-3) «Z*l .J285R2E-2)  *Z*1,  R89a23E-n 

RETURN 
ENO 


The  SUBROUTINE  "TNIJE"  CALCULATES  A set  or  MODIFIED  3E3SEL 
functions  of  oroers  i through  »n»  fqr  a given  argument 
"X"  ANO  The  value  OF  the  CORRESPONDING  ZEROTH  OROER  3ESSEL 

function  »zi*.  the  values  of  the  vessel  functions  fop 
0RDEs3  t THROUGH  »n"  are  assigned  to  the  one  DIMENSIONAL 
ARRaV  “RI". 

SUBROUTINE  inue  C* / n , ZI # R I ) 
dimension  PK3H) 

IF  (Njia,:a,i 
l Fn«N*N 
0 1 * X /FN 

IF  (A8S (xj -5.E-45 a, a, 2 

2 *0«l. 

*i»a. 

80*0. 

«1*1  . 

FI  aFN 

3 FI.FU2 
ANbFT / A8S  (X) 

A a AN  * A t * A0 
9b AN«a I *B0 
*0«At 

80»ei 
A 1 a A 
3 I »8 
T0«Q  1 
qi aA/b 

IF  (AB3C(31a501/on-I  .E-8)4,a,3 
a IF (x ) 5 r fe «a 
5 C 1 a*(j l 
8 <aN 

T 0 1 a X / CFNa  X *Q 1 ) 

R l (K)  «<Jl 

FNaFN-2 
A «K«  l 


' . 


^ --jl  •_  ....  •' 
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iF(K)a,a,7 
a ^i»zt 
00  *»  I«t,N 
f:*fi*bki] 
0 »I(I)«fi 
19  9gTU»N 
fNO 


C T.-(C  SUBROUTINE  "PLOT"  PLOTS  ThE  PROBABILITY  In  QUESTION 

C *3  A function  OF  lattice  POSITION, 

c 

SUBROUTINE  plot 
»EAL  LINE 

DIMENSION  LINE (90) , FUNC I 25 . 5) , ALPHA (13) 

COMMON  FUNC 
TYPE  I 7g 

accept  taa, scale 

Type  160 

ACCEPT  120, (AcPHACI] 

TYPE  150 
00  30  1*1,77 
30  LINE(I)«ALPHA(4) 

TYPE  110,  CLINEdl  ,1*1,771 
00  50  1*0,6 
00  35  K*2,7t, 

35  LINE(K)*ALPmA(5} 

00  40  Ja 1 , 3 

IFUNCaFuNCa.J1*T2./SCALC*.5 

*0  LINEaFUNC*2)aALPHA(J) 

5«  TYFE  130,1, (LJNEfX), Hal, 771 
100  F0PH*TflHl, A77J 
1)0  FORMAT ( 1 H ,77*1] 

120  PQPMATCfcAtl 

13a  FORMAT^*  , 12,7741) 

150  FORMAT  (2Hl  ) 

160  P0PMA7I*  SIVE  TH?  »ivc  AL.*man(jhe4IC  values  (E.S,  123,  ) l *] 
170  format C * UM*T  13  the  SCALE?"} 

180  Format (512.6) 

CnO 


on  n r»or»no  nnnnnnnnnnnnnonfi 
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rais  pa®  IS  SKI  ausLin  mcnciM* 

S«pi^““10005  — 


KPLUT 

THIS  PROGRAM  -kWLUT-  IS  0E3IGNE0  TO  CALCULATE  THE  VALUES  Fmr  T“£ 
STanOaSO  DEVIATIONS  OF  ALPWA-T  ANO  9£TA-t,  The  FOLLOWING  ST^a^LS 
"ILL  a€  USED  ThrouGHOU  T ThE  PROGRAM,  ano  "ILL  3E  GIVEN  THE 
ADJACENT  OEFTNTTTUNSj 

ALFT...THE  P9QOUCT  OF  -ALPHA"  ANO  "T". 

PTAT...THE  PPOOUCT  OF  -9ETA-  anq  -T", 

H2  ...TH£  2*Q  P *R£NT  «0h£NT. 

Eh«a...TH£  STANOASO  ngviATTCN  OF  ALPT  OIVIOEO  3Y  aL?T, 

E9P0...THE  STanoaHO  OEVJATI3N  OF  3TAT  OIvIOED  3Y  9TaT, 

I NO  X ... A LABEL  FOR  THE  24  VALUES  OF  «2. 

NOTE,  That  FOR  THIS  PROGRan,  the  oISTance  BEThEEN  aqJaCENT  SITES 

is  assumeo  to  se  unity, 
implicit  SEAL  C< , m , n) 

niMfiNSlON  F'jNC  (25,51  , aP ( 8 J 
COMMON  FUnC 
ITABLE»0 
1J  TYPE  180 

ACCEPT  150, M 
IF  CM.LT.l.)  STOP 
IM«M 

ITA6LE*ITAPLF"1 
Type  140 

type  170,ITAPL6,M 

TvPE  130 
TYPE  130 
Type  leu 
00  50  J«t  ,S 
type  130 

00  50  C«l.,24. 

RJ*J 

! sC  + 1 • 

*•0/8. -2. 

M 0 * 1 0 ,-«  * * 

A IS  EQUIVALENT  TO  ALPHA-TJ  9 rs  EQUIVALENT  to  3ETA.T, 

M3  is  EQUIVALENT  TO  the  MEAN  SQUARE  MISPLACEMENT  OF  THE 
DIFFUSING  AOATOM, 

A»M2/(2.".8«PJ5 

8*»J*a/10. 

I0»2* J-l 

1 X •«£ » J 

this  part  of  the  program  calculates  fsra  & eprp, 

N8«512,»0Y2.*A"B45t2.*P**2*282.*A**2*l?3R6.«A«3*20'LS0,i*a 
l*(l*3-*A]  #«2«  1 S80.  » ( 4«9* A 1 «*3*1<»80.  » ( 4 »3*  4 ) »»4 
N6« 1 28 , *3*2 , *A"S0. »{4. *9* A } «*2*?20. .3* { 4 . *3* A ) • 1 20 . 

1 • (4 , »9* A 3 • • 3 

N4»2.»A"32.»8Mt2.«(4«a*Al ««2 
N2»a,»C4,.3#A3 

Ka«N8-28,»N6**i2,3S . »N4**2»420,*N4»N2«*2-830,»N2*«4 

*b«N6-l5.»Ntt«M2"30, *N2**1 

K4«N4-3.«N?**2 

*2»n2 

VA«  (Fd/M*u,.A6«F2/  CM-t  .3  ♦(m>33.1  *K4«.2  / fM-1 . 3 
1*72. »M*A4»«2«,2/ fM-1 . 3 / (M.2,3 


HBi 
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l 

1 

1 

T 

t 


l *04, l. 5 «K2««4/ (M-t .3 /(M-3.1 / (H-3.1 
1»4.i*6/m-8.*(m*7.3  *K«i*4/(M»l .1 

«K2««2/ (M-t ,3  3 /36,-A««a. 

,3 ♦(*♦33.3 aK4a#3/ (M-i.) 
fM.1.3/ (**2.1 

1*44, aMa(**l  .1  •*£*•*/ (M*l.)  /(M-a.1  / 0-— 3 . ) 
l-4.*Kb/M.2.*C«*7.3  *K4*IC2/  (M.t .1 

(M.i ,) 1 /S7b.-4«*2 

STAaVA..,5/A 

3T8«VB**,5/8 

(1.1*3 »ALOGia(3T*3 
FUNCd.  Ill  »AUJGia(3T4  3 
C 

C ThIS  3TArENg»4T  CREATES  TH£  Table  ENTREES. 

C 

IOC 

34  TYPE  ll0,A,8.Ma,3TA,AT8, IC 
30  CONTINUE 

10VJ  FOPmaTCIh  itrtH  ,*4. 4H  33 
im  4omh4T(im  ,sEia.6,rsi 
120  foPhaT(aa) 

130  FU*H4T(2H  3 
140  *00*4712*1  3 

190  4O»M4T(B10,43 

160  FORMAT  ? I*  ,4X.  »ALPT*,0X,  *8TlT»,s\f , '*a*  . 10X. 'ERRA'.Sa, 
I'ERRS'.dX,  d»»0X«3 

170  FuRMATd*  .*T*0tP  n'.Il»5*«**»  *'.£9.33 

140  FQRM4TU*  , *WMAT  IS  the  number  OF  OBSERVATIONS?*} 


This  SUBROUTINE  PLOTS  ERR*  ANO  ER99  AS  FUNCTIONS  of  M2  on  a 
THREE  TIEREO  LOS-LOS  SCALE.  the  key  8EL0H  INOICATES 
Which  symbols  CORRESPOND  TO  which  values: 

1.. .ERRA  WHEN  8«A/10 

2.. .ERR8  WHEN  BaA/10 

3.. .ERRA  hHEN  is a/9 

4.. .ERR8  WHEN  9aA/9 

CALL  PLOT 
GO  TO  10 
ENO 

SUBROUTINE  plot 

real  line 

dimension  LINE (903 . FUnC  (33.51 . alPmA(101 
COMMON  FUNC 
TYPE  160 

ACCEPT  130. (ALPHACI3 .1*1.61 
TYPE  150 
00  30  Ot, 77 
LINfe(I3«AL?HAC5) 

TYPE  118. CLINE (II .1*1 ,771 
DO  50  r»l,34 
00  35  K«a,76 
HN£(4)«ALPHA(6] 
nu  40  Jal.a 

IFUNCaFuNCd,  J1  •a5-»4  9,5 
LINE (IFUNC*ai a ALPHA (J) 

I M 1 a I • l 
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TrP£  130,  1*1  , CL  I we  f K •)  ,Kat  , T71 
FUSMATClHl , A771 
ruRMATfiH  # 7 7 A 1 ] 

FURMATfjjAp 
FURMATf^H  ,ira,77AlJ 
FURMATCIHH 

FuRMATfj^  # *WR ITE  1234.  AnO  A <3lan**1 
EnO 


no  n nrinoo  or»  oo  r»r»  r»  o r»r»  r»o  r» 
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"PLOT1 

this  program  -"plot-  is  oesigneo  to  calculate  the  values  fob  the 
STanOaRQ  deviations  OP  alpha. r AHO  3ETA.T,  The  EOLLOhIHG  SYMBOLS 
HILL  9E  USED  THROUGHOUT  THE  PROGRAM,  ANQ  WILL  8E  GIVEN  TH£ 
adjacent  OEFINITTONSJ 

EHRA...THE  STANOAHO  0EvI ATION  OP  ALPT  OIVJOEO  8Y  ALPT, 

ERRB...THC  STANOARO  DEVIATION  OP  8TAT  DIVIDED  8Y  3TaT, 

INOX.. .A  LABEL  POR  T«E  24  VALUES  OP  THE  SeCONQ  MOMENT  M2. 

NOTE,  THAT  POP  This  program,  THE  OIJTance  BETWEEN  AOJACENT  sites 
IS  aSSUMEO  to  BE  UNITY. 

IMPLICIT  real  (K.M.NT 
0IH£N3I0N  PUNCf25,P] , aPC8) 
common  ptJNC 
IP  TYPE  180 

ACCEPT  150, M3 
DO  50  J»  t , 2 
TYPE  130 
RJ«J 

DO  S0  1*1,25 

A IS  EQUIVALENT  to  ALPHA-T;  8 IS  £0U r V *L£NT  to  8ETA.T, 

M2  IS  Equivalent  to  the  mean  SQUARE  displacement  op  the 
DIPPUSING  AO ATOM . 

A«M2/C2.*.a*RJ7 

R*RJ*A/10. 

I0»2*J-1 
1 1 *2*J 

this  °art  op  the  program  calculates  erra  4 errs, 

N8»512.»R*2,  *A**45  12. «8». 2*252. *A**2*!20R6,  * A *8*32 1 80  , *8 
l»(4*S>A}.,2Qt  752’.  ■ (4«a*A3  ••3*l8a0,*C4*8*4).*4 
N6*i28,»B.2.»A*b0f.(4,.8*45*.2*720.*R«(4,*8*A)*l20, 

1*(4,.U.A1**3 

N4«2.»A.32.»«*t2.*f4*d-*A)  «*2 
N2*2.*<4.»«*A) 

Ka«N8-aa.  »Nb*N3-35.  «N4*.2*42<?,.N4*N2»«2.A30.  *N2»»4 
Pe*Nd-15. *N4#N2*30. »N2«»3 

K4aN4.3.«N2«»2 

K2«N2 

VA«  (Ka/M*  16.. K6.K2/ (M.t  .)  *(M-*33.)  *54**2/ £N-t.) 
l*72..M*54,X2„ay  (M-l’.j/CM-a.) 
l*24,*M*(M*t,I*K2**4/ fM-t .) /CM-2.5 / (M.3.J 
l-8,*KB/M-a.,£M*7.3  *54*52/ £M.l .) 
l*16,*44/M*ts.*(M*r.J*K2.*2/£M.t,)3y36,.A**2. 

V8*(Ka/M*i6,.K6*5  2y  CM- 1 .1 *(M*33.1 *5»**2/(m.i.3 
1*72, *M *Ka*52**3/ CM- 1.3 /CM*a,I 
l*24,.M*(M^t ,1 *<2**4/ fM-t ,J /(M-3.)/£m.3.3 
l-2.*5d/M-2.*fM*7.3 *54*52/ CM. I ,3 
1*54/m.^m*i  .3  .52**2  / CM-1  ,1  3 /37fc‘..A«*2 
STA*VA..,5/A 
STS*V8*», 5/8 
PUNCH,  I»]  *STa 
PUNCH, II)  *378 
50  CONTINUE 

100  PORMATHM  ,of4H  ,A4,4M  )3 
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113  FuPmATCI*  ,5E13.6,Io1 
124  F(JOm  A T ( A4  J 

133  FORMAT^  ) 

1“3  PQPH4T(2H^  ) 

154  *URMAT  i-Gt2.e) 

lf»y  FOPMATtlH  ,'WHAT  is  THE  SECONO  MOMENT? ' 3 

c 

c The  SuPPOUTjnE  "BtnT”  B(_qt3  EPPA  ANO  £PPB  as  JUNCTIONS  QF  ThE 
c NUM3£N  of  DIFFUSION  INTEPVAUS  "M»,  the  KEY  BELOW  InOICaTES 

c which  s y MgOL s coppespono  to  which  values. 
c 

C 1 ...EHRA  when  3 a a / t 3 
r 2...EHP8  -*HEN  **A/14 

C 3...EHPA  ,H£N  a*A/«j 

C “...Eh »8  when  8*a/s 
C 

CALL  PLOT 
ENO 

SUEWOUT I N£  PLOT 

peal  line 

Dimension  LINE  (93) , FhnC(25,51 , alpha  Cl ai 
COMMON  FUNC 

type  163 

ACCEPT  133.  (ALPHACn  . I « 1 « b 3 
Type  153 
00  30  1*1,77 
33  LINE ( 15  * ALPHA (5) 

Type  1 13. CLINEf n , T»i ,T71 
OU  50  1*1,3“ 

DO  35  **3,7!, 

35  LIN£(KJ*ALPHAC6) 

DO  “0  J*l,“ 

IFUNC»FUNC { I , J 1 *73 , 5 

43  LINE(IFUNC*21*ALPHA(J1 

IM1«I-1 

50  Type  133. IMl. (LINEfKI ,*«1 ,771 
100  . FQPMAT  ( I H t , A771 
110  FiJPH  A T f ^ H ,77*1) 

12-4  F OPM  A T ? b A 1 1 

134  *UPmaT  f j n , 1 12 , 77 A 1 J 
154  F0SMATflHt) 

164  FOBMATIIH  , 'WRITS  1234,  ANO  A PLANK  * I 
ENO 


n * 
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TADl  £01 


M=  0.10E+03 


At  PT 

.476238E-02 
. 633100E-02 
. 846919E-02 
.U293SE-01 
. 150606F-0  L 
> 200336E-01 
. 267S19E-0l 
.2S7143E-01 
.476258E-01 
.635100E-01 
.846919F-01 
. 112938E+00 
. 150606E+00 
. 200836E+00 
. 26731 9E+00 
• 3371 43E+00 
. 476233E+00 
•635100E+00 
. 3469191+00 
. 1 12933E+01 
. 150606E+01 
. 200S36E+0 1 
. 26781 9E+01 
. 3371 43E+0.1 

» 370423E-02 
. 493967E-02 
.638713E-02 
. 37841 0E-02 
. U713SE-01 
.156206E-01 
. 2G3304E— 01 

. 2777TgE_Ai 

. 370423E-0 1 
.493967E-01 
. 6337133 -01 
.3784101-01 
. 117133E+00 
. 136206E+00 
. 20S304E+00 
. 277778E+00 
. 370423E+00 
. 493966E+00 
. 63371 5E+00 
. 37341  OE-f-00 
. 1I713SE+01 
. 136 20 6 £—01 
2033041+01 

■i  — r*  - • 

. - / ■ - i. 


a TAT 

. 47623SE-03 
» 635100E-03 
♦8469191-03 
. 112933E-02 
. 150606E-02 
. 200S36E-02 
.267919E-02 
. 3571 43E-02 
. 476233E -02 
. 635100E-02 
. 846919E-02 
.1129381-01 
. .1 30606c  —01 

• 2.00336E-01 
. 267S19E-01 
♦357143E-01 
• 476233E-01 
.6351001-01 
.846919E-01 
. 11293SE+0Q 
. 130606E+00 
.2 OQ836E+00 
. 267919E+00 
. 3571 43E+00 

• 740845E-03 
. 987933E-03 
. 131743E-02 
. 173632E-02 
.2242761-02 
. 3 12412E-02 
. 416608E-02 
. 333336E-02 
. 740S45E-02 
.9379331-02 
. 13.1 743E-01 
► 173632E-01 
. 234276E-01 
• 3124122-01 
. 416o0SE-01 
.335336E-C1 

• 740845E-0 !. 
.987933E-01 
.1317  43E+00 
. 1 73aS2S+00 
• 234276E+00 
.31 2412E+00 
. * 166081+00 

• 3.':' 735  *£-'•00 


M2 

. 133352E-01 
. 177S29E-01 
. 237137E-01 
.31 6223E-01 
• 421697E-01 
.S62341E-01 
.749894E-01 
. 100000E+00 
. 1 33352E+00 
. 1 77328E+C0 
. 237137E+00 
.316223E+00 
. 421 697E+00 
.362341E+00 
..  749894E+00 
. 100000E+01 
. 133332E+01 
. 17732SE+01 
. 237137E+01 
. 316223E+01 
.4216961+01 
. 562341E+01 
. 749394E+01 
. 100000E+02 

. 1 j^E— 0 1 
. 177S2SE-01 
. 237137E-01 
. 31622SE-01 
.421697E-01 
. 56234 11-01 
.749894E-01 
. 100000E+00 
. 133332E+00 
. 17732SE+00 
. 237137E+00 
. 31 6223E+00 
. 421697E  +00 
. 36234 1E+00 
. 74?594Er00 
. 1000001+01 
. I 33352E+01 
. 17732SE+C 1 
. 237137E+01 
♦ 31622SE+01 
. 421a96E+01 
. 362341 E+01 
. 749394E+01 
. 100000E+02 


ERRA 

. 1 1 4315E+01 
. 102253E+01 
. 922200E+00 
. 839759E+00 
. 7730°7E+00 
.720373E+00 
. 680033E+00 
. 630772E+00 
. 631614E+00 
. 621 936E+00 
. 621 52SE+00 
. 630654E+00 
.630133E+00 
.681 430E+00 

.726787E+00 

. 789409E+00 
» 3737 43E+00 
. 985S63E+00 
.1134031+01 
. 132945E+01 
. 15S72SE+01 
. 19279SE+01 
. 237906E+01 
.297744E+01 


. 14320 

. 13099 

.12109 

.11320 

. 10704 

. 10237 

.93986 

.96736 

.9531 

.93264 

.95993 

. 97767 

. 10072 

. 1050? 

.11113 

.11944 

. 13046 

.14306 

.16431 

. 18965 

.22303 

.26708 


1E+01 
11+01 
6E+01 
E+01 
3E+0 1 
0E+01 
2E+Q0 
2E+00 
5E+00 
5E+00 
3E+00 
3E+00 
61+01 
11+01 
2E+01 
2E+01 
3 E+01 
6E+0 1 
1E+01 


II  E+01 
ISE+01 
I6E+01 
.3253301+01 
. 40230 l E+01 


ERRS  . INDX 

. 359410E+01  1 
. 320939E+01  2 
. 28S832E+01  3 
.262313E+Q1  4 
. 240699E+0 1 3 
.2233861+01  6 
.2Q9845E+01  7 
. 199626E+01  3 
♦ 192359E+01  9 
. 137774E+01  10 
.133711 E+01  11 
. 1361 43E+0 1 12 
.1391911+01  13 
.1951581+01  14 
. 204561E+01  13 
.2181841+01  16 
. 2371 48E+01  17 
.263006E+01  13 
. 2Q7975E+0 1 1* 

. 34461 8E+01  20 
.4070891+01  21 
.4904521+01  22 
. 601620E+01  23 
• 749836E+01  24 

. 238271 E+01  L 
. 25741 3E+0.1  2 
.2316S6E+01  3 
. 210356E+01  4 
. 1 92963E+0 1 3 
.17899SE+01  6 
. 16798SE+01  7 
. 1595S8E+01  3 
. 133471E+01  9 
. 149393E+01  10 
. 147193E-L0I  11 
.1468021+01  12 
.1482341+01  13 
. 131709E+01  14 
. 157477E+01  15 
. 166039E+01  16 
.173 101 E+01  17 
. 194630E+01  13 
. 21 7036E+01  19 
. 247092E+0 1 20 
.2872841+01  21 
.3409231+01  22 
.41 2451 E+01  23 
. 50 77° 6 E+01  24 


y ...  L 
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AL  r T 

3 r 4 : 

7 2 

ERR  A 

E:VR3 

[N 

47623SE-02 

. 476253E-03 

. 1 3 3 3 3 2 E - 0 1. 

* j 6 ^ !.  juClt  w? C! 

. i 1 3537E+0 1 

i 

433  L 

. 0 35 1 .)0E  -03 

. L 77329E-01 

.32297  ”=>00 

, 101379E+01 

■» 

346 3 t 3E -02 

.3463l3E-0'3 

• 237 1 3 "'E  -0  i 

. 2*5.'  201  £+00 

. 9121.L3E+00 

T 

L l 2938 E- Cl 

► k U w ■ ' *- 

*7  ! 1 *> O Q ET  1 

* j iO*.  w <z.  — i. 

. 263077P+00 

. 32S107E+00 

4 

ISO iOa£-0 1 

. 150606E-02 

, 4216775-01 

. 2 4 393 7 E +00 

. '"3?59cE-00 

5 

2O0836E-O 1 

• 2003 76 £-02 

. 3o234  i 

. 22720  IE -00 

. •’046  7 4E+00 

o 

2.»78l9E-0t 

.2  ••  731  5E -02 

.749994E-01 

* v_  .*.  4 ;> 0 4 f 00 

. so  1.  j 6 0 £ + 00 

“2* 

33  ”1433  -.)1 

. 33  21 4.*:  £-02 

. ICOOOOE+O 0 

. 20S033E+00 

. 0 2° 1 43E+00 

;6 

476258E-0  1 

. 476253E-02 

. 13 33 3 2 £+00 
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TABLES  OF  POWER  SUMS  AND  k-STATISTICS  IN  TERMS  OF 
AUGMENTED  SYMMETRIC  FUNCTIONS 

In  Algorithm  Bl  of  Appendix  3,  we  outlined  a procedure  for 
finding  the  unbiased  estimator  of  a given  population  value.  We  found  that 
this  estimate  could  often  be  simply  expressed  in  terms  of  k-statistics. 

In  order  to  generate  a given  k-statistic  from  a data  set  however,  it  is 
best  expressed  in  terms  of  sample  moments  about  the  mean.  By  expressing 
a k-statistic  in  terms  of  augmented  symmetric  functions  and  then  finding 
these  augmented  symmetric  functions  in  terms  of  power  sums,  we  can  easily 
write  the  k-statistic  in  terms  of  sample  moments  about  the  mean. 

Because  the  expressions  for  k-statistics  in  terms  of  augmented 
symmetric  functions  and  augmented  symmetric  functions  in  terms  of  power 
sums  are  often  long  and  complex,  these  quantities  are  best  expressed  in 
terms  of  one  another  via  tables  such  as  Tables  32  through  B5.  For  the 
reader's  convenience,  we  have  collected  and  reproduced  many  of  the 
relevant  tables  from  the  literature  in  this  appendix.  These  tables  are 
reproduced  from  four  sources.  Table  El  is  taken  from  Kendall  and  Stuart 
[14],  Table  E2  is  from  David  and  Kendall  [13],  Table  E3  is  from  Wishart 
[11],  and  Table  E4  is  from  Abdel  Aty  [12]. 
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To  express  the  [ ] functions  in  terms  of  ( ),  read  downwards  up  to  and  including  the  mam 
diagonal,  e.g.  [41*]  = 2 (6)-2  (5)  (1)  — (4)  (2)  + (4)  (1)*.  To  express  the  ( ) functions  in  terms 
of  [ J,  read  across  up  to  and  including  the  main  diagonal,  e.g.  (4)  (1)*  = [6]  + 2 [51]  + [**]  + U1’!- 
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TABLE  E4  (continued) 
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-7560 

6720 

-2400 
9960 
- 14280 

-7920 

34560 

-6720 

6720 

-XX 

“ II824O 
I8I44O 

*3440 

“5040 

-20l60 

I848O 

“ 1680 

- 1260 

7980 

8820 

-88200 

-5040 

210 

- 1260 

-630 

12600 

-X 

1120 

- 1120 

-4200 

20160 

-1680 

1680 

- 1680 

2800 

6720 

-42000 

560 

-2240 

560 

140 

210 

-1890 

- 1260 

21000 

• 

210 

-2100 

140 

-140 

-360 

2800 

140 

- 1400 

-210 

-420 

4»o 

1980 

-8640 

1680 

- 1680 

630 

1050 

-\x 

-6lOO 

30240 

-2520 

4200 

-630 

-630 

5670 
- 1260 

“3*500 

420 

-2940 

210 

-630 

9450 

420 

-35 

-280 

280 

1890 

-9240 

560 

-560 

70 

3*5 

“595 

-2520 

14700 

560 

“35 

3*5 

3*5 

-4200 

. 

“35 

35 

280 

- 1400 

. 

-35 

350 

-210 

1050 

“35 

35 

210 

- 1050 

“35 

“35 

*75 

• 

• 

42 

-16$ 

-84 

-252 

1728 

“336 

336 

-105 

252 

630 

“5544 

336 

-672 

84 

63 

“23* 

-378 

5040 
- 1260 

336 

“21 

63 

-h 

42 

-42 

1680 

“S6 

56 

-a, 

63 

84 

-2100 

“56 

-21 

420 

140 

42 

-420 

• 

-21 

3*5 

“35 

, 

. 

42 

, 

-21 

“7 

*4 

“ *4 
21 

*4 

“35 

-n 

-168 

420 

56 

-28 

“7 

21 

21 

-210 

-28 

“7 

7 

28 

-140 

. 

“7 

• 

70 

“7 

35 

“7 

1 

2 

-2 

-6 

24 

-8 

8 

-2 

“J 

5 

12 

-60 

-8 

X 

-3 

“3 

30 

« 

— 1 

“4 

20 

3 

1 

X 

- 10 

3 

4 

X 

“5 

*5 

10 

10 

5 

X 

X 

1 

- 1 

X 

8 

8 

• 

1 

1 

2 

3 

»4 

J1 

• 

8 

1 

9 

6 

1 

• 

36 

9 

9 

108 

36 

36 

9 

1 27 

120 

45 

45 

360 

120 

120 

45 

OO 

„8S 

1320 

*65 

495 

7920 

7920 

3960 

792 

475 

2970 

hi 


I Weight  1 2 (vi)  , 


Jf® Is  BEST  quality  practicable 

COPY  PURNISHK.ri  ejjq 


TABLE  E4  (continued) 


!i  I 


ail*l/n'“ 
a*i»  /n<“ 
a*i*  /«<•• 
a*i*  /»'•> 
a‘i 

a*)  /nw 
3i  •]/*"•' 
3a  i 

3a*i)/n‘«' 

3,ai4]/n»t‘ 
3*a‘ i *)/»<•» 
3,a,J/*»u' 

3,i,1/"'* 

J’aij/ir"’ 


4al•|/»«,•, 

4a*«*|/«‘*> 
4*41/"<4' 
43»M/"<’' 
43*  iW*' 
43»,«  /»<*> 
43,l,]/n<*> 


-♦53000 

-SSiS 


- 162880  I 
18)4400  , 


43*  iW*' 
43»,«  /»<*> 
43,l,l/’><*> 

♦»i  •(/»<•> 

4,i,i/m« 

4*Jl]/«“» 

4,)l«l,> 

Sl’J/it'*’ 

JH  •)/*'” 

Si'il/n'*’ 

S3I*]/"'*’ 

sa2i,)/«“> 

34l,)l»l‘> 

J41l)/»,*> 


(6a  1 •}/«'•> 

63ail/ll<« 

«3*l/4'*' 

641 

642]/n^•, 

fell]/*'*' 

71  *]/■<" 

7a*i]/»'4' 

73l,V"<*’ 

73*]/n,,> 

74>V«'*' 

75]/"(*> 


oin/iio* 

flu)/"11' 

931/*'*1 

to,  aVu"1 

II.  I ]/■«•* 

I2j/« 


- 2520OO  i 

302400 

-37800 


362880 

- 2177280 
4989600 

- 5443200 
2815000 
- 589680 

22680 

604800 

-2419200 

3326400 

- 1814400 

302400 

252000 

-554400 

340200 
-37800 
(6800 
- I68OO 


- 100400 

176400 

-75600 

-4*00 

4200 

6450 


3628800  ! 

- 1 09 58400  | 

39916800  ( 
-34927200  i 
12474000 

- 1247400 

6652800 

- 23284800 


-8316000 

415800 

33264OO 

-5544000 

1663200 


332640 
-831600 
498960 
-41580 
221760 
- 166320 


-,?3C 

-41580 


- 399  16800 
239500800 

- 538876800 
5S8835200 
-261954000 
44906400 

- 1 247400 
-79833600 

31933440C 

-41912640c 

19958400c 

-24948000 

-46569600 

99792000 

- 49896000 

3326400 

-7392000 

4435200 

-92400 

19958400 

-69854400 

74844000 

- 24948000 

1247400 

19958400 

-33264000 

997Q200 

3326400 

-831600 

- 2079000 

24948 00 
-31*850 
— 415800 

11550 

- 3991680 
« 1975040 

- 9979200 
1995840 

- 3326400 
3qqi6So 


I 

1 

I 

1 

1 

! 


THIS  PAQE  IS  BEST  QUALITY  PRACTICABLE 

FROM  COPY  FURNISHED  TO  DDC 
APPENDIX  F 


FORMULAE  FOR  POWERS  AND  PRODUCTS  OF  k-STATISTICS , ORDERS  2-8. 11 


2nd  order 

if  33  i3/^  4*  *u 
3rd  order 
*2*1  * kjn  + k22, 

*1  = *3/ft2  ■+■  3k2Jn  4-  km, 

*11*1  = 2kn/n  + knv 
4 th  order 


*3  *1  — kjn  + i3l, 
if  = kjn  4-  (n  4-  1)  AM/(n  - 1), 
i3  k\  — kjn2  4~  2k3Jn  4~  k^jn  4-  k^n, 
if  = kjn3  + 4 kn/n2  + 3 k2Jn2  + 6i2ll/n  + iuu. 

i3i*i  = k3Jn  + k2Jn  + £jh> 

*2*11  = 2 k3Jn  — 2k2Jn(2)  + fcni, 

inif  = 2k3Jn2  + 2k2Jn2  + §k23Jn  + kmi> 

*111  *1  = 3i2  n/n  + *uu> 

*11  — ‘Zk&ln1®  + 4k21Jn  + klUv 


5<A  order 

k2k3  = kjn~^k22, 

k3k2  = kjn  4-  (»  4-  5)  i32/(w  - 1 )> 

i3if  53  kjn2  4*  2knjn  + k32jn  + k332, 

ifix  = kjn2  + ktjn  + 2{n+  l)  k3J  ni2)  + (n  + l)knJ(n—  1), 

*2*1  “ k$jn2  4-  3 k2Jn2  •+■  ^k32jn2  -¥  3 k333]n  4*  3k22Jn  + k2m, 
if  =*  ij/n*  4-  5 iu/n3  + \0k3Jn3  4-  I0k3ujn2  + 1 oknJn2  + lOk^Jn  + &mn- 

*21  *1  “ *4i/w  4-  iM/  Ji  + i3u , 

*3*11  = 2i41/n  — Bijj/n1*' + i311, 

*22*1  = 2iM/ra + ij2i, 

*21  *2  — ^nin  4-  (ra  — 3)  i32/W(2)  + (n  + 1 ) k22J(n  — 1 ), 
i'21  if  “ k2Jn2  4-  3k3Jn2  4-  2i311/n  4-  3k22Jn  + k3m, 

*2*11*1  =*  2*«/nS  + 2(»  - 3)  i3j/(7i7i‘4')  4-  3i3ll/n  4-  2(rc  - 2)  ijsJn®  4-  ijm, 

*11*1  = 2*«/n3  + Qk3Jn3  4-  1k3lJn2  4-  I2k222ln2  4-  9*2m/w  + *mii> 

*111*1  = *sn/n  + 2*22l/71  + *2111> 

*21  *11  = 2*32/n(1)  4-  2 k3XJn  4-  2(71  — 2)  k32Jn>'1)  + *2m> 

*1*111  = 3*3ii/W  — fiknJnif>  4-  ilm, 

*111*1  = 3i311/n*  4-  3knJn2  4-  ~k211Jn  4-  *nm. 

*11  *1  * 4*ai/ (»»*)  + 4*in>*  + 2(5n  - 4)  *„i/(n«(t>)  + 8 i,m/»  + *iuu» 
*1111*1  = 4*im/n  + *1U11> 

*111*11  “ 6*j»i,2»(n  4-  6*nn/n  "P  *11111- 


I 


oth  order 

|ks*'l  = k<yV,  — k Sl, 

*4*.  = kvn  — (n-r~)kt2/(n—  l)  -4  - l), 

£ij  » £,/A  — 9*42,'(a  — 1)  - (71  — 3)  ^33/  (71  — 1)  - 0A*,.*,  (71  - I)1*’, 

*v*i  = *,,  7i-  4-  2£5l/n  — £4J,  7i  4-  £411, 

£,*, £t  = £„ 71-  * 4Jl( n-(n-5)  £w/ n'a  - ( a 4-  5)  k^i  a'21  - :a  - 5)  -fcja(  (a  - 1), 

<3  » £„  a2  - 3(n  _ 3)  Jci2l  n-2)  - 4a(a  - 2)  &M(  ( A'21)2  +■  (n  -r  i)  (a  ~ 3)  k^in  - l)2. 

*3^1  * *3/  7l3  “4  3<,i/  A2  4-  3*42)  71“  4-  *33/A"  T t 3*333, A — *3133, 

*J  *?  = £„  a3  -r  2*si/a2 -(3n  4- 1)  £42/(aa‘2>)  - 2(n  - 1 ) £33/(71 «‘2>)  - £Wj/n  r4(jirl) kJtvnrJ 
4-(71t-  1)*322/A(2,-t(A4-  l)£sail/(7i-  1), 

*t*$  = £„  »*  t 4kfljn3  - 7£44,a3  j.  ik3Vn3  4-  3£111( 7i2  - 1 BAjjj.  a2  - 34^  n*  - 4£1UI, ti 

+■  0*221.1/  a ~ £21111- 

£^  =■  £,( n4 4* 6i41y'n4  4-  15£42/a4 ~ IOAjj/tv* -4  15£4U, n3  ^60£j21/ti3  -4  \5kKiln3 
— ■20k^livni  -4  -to£„in,  ti  2 — l ok.inlln  — £lulll- 

''u^l  = £51/ 72  ■*.  icitj  ti  — -cui, 

£4  £33  * -^33;  71  — 3£42/7l‘°  — 0£33,  A'J  — £333, 

All  A3  = £«/  71  T £33;  71  4"  £331, 

£33  £3  = ^31/ 72  — “£33;  A'a  — £33/71  — (71  — 5)  £7323/ ( 71  1;, 

£3  £33  = £,3.  71  4.(73  - 4)  £44/7ll2’  - jAjj, A'27  4 (71  -5)  £333/(71  - 1)  - Dn£525/ 7l'3’, 

£3. £3  * 2*,*.  71  - -£33;  A^27  - 171  - 3)  £..233  (71  - I), 


£33 £?  = £33,  71*  ~ 2£42/ 71-  — £33/71-  — 2£4ui  71  — 3 £333/  71  4-  £3333, 

£3 £31  £3  = -£33,  A2  4-  2(  71  — 4)  £ 43j  ( AA'27)  — o£*33y  ( 71 A'*7)  — 3£411;  A — 2(  A 4)  £333/  A'27  — <3333, 

A'iaA'i  = -£«»/ 7,2  ~ -£33/  a2  +•  4 £333;  ti  — £333,  ti  — £32  u , 

£^£,£3  = £j3/7»2  4. 2(n  - 2)  £42;'(7171'2>)  -4  (71  - 3)  k^Knn^)  4-  £4U/a  -4  <3a  - l)  k^njrk-) 

* (A  4-  1 ) £322/ A'27  — (A  -r  1 ) £3033/  (A  — 4), 

£* £31  =*  -£5l/7i2  - 3£42/;(nA<27)  - 2(ti2  - 271 4-  3)  £33,  (a'27)2  4- £4u/ti  -4  4(n  - 1)  £,23,  a'27 

- 471(71  4- 1 ) £35*/ ( a(27)2  -4(ti—  1)  k«3l3j(n—  1), 

£'21  £\  * £51)  71 3 4"  4<C4t(  71 3 — 3£334  71 3 t-  3£41i;  71“  4*  1 3 £ 323 / 71 2 — 3-'' 222/  71 ""  — 3 -C33  33,  *1 

0 £323  37  71  £2U11> 

£,£31  £2  « 2£,i;»J  - 4(7i  - 2)  £*„' (»*»'*)  - 2(71  - 3)  £33/(7i2?i(a)  -4  5*333/ »*  “ 41 3>l  ~ 5)  kmt(nnS*>) 

4-  2(71  — 2)  £335/  (7171' °)  4-  4£3333(  TV  -4  ( O'TV  — i ) £3.333)  7l'a  -4  £33333, 

£33  £t  =-  2£,i/tv4  - S£42)  ti1  4. 6£M;  tv2  - 9£411,  tv3  - 44£,21/7i3  -4  1 2.1-322/  a3  4- 1 0£,1U)  a2 
4-  39£3233)'  A2  -4  1 4£23333)'  A — £333333, 

£311  £'l  = £vu/22+-2£.4d'7VJ-£ji33, 

^31  £'ll  ™ -£12/  7ll2)  -4  2£ul|  TV  4.  2(  A — 4)  £323/  A'-’  -4  £3333, 

£221  £1  = -£321/ 72  — £2221 72  4-  £3233, 

£32  4.33  = “£3 37  A^a  — 4 <323)  A *^£333)  A'a  £3333, 

£3  £333  ~ 3£433,  A — 1 ^£321/  A'"7  — l -^3227  A'3*  ^3 111 ’ 

£'211  £3  ~ £411/ 72  4“  2(A  — 3)  £321/ A'J  — -£322/  A'”7  — ( A 1 ) £3211/  (A  — 1), 

£?1  * £«27  A'8  T £j3)'A'a  - £411/  a - 2(  A - 3)  £,21/  A'21  - (A2  - A - 4)  £„„  A'3’ 

+*(A4-  l)  £2211/  (A—  l), 

£’211  £’l  m £'UU  A2  — 0£j2l/  A2  4-  2 £322)  A-  — 2 £3333,  A — 3 £3233)  A — £33333, 


£23 iu £3  = 2£,2)  ( AA'a)  - 2£j2/(aa'2))  - 2£U1/A2  - 2( 3a  - 8)  ^ ( nA'a)  - 2(a  - 2)  £32*)  (aa'21) 

~ ^*3112,  A — ( 5a  — 1 ) £3233,  A1®  — £'2iiU> 


^2^111^1  — 3*n  , / 71 2 * 6(71  — 3)  *^a,  i ( 7171'-*)  0*o.,<.,  t 
£s£fi  = 4*t2,  (»»®)  - ^33,  (re'a>)2  - 4£111;  re-  - S( 
+■  ■AAjiu/re  7"  — -)  Aojn;  re'-’  — A,llu, 


3)  A311;  ( rere'21)  - 6*222,  (rere'2*)  4 4*3ln;  re  * 3(  re  - 3)  iMU,  re'a)  4 i21111, 
3*  (»*)*  ~ -V*4ii/  re--S(re-3)  fcJ21(  (rere'2’)  - 2(re2  - re  -r  4)  £,M,  (re,2))2 

, _ OS  L 4(2)  _ I. 

* — ) Axiom/  10  'Von  i o , 


*hl£i  — 
Aix  A*  = 

*2111  *i  = 
*211  An  = 
*21*111  = 
*2*1111  = 
*1111*1  = 
*111*11*1  = 


1 4*12/  (rere'2’)  - 4^  (re'2’)2  - 4&,11(  re2  - 5(re  - 3)  *321,  (rere'2’)  - 2(re2  - re  - 4)  *222,  (re'21)2 
4-  4*3  nxjn-r  4(re  — 2)  fc2J11;  re'-’  — A,llu, 

3*4n,  re3  4 13&3J4(  re3  -f  6/Cjjj,  re3  -r  10*3m,  re2  — 2 1 *3344,  re-  -r  12*34414,  re  4 £444444, 

' iiu,  (reV21)  -r  4^33,  (re2re'2’)  - 4*4U/re3  - 8(4re  - 3 )k3Ul(n2n'2))  -Hon-  4)  *222/  ( re^re'2’) 
-r  I2A3444 ; ft3  7-  2(  1 1 re  — 16)  *2211/ (rere'-’)  — 13*31444,  re  4-  '■‘414444. 

*3111 Z72  4-  3fc2211)  n 4 £34444, 

4A334,  re'-*  4 2*333,  re'2*  4 2^3444,  re  4-  2(2re  — 3)  *22n/  re'2’  4-  *21144, 

6*321/  re'2*  — 6A222/re'3*  4-  3*3lu,  re  -r  3(re  — 3)  *2344/ re'-’  4-  kinil, 

4^3444/  re  — 1 2&22u/'re(a*  4-  A21444, 


8*321/  re'2’  — 6A222/re'3’  4-  3A3444,  re  -r  3(re  — 3)  £3344,  re'-’  4-  £34444, 
■*£3111/ 71  — ^ -*wii /re'2*  4-  £34444, 

^£3111/ n2  4-  l^ijjji/n2  4-  9£3U44,  re  4-  kullu, 

12*3,4  / (rere'2>)  4-  6*333 /(rere'2’)  - 6*3444,  re2  4-  8(4n  - 3)  *,344 /(rere'2’) 

1 1 Z»  /■*>  -1—  Z* 


“3111/ 71  + ^-^wn/71  + 9k»iiu>  n + ki 
12*334,  (rere'2>)  4-  6*2M,  (rere'2’)  - 6*,U1, re-  - o^-tre  - 0/  /C3344,  ire 
4-  1 1 *21111/ 71  4"  £411111, 

4*33,  (re'2*)2  4-  24*334/ (rere'2*)  ~ 3(71  ~ -)  *««/’(  t^2’)2  4 8*3444/ 712 
4-  6(5re  — 4)  *3344;  (rere'21)  — 12*34444/  re  — *144414, 

5*siin/n  f £111111, 

1 -£3211/  ^2)  4-  6£2iiii,  re  4-  *111111, 

6£233/re'3’ 4-  ISijjn, n<2> 4-  9£»llll/n 4" *111111- 


£'1111  £11  ~ 


TiA  order 

£,£4  = 
£*£j  = ■ 

£4  £3  = ■ 
*»*!  = 
k^Jc^Jc^  —— 

£3**1  = . 


£.£?  = 
£*£*£!  = 


kfjn  4-  *ji, 

*7/re+(re4-9)*S2/(re-l)4-20*43/(re-  1), 

*7, 're  4-  12isa/' (re  - 1)  4 (re  4-  29)  *33/(71  -1)4  3*re*221/(n  - 1)<», 

*7/n*  4-  2*34/ re  4-  *Sj/n  4-  *511, 

*7/re2  4-  *,Jre  4-  (re  4-  7)  *Sj/re'2*  4-  (re  4- 19)  *«,  re'2*  4-  (re  4-  7)  *m/(re  -1)4-  8*jjU(re-  1), 
iT/re2  4*„/re  4 9*53/ re'2*  (2n  4-  25)  i^/re'2’  4 9 i„i/(n  - 1)  4 (re  4 8)  *334,  (re  - 1) 

4 18*3jj/  (re-  1)'2*4  5re*,23i/(re-  l)'2*, 

*7/re2  4 2(re  4 7)  kitjni2)  4 (re2  4 22re  - 35)  *w,  {(re  - 1)  re'2’} 

4 (re 4 5)  (re 4 7)  *333/(71-  l)2, 

*7/re3  4 3*ji/re2  4 3*33/ re2  4 *43,  re2  4 3*5u,  re  — 3*431/71 4 *im , 

*-/re3  4 2fc,4/  re2  4 2(re  4 2)  *53/  (Tire'2*)  - 3(re  4 5)  *43/ (rere'2’)  4 *511/71 
4 2(»  4 5)  *431  /re'2’  4 2(re  4 5)  *334, 're'2*  4 (re  - 5)  kznjn™  4 (re  4 5)  fcjJU/  (re  - 1 ), 

*-/ re3  4*34/ re2  4 3(ti4  3)  *52 /'(rere'2’)  -(3re24  lire  — 25)  *43/  (re'2’)2  4 3(re  4 3)  kitl,n{2) 

4 4(re  - 2)  *334/{(re  - 1)  re'2*}  4 3(re  4 1)  (»4  3)  *333/;{(»- 1)  *2'} 

+ (t»  4 1)  (re  4 3)  *333i/(re—  l)a, 

*7,re4  4 4*34/ re3  4 6*53,'re3  4 5 *43/ re3  4 6fcsli/  re2  4 12*121/rea  4 4*334,  re2  ^ 3*333/' re2 

~ **£411 lln  + 3*3211/  71  ~ *31111’ 

k.j  n*  — 3*,i/re3  4 ( ore  — 1)  *52/(71 2re'2*)  4 (7re4  5)  *13/!(re2re'2’)  — 3kiujn2 
4 3(3re  4 1)  *434/ (rere'2’)  - 6(n  4 1)  *334,  (7»re‘»)  4 7(re  4 1)  *,»(rere'2’)  4 *uu/n 
-r  6(re  4 1)  iaju/re'2*  4 3(re  - 1)  *,»,/re'2>  4 (re  - 1)  *Mm,  (ti  - 1), 

*7/re3  4 5*,i/  re*  4 1 1*52/' re*  -1- 15*43/71*  - 10*5n/n3  4 35*42i/rea  4 20*Ml/'re3  - 25*333/ re3 

4 10*4411/ re2  — 40*3211,  71,2  7"  15*2331,  re2  7"  '5*31111/  71  7"  1 0*22111/  71  + *311111, 

fc7/re*  4 7*,i,'re3  — 21*52,' re3  ~ 3 5*  43/  re3  4 21*S11/  re*  - 105*m,  re*  4 70*334, 714  “ 10o*32,,  re* 

4 35*4141,  re3  -f  2 1 0*3244, 732  7“  ^ 05*2221/  re3  4 35*3im,  re-  — 1 05*33444, 77 ' 

4 21*344444/71  4 *1444141- 


^ * 4 
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S th  order 
k2kx  — + 

ktJc,  = ifcg/»  + (n+ll)fc„/(n-  1)  30*53/(71  — l)  + 20fc„/(n-  1), 

ksk3  = kjn+  15I82((ti-  l)  + (n  + 44)fcM/'(n- l)  + 30fcM/(n-  l)  + 80wi122/'(n-  1)® 

+ 90nkiaii(n-  1)®, 

k\  = fcg/n+  16 Jfcgi/(n-  1)  + 48i-M, (n - l)  + (n  + 33) kuj(n-  1 ) + 72nfc42g/(n - 1)® 

4 144niM2/(n-  1)®  4 24n(/i  4 l)k2222l(n-  1)®, 
l'» 11  = l’s/»2  + 2A7i/»  4 i-,2/n  4 l-gu, 

i-s  i*2  ij  = ifcg/7i2  -f  t7l//i  4-  (n  4 9)  Jfcg,/7i®  + (n  4 29)  £53/71®  + 20£44/7i®  + (n  4 9)  kS21j(n  - 1) 

+ 20kmJ(n-  1), 

Ms*!  = fc8/*2  + fc7i/n  + 1 2 tjn*  4 (n  4 41)  *„/»»  4 (»  + 29)  fc44/n®  4 1 2JfcMl/(n  - 1 ) 

+ (n  + 29)  ki31/(n  — 1 ) + 36kitt/(n  — I )®  4 72£M2/(n - l)®4  3&nk3221/(n  - 1)®, 

£4/fc|  = ka/n2  + 2(n  + 9)kt2/nti>+8(5n—7)kS3/{(n- l)n(2)}  + (n2  + 2&n-39)kuJ{(n-l)n(2>} 

+ (n-r  7)  (ti  4 9)  k4tt/(n  — 1)*  4 12(n  + 9)i-M2/(n-  l)2, 

£2£a  = fc8/n2  + (n  + 20)  fcg2/7i®  4 2(n2  4 22n  - 32)  fcK/{(n  - 1)  ti®}  + 9(3n  - 5)  kul  {(n  - 1)ti®} 

+ 9(ti2  4 9n  — 20)  ki22j{(n  — l)2(»-2)} 

-i-(n3+  17n2+  104n-  320)  £M2/{(n  — l)2(w-  2)}  + 6n(n  + 5)  k**,J{(n—  l)2(n-  2)}, 
k3k\  — ksjn3  4 3£717i2  4 3£g2y7i2  4 k33l  ri~  + 3knijn  4 4 

/fc4£2£f  = kjn3  + 2k11ln2  + 2(n  + 3)kttl(nni2))  + 2(n  4 13 ) ki3j  (nn(2))  + (n~  19)  kuj(nnm) 

4 knlln  4 2 (n  4 7)  k521lnw  4 2 (n  4 19)  kui  ,7i®  4 (to  + 7)  ki22jnm  4 6£m2/?i® 

+ (w  4 < ) &42ii/  (w  — 1)4*  &k331J  (n—  1 ) , 

£§A1  = fcg/n3  + 2jk71/n2  + (n+8)iM/(nn‘2,)  + 2(n+  17)  kM;(nnn))  4 [2n  4 25)  kuj(nn{2)) 

+ jfcgu/n+  18£521/'ti®  4 2(2n-i-25)  ki3ljn(2)  4 9A422/(n—  I)®  4 (n2+8n  4 20)k332/n(3) 

+ 9kA2u/(n-  l)  + (n  + 8)kmil(n-l)-r36k32ilj(n- 1)®46 k22MI{n-  1)® 

4-  8nk222nl(n  — 1)®, 

k3kikx  = ifcg/n34-fc71/n24-2(n4-7)iM/(nn(2>)4-4(n24- 14»  — 21)  £S3/(ti®)2 

+ (n*  4-2271  — 35)kul(nw)2  + 2(n  4-  7)  ks21jn(2)  + (n24  22n-  35)  ki3lj{(n-  1)ti®} 

4-  (ti  4-  5)  (n  4-  7)  ki22i{{n  - 1 ) n<2>}  4-  2(n  4-  5)  (ti  4-  7)  kM2j{(n  - 1 ) re®} 

4-  (t»  4-  5)  (ti  4-  7)  k3221j(n  — l)2, 

A:2  = jfcg/n3  4-  4(n  4-  5)  kt2l(nn{2))  4-  32(/i  - 2)  iM/(7i(2))2 

4-  (3ti3  4-  23ti2  — 63n  4-  45)  kul{(n—  1)  (ti(2))2}  4-  6(n  + 3)(n~  5)  ki2ij{(n—  1)  ti<2>} 

4-  16(n-  2)  (714-5)  ArM8/{(7i-  l)27i(2)}-t-  (714- 1)  (n4-  3)  (n  + 5)  k2222l(n  - l)3, 
fc4ii  = fc8/n4  4-  4 i71/re3  4-  6kt2jn3  4-  4 &53/n3  4-  ktJn3  4-  6knljn2  4-  12£ui/na  + 4/fc431/7»* 

4-  3 kU2/n2  4-  4 kiluln  4-  6kMn/n  4- 

^ iiiii > 

k3k2kf  = i9/n4  4-  3knjn3  4 2(2n  4- 1)  kt2jln2n{2))  4 (5ti  + 19)  kMl (n2ni2))  4 3(n  4 5)  fc14/  (»iV) 

4 3fcgll/n2  4-  6(n  4-  2)  fcM1/(7i7i®)  4-  9(ti  4-  5)  4m/ (7i7i<2’)  3(n  4 5)  fc422/( nn<2)) 

+ 4(ti  4 5)  fcm/  (nn<2))  4 isui/7i  4 3(ti  4 5)  fc12U/'n(2)  4 3(ti  4 5)  fcM11/ ti® 

4 3(71 4 5)  kmilnlt>  4 (7»4  5)  k32lnl(n  - 1), 

= ktln*+  2knln3  + 4(n  + 2) kt2l (nW23)  4 2(3ti*4  10n—  17)  {w( 77(2>)2} 

4 (37i2  4 14ti  - 25)  l-44/{7i(7i(2’)2}  4 ktujn2  4 6(«  4 3)  ki21j(nn{2)) 

4 2(3ti2  4 14ti  - 25)  Jfcttl/(  t»<2>)*  4 6(ti  4 3)  l-122,;{(n  - 1 ) ti®} 

4 2(37i 2 4 14ti  4 5)  fcM2/(7i®)2  4 3(7i  4 3)  ki2njn™  4 4(n  - 2)  kmil[(n-  l ) ti®} 

46(714  1 ) (ti  4 3)  k322i[{(n  — 1)  ti®}  4 (ti4  l)(n-r3)k2222l{(n-  l)n®} 

4 (»  4 1)  (ti  4 3)  k222llj(n  — 1)*, 
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*3*1  = *8/'«s  + 5A71/tt4  + 10fc82/n4  + 1 1^53/' n4  + oA^/n4  + 10t811>3  + 30A521,n3  + 25  A431/n3 

+ 15A122/n3  + 1 9AM2/n3  + 10Asul/ra*  + 30kiilljn1+  l0kMUj  n2  4- 1 5A3221/n2  -t-  5A4llll/n 
+ 10iS4111/n4-  A3UU1, 

*a*i  = *s/  + 4*7i/n4  4-  4(  2n  — l)*6J/(nV»)4-4(3n+  l)i53/(nan<»)  + (7n4-5)i44/(M3n<») 

+ 6*«ii/»3  + 4(5n  - 1)  4-  4(7n  4-  5)  ^/(nV)  4-  2(8m  4-  5)  *4M/(n*n(») 

4-  20(n  4- 1)  W(n*nW)  4-  4Asm/n2  4-  6(3n  + 1)  kmij(nn™)  4-  12(n  4- 1)  *„„/(nn«>) 

4-  28(n  + 1 )kamH*nf»)  4-  3(»  4-  l)A2M2/(nn‘»)  4-  Aluu/»  4-  8(n  + 1)  AS2111/n<» 

4-  8(n  4- 1)  *njn/n(S>4-  (n  4- 1)  *Mu  u/(»  - 1). 

*2  *i  = *8/»*  + 6An/n5  4-  16A82/n3  4-  26A53/ni4-  15/t^/n5  4-  15A811/n4  4-  66A521/n4  + 90  A4il/n4 
4-  60  A421/n4  4-  70AM2/n4  4-  20 *slll/w3  4-  105 A42u/ra3  4-  60A:33u/n3  4- 150  *3221/n3 
4- 1 oktnt/n3  4-  15A41111/ti24- 80A3211j/re2  4- 45A22211/ft2  4- 6A31lm/ra 
+ 15***iiii/to  + *11111111 

*?  = *s/*7  + 8*71/n*  4-  2Sknjna  4-  56AM>8  + 35A44/n8  + 28A811/w3  4-  168A521/«3 

+ 280A431/n3  4-210A4ni'n3  4-  280fc33a/»*  4-  56  Asm/n4  4-  420 A42U/«4  4-  280Amu/ti4 
+ 840A-3221/n4  4-  105fc2222/n4  4-  70  *41111/n3  4-  560 *321u/»3  4-  420Atmi/n3 
4-  66  A3ulll/n2  4-  210AMull/n2  4-  28A2uuu/»  4-  Aninm- 
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